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Chapter  1:  Titanium  Dioxide  Nanoparticles  and  Electrons  in  the  Solid  State 


1.1  Introduction 

The  properties  of  nano  particles,  aggregates  of  material  at  dimensions  less  than  100  nm  (1 
x  10'7  m),  have  been  of  much  recent  interest.  The  reduction  in  particle  size  to  the  nanometer 
scale  has  been  observed  to  significantly  change  the  properties  of  some  substances.  This  shift  in 
chemistry  as  a  function  of  size,  challenges  the  idea  of  a  material's  chemical  reactivity  as  an 
intensive,  i.e.  size  independent,  property. 

Gold  presents  an  example  of  a  size-mediated  shift  in  chemical  reactivity.  The  chemical 
inertness  of  gold  is  one  of  the  properties  that  contributed  to  its  high  valuation,  in  addition  to  its 
aesthetic  appeal.  For  most  of  human  history,  gold  has  been  perceived  as  eternal  and  unchanging. 
Although  it  is  not  completely  inert,  gold  in  bulk  form  does  not  easily  enter  into  chemical 
reactions.  Hammer  and  Norskov  described  gold  as  the  "noblest  of  all  metals",  and  explained  that 
this  low  reactivity  is  the  result  of  the  stability  of  gold's  electronic  structure.1  Gold  atoms  have 
filled  d-orbitals  and  relatively  high  ionization  potentials.  When  compared  to  the  other  coinage 
metals,  silver  and  copper,  gold's  first  valence  electron  is  more  difficult  to  remove  in  a  chemical 
reaction  (9.22  eV  vs.  7.58  and  7.73  eV  respectively).2  Gold  does  not  oxidize  under  ordinary 
environmental  conditions.  When  this  resistance  to  forming  oxide  coatings  is  considered  with 
gold's  relatively  low  electrical  resistivity  (2.255  x  10"8H  •  m  at  298  K),3  the  utility  of  gold  for 
protective  coatings  in  electrical  contact  surface  applications  is  obvious. 

Because  of  its  low  chemical  reactivity,  bulk  gold  was  regarded  as  a  poor  catalyst.  Bond 
described  bulk  gold  as  "the  least  catalytically  useful''.4  In  his  review  of  the  history  of  gold 
catalyst  research,  Haruta  observed  that  bulk  gold  was  regarded  as  almost  catalytically  inactive.5 
This  behavior  changes  significantly  when  gold  enters  reactions  in  particle  sizes  of  approximately 
5  nm  (5  x  109  m)  diameter.  On  this  size  scale,  gold  becomes  very  catalytically  active.  Gold 
nanoparticles  placed  upon  titanium  dioxide  supports  have  been  observed  to  dissociate  molecular 
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hydrogen  at  300  K6  and  catalyze  the  reduction  of  carbon  monoxide7.  The  application  of  catalyst 
systems  comprised  of  gold  and  metal  oxide  support  nanoparticles  to  the  detection  and 
decomposition  of  chemical  warfare  agents  (CWAs)  and  similar  materials  remains  an  active  area 
of  research.8"14 

Although  titanium  dioxide  nanoparticles  first  found  use  as  a  support  for  the  dispersal  of 
smaller  particles  of  noble  metals  such  as  gold  and  platinum,  the  titanium  dioxide  itself  also 
exhibits  unique  catalytic  properties  apart  from  the  supported  metal  particles.  The  application  of 
titania  nanoparticles  to  catalysis  is  appealing  due  to  its  relatively  low  material  cost  and  the 
numerous  routes  to  titanium  dioxide  nanoparticle  synthesis. 

Several  novel  applications  use  titanium  dioxide  nanoparticles  where  the  cost  of  equally 
effective  metal  catalysts  would  be  prohibitive.  With  alternative  approaches  to  supplying  the 
world's  growing  energy  needs,  titania  catalysts  combined  with  photoactive  dyes  have  attracted 
attention  for  use  in  photovoltaic  cell  applications.15"21  Photolysis  of  water  for  the  production  of 
hydrogen  fuel  is  another  energy  application  for  titanium  dioxide  nanoparticles.22'27  In  addition  to 
decomposition  of  chemical  warfare  agents,  titania  nanocatalysis  is  being  investigated  for 
treatment  of  environmental  pollutants.28'33  Titanium  dioxide  nanoparticle  catalysts  are  already 
commercially  produced  and  sold  under  the  brand  name  Aeroxide®  P-25  for  industrial  application 
by  the  Degussa  subsidiary  of  Evonik  Industries  A. G. 

Titanium  dioxide  nanoparticle  catalysis  is  dependent  upon  the  movement  of  electrons  to 
and  from  the  particle  surface,  where  adsorption  occurs.  The  factors  affecting  this  transport  of 
electrons  within  the  nanoparticle  derive  from  its  electrical  properties  and  the  extent  that  titanium 
dioxide  can  use  the  thermal  and  electromagnetic  energy  of  its  surroundings  to  cause  electric 
charge  separation.  The  most  easily  separated  electrons  within  titanium  dioxide  nanoparticles  are 
located  in  structural  features  known  as  shallow  trap  states.  These  sites  are  primarily  located  at 
defects  in  the  crystal  lattice  of  the  nanoparticle.  Within  this  thesis,  I  will  describe  an 
experimental  approach  that  applies  infrared  spectroscopy  to  the  measurement  of  the  energy 
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required  to  make  shallow  trap  state  electrons  available  for  oxidation-reduction  chemistry  with 
adsorbates  at  the  nanoparticle  surface.  That  energy  is  the  shallow  trap  state  energy. 
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e'  shallow  ti 
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Figure  1:  Semiconductor  band  structure  and  trap  states  within  a  nanoparticle.  Electrons  are 
represented  by  e  and  electron  holes  by  h+.  Adapted  from  Beydoun  etal.34 

1.2  E  lectrons  in  the  C  rystal line  Solid  State 
1.2.1  Probability  Distributions  for  Particle  E  nergies 

My  approach  to  calculating  the  energy  of  shallow  trap  states  in  crystalline  titanium 
dioxide  nanoparticles  requires  the  use  of  the  Fermi-Dirac  probability  distribution  for  electrons. 
Because  early  models  for  the  movement  of  electrons  within  solid  -  phase  matter  were  based  on 
energy  probability  distributions  derived  from  classical  Newtonian  mechanics,  I  will  review  the 
development  of  the  most  successful  classical  model,  the  M  axwell  -  Boltzmann  distribution. 

The  development  of  an  accurate  quantitative  description  of  the  interaction  of  energy  and 
matter  has  been  a  principle  goal  of  physical  science.  The  early  efforts  preceded  knowledge  of 
quantum  mechanics,  during  the  nineteenth  century.  Investigators  directed  their  work  at  a 
mathematical  model  that  could  calculate  the  energy  of  every  particle  comprising  a  closed  system. 
This  effort  began  with  an  assumption  that  Newtonian  mechanics  were  sufficient  to  derive  a 
comprehensive  model.  Because  of  the  large  number  (approximately  1020)  of  particles  in  even  a 
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miniscule  sample  of  solid  or  liquid  matter,  it  was  logical  to  first  examine  the  least  dense  state: 
gases. 

Over  the  course  of  the  nineteenth  century,  acceptance  of  the  atomic  theory  of  matter  was 
far  from  universal.  The  idea  that  matter  was  comprised  of  irreducible  particles  that  defined  its 
properties  is  known  to  have  existed  since  ancient  Greece  and  is  usually  attributed  to  Democritus 
(460  -  370  B.C.).35  The  nature  of  atoms  themselves  remained  a  mystery  until  the  early  twentieth 
century  and  in  the  absence  of  more  evidence  many  physicists  were  skeptical  as  to  their  very 
existence.  As  late  as  1897,  Ernst  M  ach  openly  declared  his  disbelief  in  the  existence  of  atoms 
during  a  presentation  before  the  Viennese  Academy  of  Sciences.36  Physicists  used  the  term 
'particle'  in  reference  to  atoms  and/or  molecules.  This  substitution  was  motivated  by  their 
reluctance  to  fully  accept  atomic  theory  and  avoidance  of  tedious  repetition  of  the  phrase  'atoms 
and/or  molecules'  in  settings  where  either  might  be  applicable.  I  continue  the  usage  of 
'particles'  in  this  section  for  the  latter  reason. 

The  starting  point  for  consideration  of  the  behavior  of  gaseous  matter  is  the  kinetic  theory 
of  gases.  This  model  is  comprised  of  several  axioms.  Included  are  that  gases  are  comprised  of 
particles  whose  sizes  are  greatly  smaller  than  their  volume  of  confinement.  The  disparity 
between  particle  size  and  sample  size  is  great  enough  to  render  the  sizes  of  the  particles 
insignificant  and  their  influences  to  those  of  point  entities.  Another  aspect  of  kinetic  theory  is 
that  confined  gas 'particles'  travel  at  high  rates  of  speed  and  make  collisions  with  other  particles 
and  the  walls  of  their  container.  The  result  of  these  collisions  is  that  kinetic  energy  exchanges 
between  the  colliding  entities.  If  two  particles  collide,  they  leave  the  collision  with  different 
velocities.  If  a  particle  collides  with  a  container  wall,  the  result  is  a  miniscule  movement  of  the 
wall.  Gas  pressure  is  the  collective  result  of  numerous  instances  of  particle-wall  collision. 

Another  assumption  of  the  kinetic  theory  of  gases  is  that  for  a  sample  of  gas  in  thermal 
equilibrium,  the  gas  particles  will  exist  with  a  non-uniform  distribution  of  speeds.  This 
statement  is  equivalent  to  saying  that  the  particles  will  exist  in  thermal  equilibrium  with  an 
unequal  distribution  of  kinetic  energies.  This  partitioning  of  kinetic  energies  in  gases  is  evident 
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in  the  Earth's  atmosphere,  and  is  described  in  an  example  given  by  Gilbert  Castellan  in  his  text 
Physical  Chemistry.37  The  argument  is  that  if  all  atmospheric  gas  particles  possessed  the  same 
upward  component  of  velocity,  then  at  ground  level  all  particles  would  be  limited  to  a  maximum 
altitude  which  they  could  reach  under  the  influence  of  gravity.  This  situation  would  lead  to  the 
Earth's  atmosphere  possessing  a  sharply  defined  boundary,  beyond  which  no  gas  molecules 
could  travel.  In  addition,  the  density  of  the  atmosphere  would  increase  with  altitude  as  particles 
near  the  end  of  their  upward  travel  would  slow  down  and  accumulate  at  high  altitude.  Neither  of 
these  two  conditions  are  observed  in  the  Earth's  atmosphere.  Atmospheric  density  decreases 
with  altitude  and  this  property  is  directly  observable  during  an  ascent  up  a  mountain.  Another 
consideration  is  that  the  atmosphere  does  not  have  a  well-defined  boundary,  as  is  demonstrated 
by  the  existence  of  significant  atmospheric  drag  effects  upon  orbiting  spacecraft.38,  39  The 
conclusion  is  that  thermally  equilibrated  gas  particles  exist  with  an  unequal  distribution  of 
kinetic  energies. 

Several  researchers  had  addressed  the  question  of  how  confined  gas  particles  distribute  a 
fixed  quantity  of  energy.  Between  1738  and  1845,  Daniel  Bernoulli,  John  James  Waterston,  and 
John  Herapath  each  proposed  models  for  gas  behavior  that  anticipated  parts  of  the  kinetic  theory 
of  gases,  although  their  efforts  were  either  forgotten  or  ignored  by  the  scientific  establishment  of 
the  nineteenth  century.36 

The  Scottish  physicist  James  Clerk  Maxwell  published  a  theory  for  the  distribution  of 
kinetic  energy  within  a  sample  of  confined  gas  in  1860  and  followed  this  paper  with  another  in 
1867  where  he  offered  more  justifications  for  his  approach.36  The  derivation  of  the  M  axwell 
distribution  is  based  upon  Newtonian  descriptions  of  the  component  velocities  of  a  particle  in 
each  of  the  Cartesian  axes  of  space.  M  axwell 's  equation  for  the  distribution  of  number  nc  of  gas 
particles  with  speed  between  c  and  c  +  dc  is  given  by  the  equation: 
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where  N  is  the  total  number  of  particles,  m  is  the  particle  mass,  NA  is  Avogadro's  number  (6.022 
x  1023  particles  •  mol"1),  R  is  the  gas  constant  (8.313  J  •  mol"1  •  K )  and  T  is  the  thermodynamic 
temperature  of  the  confined  gas  system.  Recognition  that  the  ratio  R/NA  is  equivalent  to 
Boltzmann's  constant  kB  (1.381  x  10"23J  •  K"1)  and  that  E kinetic  =  Vimc2  allows  Equation  (1)  to  be 
rewritten  as: 
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The  M  axwell  distribution  is  frequently  rearranged  for  convenience  to  the  probability  distribution: 
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The  terms  on  the  left  side  of  the  equation  represent  a  probability  that  a  particle  of  mass  m  will 
have  speed  c  when  thermally  equilibrated  to  temperature  T.  This  probability  function  has  the 
units  of  sec  •  m"1  and  is  represented  by  the  function  f(c): 
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The  use  of  the  term  'distribution'  in  the  Maxwell  distribution  is  a  reference  to 
classification  of  particles  into  categories  or  intervals.  In  this  case,  the  intervals  are  particle 
speeds  and  the  width  of  those  categories  is  taken  to  be  the  infinitesimal  increment  dc.  While  it 
may  be  implicit  that  an  expression  for  particle  velocities  based  upon  Newtonian  mechanics  could 
be  used  to  calculate  the  velocity  and  kinetic  energy  for  every  constituent  particle  of  a  system, 
such  comprehensive  characterizations  of  a  system  are  impractical.  As  an  alternative,  statistical 
characteristics  for  a  gas  particle  system  are  derived  from  the  M  axwell  distribution. 
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The  mean  speed,  cavg,  for  a  system  can  be  derived  from  Equation  (4)  and  is  given  by  the 
expression: 


C 

avg 
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where  M  is  the  molar  mass  of  the  particle.  Another  useful  statistic  that  can  be  derived  from 
Equation  (4)  is  the  root  mean  square  speed,  crms- 
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The  most  probable  speed,  c*,  is  derived  from  Equation  (4)  by  evaluating  the  first  derivative  of 
f(v)  with  respect  to  c,  and  finding  the  speed  at  which  the  derivative  is  zero  (excepting  the  cases  c 
=  0  and  c  =  oo).  The  resulting  expression  for  c*  is: 
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The  relative  mean  speed,  crei,  is  the  speed  at  which  two  particles  approach  each  other.  This 
speed  is  derived  from  the  mean  speed  (Equation  (5))  and  for  two  particles  of  equal  mass  it  is 
given  by: 
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For  two  particles  of  unequal  mass  the  crei  is  given  by: 
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where  k  and  T  have  their  previously  describe  values  and  the  term  //  is  the  reduced  mass  of  the 
two  particles,  given  by  the  equation: 
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The  terms  ma  and  mb  are  the  respective  masses  of  the  two  particles  being  described.  The 
frequency  of  particle  collisions,  z,  is  given  by  the  equation: 


z  = 


(7C  ,N 

re I 

V 


(11) 


where  N  is  the  number  of  particles  and  V  is  the  volume  of  their  confinement.  The  term  a  is  the 
collision  cross  section  of  the  particle,  given  by: 


a  =Kd 2  (12) 

where  d  is  the  particle  diameter.  All  of  these  derivations  from  Equation  (4)  lead  to  an  expression 
for  the  average  distance  traveled  by  a  particle  between  collisions,  the  mean  free  path  1: 

(13) 

z 

M  ean  free  path  descriptions  will  become  important  in  the  development  of  models  for  conduction 
of  electrons  in  metals.  An  example  of  the  M  axwell  speed  distribution  for  nitrogen  gas  atT  =  298 
K  and  T  =1298  K,  with  c,  cavg,  and  crms  identified  at  each  temperature  is  given  in  Figure  2. 
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Figure  2:  The  Maxwell  distribution  for  nitrogen  atT  =  298  K  and  T  =3296  K. 


Ludwig  Boltzmann  published  a  paper  in  1868  that  provided  additional  justifications  for 
the  Maxwell  distribution  and  extended  it  to  include  not  just  kinetic  energy  but  other  forms  as 
well  (e.g.  potential  energy).  Because  of  Boltzmann's  work  in  validating  Maxwell's  previous 
work  and  his  extension  of  the  theory,  the  energy  partition  probability  is  referred  to  as  the 
M  axwell  -  Boltzmann  distribution  and  it  is  given  by  the  equation: 
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The  terms  of  the  Equation  (14)  have  their  previously  described  meanings  and  values.  Potential 
energy  is  described  as  the  energy  of  configuration  for  the  particles.40  An  example  of  potential 
energy  is  given  by  considering  the  upward  trajectory  of  a  particle  in  a  gravity  field.  At  the  start 
of  the  particle  motion  at  ground  level,  all  of  the  particle  energy  is  kinetic.  During  the  upward 
travel  of  the  particle,  the  energy  is  kinetic  and  potential.  At  the  maximum  height  of  travel,  when 
the  particle  upward  motion  stops  and  speed  is  equal  to  zero,  all  of  the  particle  energy  is  potential. 
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During  downward  travel  the  particle  again  possesses  a  combination  of  kinetic  and  potential 
energies,  and  at  impact  upon  the  ground  all  energy  is  once  again  kinetic. 

L22  Classical  Free  Electron  Models 

Researchers  anticipated  the  existence  of  a  fundamental  charge  carrier  for  electric  current 
for  centuries.  The  ability  of  objects  to  accumulate  static  electric  charges  implied  that  something 
was  present  on  their  surfaces  which  carried  the  electric  charge  and  that  entity  could  be 
transferred  off  the  charged  object  upon  discharge.  With  the  discovery  of  the  electron  in  1897, 
John  Joseph  Thomson  established  the  particle  nature  of  electric  charge.41 

Thomson's  experiments  leading  to  the  discovery  of  the  electron  were  based  upon  the 
deflection  of  cathode  rays  by  electric  and  magnetic  fields.  Cathode  rays  are  actually  a  stream  of 
electrons.  One  means  of  imparting  enough  energy  to  the  electrons  of  a  sample  of  metal  is  to  heat 
it.  I  n  the  presence  of  a  reactive  gas  such  as  oxygen,  the  metal  is  oxidized  by  the  loss  of  electrons 
and  the  reacting  gas  is  reduced  through  the  gain  of  those  departing  electrons.  In  a  vacuum,  when 
the  heated  metal  has  an  electric  potential  applied  to  it  with  respect  to  an  anode,  an  electric  current 
flows  from  the  cathode  to  the  anode.  A  schematic  of  the  Thomson  cathode  ray  deflection 
experiment  is  given  in  Figure  3.  This  effect  is  the  basis  for  operation  of  thermionic  vacuum 
tubes  and  cathode  ray  tube  displays  that  were  the  basis  for  most  modern  electronic  devices 
during  the  twentieth  century. 
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Figure  3:  The  Thomson  cathode  ray  deflection  experiment.  The  cathode  is  the  heated  filament  F  and 
the  anode  is  A.  E  indicates  lines  of  electric  force  and  B  indicates  lines  of  magnetic  force,  perpendicular 
to  the  page.  HV  is  a  high  voltage  power  supply  and  S  is  a  phosphorescent  screen.40 


In  the  course  of  his  experiments  leading  to  the  discovery  of  the  electron,  Thomson  noted 
that  the  ability  to  emit  electrons  was  not  particular  to  any  specific  metal  provided  sufficient 
energy  was  supplied.  In  1900,  Paul  Drude,  A  contemporary  of  Thomson,  proposed  that  electrons 
might  be  a  component  of  the  metal  atoms  themselves  rather  than  being  solely  supplied  by  the 
circuit  providing  electric  potential  to  the  cathode  ray  experiment. 

Drude  proposed  a  model  for  the  structure  of  metals  at  the  atomic  scale  that  accounted  for 
some  observed  behaviors.42  The  Drude  model  for  metals  was  based  upon  several  assumptions 
that  include: 

•  Metals  possess  a  population  of  delocalized,  mobile  electrons  that  behave  like  a  sample  of 
confined  gas  molecules. 

•  The  population  of  mobile  electrons  is  formed  by  the  donation  on  one  electron  by  each  metal 
atom,  leaving  behind  a  stationary,  positively  charged  atom  core. 

•  The  delocalized,  mobile  electrons  are  responsible  for  the  metal's  ability  to  conduct  an  electric 
current. 


11 


•  The  delocalized  conduction  electrons  move  with  a  M  axwell  -  Boltzmann  root  mean  square 
speed  (Equation  (6)). 

Drude's  assumption  that  mobile  electrons  in  a  metal  behave  like  a  gas  conforming  to  M  axwell  - 
Boltzmann  probability  distribution  requires  that  the  energies  of  the  electrons  are  totally  kinetic 
and  the  electron  gas  as  a  whole  has  a  total  energy  E,  given  by: 

E  =—nkT  (15) 

2  B 

where  n  is  the  number  of  conduction  electrons  present  in  the  bulk  sample  of  metal.  The  classical 
mechanics  derived  heat  capacity  of  the  electron  gas,  Cci  would  be  given  by  the  expression: 

C,  =^nkB  (16) 

One  implication  of  Equation  (16)  is  that  metals  should  have  a  heat  capacity  almost  double  of 
what  is  actually  observed.  Another  fallacy  of  Equation  (16)  is  that  the  heat  capacity  of  the 
conduction  electron  'gas'  is  independent  of  temperature,  where  the  observed  trend  is  that  the 
specific  heat  is  linearly  dependent  upon  temperature. 

The  Drude  model  results  in  a  derived  thermal  electrical  conductivity  a  for  metals  given 
by  the  equation: 


J3m  k  T 

where  e  is  the  elementary  charge  (1.6022  x  10"19  C),  2  is  the  electron  mean  free  path  (Equation 
(13)),  and  me  is  the  rest  mass  of  the  electron  (9.1094  x  10'31  kg).  Equation  (17)  requires  that 
metal  electrical  conductivity  vary  inversely  with  T1/2,  which  is  contradictory  to  the  observed  T  1 
dependency  of  electrical  conductivity.  D  rude  argued  unsuccessfully  that  the  mean  free  path  2  for 
conduction  electron-positive  atom  core  collisions  must  have  a  T  1/2  dependency. 
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In  1905  H.  A.  Lorenz  suggested  a  modified  form  of  the  Drude  free  electron  model.43 
Lorenz  considered  that  the  Maxwell  distribution  of  the  conduction  electrons  would  have  its 
symmetry  disturbed  by  the  effects  of  either  an  applied  electric  field  or  a  thermal  gradient.  The 
Lorenz  model  led  to  a  derived  electric  conductivity  a  for  metals  of: 


4nen  A 


3  J 27lm  kKT 

V  e  B 


(18) 


The  Lorenz  conductivity  leads  to  a  result  that  is  0.92  of  the  Drude  conductivity.  By  considering 
the  perturbation  of  conduction  electron  scattering  by  an  external  field,  the  Lorenz  model  predicts 
an  expression  metal  conductivity  in  a  magnetic  field  of: 


JxB  |2 

— *  — * 

where  J  is  the  magnitude  of  the  conduction  electron  path  vector  J  and  B  is  the  magnetic  field 
vector.  The  term  RH  is  the  Hall  coefficient,  which  represents  the  magnitude  of  the  transverse 
electric  field  in  a  conductor  that  is  subjected  to  a  magnetic  field.  The  Lorenz  model  leads  to  a 
Hall  coefficient  of: 


J2  +  (<jRh)2 


R 


H 


3n 
8  ne 


(20) 


The  Drude  model  predicted  that  all  electrons  would  have  the  same  thermal  velocity  and  this 
assumption  would  result  in  no  magnetoresistive  effect.  The  Lorenz  model  predicted  a 
disturbance  of  M  axwell  -  Boltzmann  distribution  symmetry  that  would  result  in  electrons  being 
subjected  to  forces  that  cause  their  departure  from  straight  conduction  paths  under  the  influence 
of  a  magnetic  field.  This  prediction  of  magneto  resistance  in  metals  was  an  improvement  over 
the  Drude  model. 
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The  inability  of  free  electron  conduction  models  for  metals  to  accurately  predict  electrical 
conductivity  or  specific  heat  led  to  continuing  efforts  to  improve  them.  One  large  failing  of  both 
the  Drude  and  Lorenz  models  was  their  neglect  to  address  the  magnetic  susceptibility  of  metals. 
Pauli  exclusion  spin  pairing  of  electrons  was  unknown  when  Drude  and  Lorenz  proposed  their 
metal  conductivity  models.  One  would  expect  that  the  large  population  of  mobile  charge  carriers 
that  free  electron  models  assume  to  be  present  would  lead  to  strong  magnetic  susceptibilities  in 
all  metals.  That  conclusion  is  not  consistent  with  the  actual  observed  properties  of  metals,  which 
vary  considerably  in  magnetic  susceptibilities. 

LZ3  Quantized  Free  Electron  Model 

Working  independently  in  the  early  twentieth  century,  Enrico  Fermi  and  Paul  Adrien 
Maurice  Dirac  (P.A.M.  Dirac)  derived  an  expression  for  the  probability  that  an  electron  will 
occupy  an  energy  state  within  a  substance.  That  probability  p(E)  is  given  by  the  equation: 


p(E)  = 


1 


f  E-E  ^ 
1  +  Exp  F 


V  V  , 


(21) 


The  term  E  is  the  energy  of  the  electron.  The  term  EF  is  the  Fermi  energy  of  the  material  under 
consideration.  The  Fermi  energy  is  defined  as  the  energy  of  the  highest  occupied  state  for  a 
material  at  absolute  zero  temperature  (T  =  0  K).  Figure  4  shows  the  energy  increase  with  the 
increasing  number  of  available  electron  states,  ns(E)  for  a  system  atT  =  0  K.  By  definition,  no 
states  above  the  Fermi  Energy  EF  can  be  filled  atT  =  0  K. 
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Figure  4:  Filling  of  electron  states  at  the  Fermi  energy  (E  F ). 40 


In  1927,  Arnold  Sommerfeld  applied  the  F ermi-D i rac  probability  function  to  the  free 
electron  model  for  metals.42  The  Sommerfeld  quantized  free  electron  model  is  based  on  the 
observation  that  only  the  electrons  in  a  metal  that  are  energetically  close  to  the  Fermi  energy 
have  sufficient  freedom  to  be  scattered  to  other  energy  levels.  Electrons  occupying  states  at 
energies  significantly  lower  than  the  Fermi  energy  are  effectively  immobilized  due  to  the 
absence  of  available  empty  states  in  their  vicinity.  The  effect  of  this  modification  on  the  Drude 
and  Lorenz  free  electron  models  was  to  greatly  reduce  the  number  of  electrons  available  to 
participate  in  the  'electron  gas'  that  is  a  key  component  of  free  electron  models. 

With  fewer  electrons  contributing  to  the  bulk  properties  of  the  metal,  the  Sommerfeld 
model  was  able  to  more  accurately  estimate  the  electronic  component  of  a  metal's  specific  heat: 


rm2k 2 

C  = - 2. T 

e  eE, 


(22) 


The  Sommerfeld  model  predicted  a  magnetic  susceptibility  x  for  a  non-degenerate  electron  gas 
given  by  the  equation: 


15 


(23) 


nA  1 

3 *bT 


The  term  gB  is  the  Bohr  magnetron  (9.727  x  10"24  J  T"1).  The  most  notable  feature  of  Equation 
(23)  is  its  inverse  dependence  upon  temperature  and  the  absence  of  a  magnetic  field  dependency. 
This  type  of  behavior  is  not  observed  for  metals.  Wolfgang  Pauli  derived  an  expression  for  the 
total  magnetic  susceptibility  xtotai  that  removed  the  temperature  dependency  and  was  successful 
in  predicting  the  behavior  of  some  alkali  metals: 


X 


total 


(24) 


Sommerfeld's  model  resulted  in  an  expression  for  electrical  conductivity  given  by: 

_  ne~X(EF) 
ms(EF) 


(25) 


where  X(EF)  is  the  mean  free  path  length  of  electrons  at  the  Fermi  energy  and  s(EF)  is  the  speed 
of  electrons  at  the  Fermi  energy.  Equation  (25)  has  no  temperature  dependency,  which 
contradicts  the  observed  temperature  dependent  conductivity  behavior  of  metals. 

The  Sommerfeld  free  electron  model  for  the  behavior  of  metals  was  able  to  predict  the 
heat  capacities  and  magnetic  susceptibilities  with  some  limited  accuracy.  The  application  of 
Fermi-Dirac  statistics  to  limit  the  number  of  electrons  participating  in  the  'electron  gas'  of  a 
metal  was  anticipatory  of  the  presently  accepted  model  for  solid-state  electrical  behavior,  band 
theory.  Flowever,  the  inability  of  the  Sommerfeld  model  to  account  for  the  observed  temperature 
dependence  of  electrical  conductivity  in  metals  cast  doubt  upon  whether  it  was  an  accurate 
depiction  for  electron  transport  processes  in  metals. 
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LZ4The  Band  Theory  Model 

By  the  late  1920's,  Felix  Bloch,  a  graduate  student  of  Werner  Heisenberg,  began  work  on 
the  shortcomings  of  the  Lorentz  free  electron  gas  model  as  the  topic  of  his  thesis.  While  the 
Sommerfeld  model  for  electron  behavior  in  metals  had  introduced  quantum  mechanics  rather 
than  the  classically  derived  Boltzmann  probability  distribution,  Sommerfeld  clung  to  the  idea 
that  electrons  in  metals  moved  freely  as  do  gas  molecules  in  a  container.  Bloch  was  skeptical  of 
this  point  of  the  Sommerfeld  model.44 

The  Bloch  model  for  behavior  of  electrons  in  metals  shortly  followed  the  publication  of 
Erwin  Schrodinger's  wave  equation  for  the  wave-particle  duality  of  matter.  The  time- 
independent  form  of  the  wave  equation  is: 


tr 

- V2  +(V  (?) 

2  m 


\\f  =  c«5\\f 


(26) 


The  term  h  is  Planck's  constant  divided  by  2 n  (1.055  x  10"34  J  •  s),  m  is  the  particle  mass,  and 
fy^Cr)  is  the  particle  potential  energy  function.  The  term  y/  is  the  wavefunction  for  the  particle 
described  by  the  equation,  and  \72is  the  second  order  differential  operator.  The  significance  of 
the  wavefunction  y/  is  that  the  square  of  the  wavefunction  (or  its  square  modulus,  if  y/  is 
complex)  is  proportional  to  the  probability  of  locating  the  electron  at  position  r.  Generally,  the 
process  of  solving  a  Schrodinger  wave  equation  is  one  of  determining  the  correct  functions  for  y/ 
and  (V(~r).  The  first  and  leftmost  set  of  terms  in  the  brackets  of  Equation  (26)  represent  the 
kinetic  energy  of  the  particle  under  consideration  and  ^  represents  its  total  energy.  The  vector 
T  describes  the  Cartesian  coordinate  position  of  the  electron  within  the  crystal  lattice. 

Bloch  sought  to  apply  the  Schrodinger  equation  to  the  description  of  electrons  within 
metals.  The  Bloch  model  for  metals  as  described  by  Blakemore43  was  based  upon  several 
assumptions: 

•  The  Schrodinger  wave  function  y/  is  calculated  for  a  perfect  periodic  lattice,  with  no 
vibrations  or  defects. 
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•  The  Bloch  model  is  a  description  for  a  single  electron.  This  electron  'sees'  an  effective 
potential  of  v(~r)  for  everything  else  in  the  lattice.  This  effective  potential  is  not  identical  to 
the  Schrodi nger  potential  energy  (V'fr ). 

•  The  electron  is  described  by  a  one-electron  Schrodi  nger  equation  (previously  given  as 
Equation  (26) ). 

Bloch  reasoned  that  the  potential  v(~r)  experienced  by  the  electron  was  the  sum  of  two 
components. 

The  first  part  of  v(?)  is  the  electrostatic  potential  of  the  positive  ion  cores  that  remain 
fixed  in  the  crystal  lattice  following  loss  of  an  electron  to  the  delocalized  population  of  electrons. 
This  potential  is  periodic,  repeating  with  the  frequency  of  the  unit  cells  in  the  metal  crystal  lattice 
structure.  The  second  part  of  vCr)  is  the  electrostatic  potential  of  all  other  outer  shell  electrons. 
Like  the  positive  ion  core  potential,  this  force  is  also  presumed  to  repeat  with  the  periodicity  of 
the  unit  cells  of  the  lattice.  The  total  potential  v(~r)  is  related  to  the  Schrodi  nger  potential  energy 
function  (¥'(?)  by  the  equation: 

(V  (r)  =-ev(r)  (27) 

where  eis  the  elementary  charge  (1.6022  x  10  19  C).  The  Bloch  solution  for  the  lattice  electron 
wavefunction  yj  is: 


y/k(r)  =  Uk(r)exp(ik»r)  (28) 

___  — ^  ( 

The  term  Ac  is  a  wave  vector  that  describes  the  direction  and  periodicity  of  the  wave  associated 

__  — ^ 

with  the  electron.  The  wave  vector  k  is  related  to  electron  momentum  described  by  the  equation 

k  =£  (29) 

ti 

.  —* 
where  p  is  the  momentum  vector  for  the  electron.  The  Cartesian  scalar  components  of  k  are  kx, 

ky,  and  kz.  These  scalar  components  comprise  a  spatial  system  known  as  k-space.  The  term  /  in 


18 


Equation  (28)  is  imaginary  root  of  negative  one.  Lastly,  the  function  UkCr)  is  a  function  with 

the  periodicity  of  the  crystal  lattice,  whose  form  is  dependent  upon  the  value  of  r. 

The  solution  of  the  Bloch  wavefunction  given  in  Equation  (28)  describes  a  partition  of 

the  electrons  within  a  bulk  solid  into  zones  that  are  referred  to  as  bands.  Imaginary  values  of  the 
— ^ 

wave  vector  k  lead  to  real  values  for  the  exponential  term  of  Equation  (28)  and  the  equation 

_ _  — A 

diverges.  This  circumstance  describes  a  forbidden  zone  for  electrons.  The  converse,  when  Ac  is 
real,  leads  to  an  imaginary  value  for  the  exponential  term.  Under  this  circumstance,  Euler's 
formula  is  applicable  and  Equation  (28)  becomes: 

y/k(r)  =  £/,  (r)[cos(£*r)  +  /sin(A*r)j  (30) 

Equation  (30)  describes  a  wave  that  is  capable  of  infinite  propagation  in  a  perfect  crystal  lattice. 
This  equation  describes  an  allowed  band  for  the  electron. 

The  partitioning  of  electrons  in  a  metal  can  be  demonstrated  by  examination  of  a  one¬ 
dimensional  lattice.  A  one  dimensional  lattice  model  would  consist  of  a  large  number  of  ions, 
spaced  periodically  and  arranged  in  a  line  as  shown  in  Figure  5: 


The  ion  elements  of  the  one-dimensional  lattice  acquired  a  positive  charge  through  donation  of 
an  outer  shell  electron  to  the  delocalized  conduction  electron  population 
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From  the  perspective  of  a  conduction  electron,  each  ion  is  the  center  of  an  electrostatic 
potential  field.  In  order  for  an  electron  to  'escape'  from  an  ion  potential  well,  either  it  must 
acquire  energy  greater  than  r%"or  it  must  pass  through  the  energy  barrier  separating  it  from  an 
adjacent  ion  by  a  quantum  tunneling  mechanism. 

The  calculation  of  a  wavefunction  for  a  one-dimensional  lattice  such  as  Figure  5  can  be 
greatly  simplified  by  assuming  the  wells  are  rectangular.  Kronig  and  Penny  described  Bloch 
wavefunction  calculations  for  a  series  of  rectangular  potential  energy  wells  in  1931. 43, 45, 46 
While  the  Kronig-Penny  model  surely  is  an  oversimplification,  it  does  demonstrate  the  evolution 
of  band  and  band  gap  structures  for  delocalized  electrons  in  crystal  lattice  structures.  While  the 
example  that  follows  is  limited  to  a  one-dimensional  array,  these  techniques  can  be  extended  to 
three  dimensions  for  the  description  of  electron  behavior  in  solids. 

The  Kronig-Penny  model  consists  of  rectangular  wells  of  width  a,  separated  by 
rectangular  barriers  of  width  b  and  potential  energy  height  <5%T  The  arrangement  of  wells  and 
barriers  is  periodic  and  while  Figure  6  only  shows  four  wells  and  barriers,  a  larger  number  is 
assumed  to  comprise  the  one-dimensional  array. 


barrier 


well 

a  +  b  2(  a  +  b)  3(a  +  b) 

Figure  6:  The  Kronig-Penny  model. 


4(a  +  b) 
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Like  the  ion  lattice  of  Figure  5,  the  available  pathways  for  an  electron  to  leave  a  well  in  Figure  6 
are  either  to  exceed  the  well  potential  energy  or  to  tunnel  through  a  barrier. 

If  the  barriers  in  the  Kronig-Penny  model  are  thin  enough  to  allow  unimpeded  passage  of 
electrons  between  wells,  the  energy  £  of  an  electron  is  given  by  the  dispersion  relationship: 


E  = 


n2k 2 


2  m 

e 


(31) 


The  terms  h  and  me  have  their  previously  described  values  and  the  term  k  is  the  magnitude  of  the 

— } 

electron  wave  vector  k.  At  the  other  extreme,  when  the  barriers  of  Figure  6  are  too  thick  to 
permit  any  quantum  mechanical  tunneling,  the  energy  of  an  electron  trapped  in  a  well  is  given 
by: 

nn2tr 

E  = - —  where  n  =1,2,3,---  (32) 

2  ma 

— ^ 

For  the  intermediate  case,  where  barriers  allow  only  some  tunneling,  the  electron  wave  vector  k 
and  electron  energy  are  related  by  the  equation: 


cos  (ka)  =  cos  (aa)  +  P 


sin(aa) 

aa 


(33) 


The  term  a  is  related  to  electron  energy  by: 


2m  E 

e 

ti1 


1/2 


(34) 


The  term  P  in  Equation  (33)  is  a  descriptor  for  the  resistance  of  the  Figure  6  barriers  to  electron 
tunneling: 


P  = 


abm  V 


n2 


(35) 
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If  the  term  P  in  Equation  (33)  is  set  at  a  finite  value,  the  equation  will  alternate  between  real  and 

imaginary  values  of  the  wave  vector  magnitude  k.  Values  for  cos(/ca)  greater  than  lor  less  than 

— > 

-1  correspond  to  imaginary  k values,  and  the  wave  vector  k  describes  a  forbidden  energy  state 

for  the  electron.  Alternately,  when  the  magnitude  of  cos [ka)  is  less  than  or  equal  to  1,  Ar is  real 

— ^ 

and  the  corresponding  wave  vector  k  is  associated  with  an  allowed  energy  state  for  the  electron. 
The  regions  of  allowed  energy  states  for  the  electron  are  analogous  to  electron  bands,  and  the 
forbidden  zones  are  the  band  gaps: 


If  the  values  of  cos(fec)  in  Figure  7  that  fall  between  1  and  -1  are  solved  for  /cand  plotted  against 
their  corresponding  energy  values  from  a,  the  band  structure  for  the  Kronig-Penny  model 
becomes  apparent.  A  plot  for  an  unrestricted  electron  (P  =  0)  is  indicated  by  a  broken  line  for 
comparison: 
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Figure  8:  Electron  band  structure  for  the  Kronig-Penny  model  for  a  value  of  P  =  2.  The  broken  line 
represents  electron  energies  for  a  free,  unrestricted  electron  (P  =0). 

An  alternate  description  of  the  formation  of  band  structure  arises  from  consideration  of 
atomic  and  molecular  orbitals  in  the  solid  state.  The  tight-binding  approximation  used  in  solid- 
state  physics  begins  with  the  assumption  that  a  macroscopic  sample  of  solid  phase  material 
behaves  as  a  single  macromolecule.  The  overlap  of  atomic  orbitals  within  heteronuclear  solid- 
state  materials  leads  to  the  formation  of  discrete  molecular  orbitals.  As  the  number  of  atoms  in  a 
bulk  solid  phase  sample  increases  (  >  104  atoms  )47,  atomic  and  molecular  orbitals  become 
confined  within  a  physical  space  comparable  to  the  spatial  dimensions  of  their  associated 
wavefunctions.  The  result  is  that  the  orbitals  overlap  and  discrete  orbital  structure  merges  into 
delocalized  band  structure  for  the  whole  solid.48  Just  as  no  two  electrons  within  an  atom  may 
have  the  same  set  of  quantum  numbers,  no  two  electrons  existing  within  the  band  structures  for 
the  bulk  crystalline  sample  may  have  exactly  the  same  energy  state.  This  energy  state 
exclusivity  within  the  band  structures  for  electrons  in  crystal  lattices  causes  a  Boltzmann-like 
distribution  of  energy  states  within  the  bands. 
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LZ5  Band  Structure  Classification  of  Solid  Phase  Matter 

Solid  phase  matter  can  be  classified  based  on  its  electrical  conductivity.  A  material 
whose  electric  conductivity  varies  inversely  with  temperature  is  a  conductor.  Solid  phase 
conductors  are  characterized  by  partially  filled  valence  bands  that  allow  electron  mobility,  as 
shown  in  Figure  9a.  As  the  temperature  of  a  conductor  is  increased,  lattice  vibrations  cause  an 
increase  of  collisions  between  charge  carriers  and  the  lattice  ion  cores.  The  result  of  this 
increasing  frequency  of  carrier-ion  collision  is  that  the  charge  carriers  are  scattered  from  their 
paths  under  the  influence  of  an  externally  applied  field  and  conductivity  decreases. 

A  semiconductor  is  a  material  characterized  by  a  directly  proportional  increase  in  electric 
conductivity  with  increasing  temperature.  Semiconductors  are  materials  whose  valence  bands 
are  completely  filled,  with  conduction  bands  sufficiently  close  for  excitation  of  electrons  to  the 
conduction  band  to  be  feasible.  Figure  9b  and  c  give  the  corresponding  diagrams  for  a 
semiconductor  before  and  after  electron  excitation  to  the  conduction  band,  respectively. 
Following  excitation,  two  modes  of  electric  current  conduction  become  available.  The  electron 
excited  to  the  sparsely  populated  conduction  band  has  freedom  of  movement,  and  the  vacancy 
created  by  the  departed  electron  in  the  valence  band  provides  a  'hole'  for  adjacent  electrons  to 
move  to.  The  overall  effect  for  the  excited  semiconductor  in  Figure  9c  is  that  electric  current  is 
carried  by  both  the  valence  and  the  conduction  band.  Because  thermal  excitation  of  electrons 
into  the  conduction  band  is  a  possibility  for  some  narrow  band  gap  semiconductors,  the 
conductivity  of  semiconductors  as  a  class  of  materials  is  described  as  increasing  directly 
proportional  with  temperature. 

A  third  possibility  for  solid  phase  matter  is  shown  in  Figure  9c.  FI  ere,  the  valence  band  is 
filled  as  was  the  case  with  a  semiconductor,  but  the  band  gap  is  too  large  for  significant 
excitation  of  electrons  to  the  conduction  band  to  occur.  These  materials  are  referred  to  as 
electric  insulators,  although  they  are  actually  large  bandgap  semiconductors.  Crystalline  sodium 
chloride  is  an  example  of  an  insulator,  with  a  band  gap  of  approximately  7  eV.49 
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Figure  9:  Classification  of  solid  phase  materials  from  band  structure.  VB  is  valence  band  and  CB  is 
conduction  band. 


The  diagrams  of  solid-state  band  structure  shown  in  Figure  9  are  simplifications.  A  more 

accurate  diagram  of  solid-state  band  structure  is  the  plot  of  electron  energy  versus  electron  wave 
— ^  — * 

vector  (i.e.  k).  The  proportionality  of  k  to  electron  momentum  was  previously  shown  in  the 

— # 

wave  vector  definition  given  as  Equation  (29).  Because  of  this  proportionality,  k  vector  space  is 

— ^  t 

alternately  referred  to  as  k  momentum  space.  An  example  of  semiconductor  band  structure  in  k 
momentum  space  is  shown  in  Figure  10.  The  parabolic  shape  of  both  the  valence  and 
conduction  bands  is  a  result  of  the  expression  for  electron  kinetic  energy: 


(36) 


where p  is  the  electron  momentum  and  me  is  the  electron  mass.  The  downward  concavity  of  the 
valence  band  parabola  symbolizes  the  decrease  in  energy  an  electron  experiences  when  it  is 

— t 

accelerated  from  maximum  energy  at  k  =  0  in  a  partially  filled  valence  band. 
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Figure  10:  Semiconductor  band  structure  for  the  electron  energy  as  a  function  of  the  k  electron 
momentum  vector. 


Photon  absorption  processes  (i.e.  inter-band  or  intra-band  electron  excitation)  should 
conserve  the  momentum  of  the  electron.50  When  absorption  of  electromagnetic  radiation  occurs 
in  a  material,  the  absorption  can  be  expressed  as  a  function  of  the  radiation  energy: 

A  =  a(/iv)  (37) 


The  absorption  coefficient  a  is  the  rate  of  change  for  radiation  energy  as  a  function  of  path 
length,  L  along  the  axis  of  radiation  propagation  x 


a  = 


l(/zv)  dx 


(38) 


The  equation  of  the  absorption  coefficient  expression  defines  the  type  of  inter- band  electron 
transition  possible  for  a  semiconductor. 

Both  the  valence  and  conduction  bands  exist  as  a  continuum  of  energy  levels.  For  the 
valence  band  there  will  exist  a  highest  energy  occupied  electron  state  and  there  will  be  a 


26 


corresponding  lowest  energy  unoccupied  state  in  the  conduction  band.  The  transition  path 

between  this  highest  occupied  valence  state  and  lowest  unoccupied  conduction  band  state  can  be 

— ^ 

described  as  either  direct  or  indirect,  depending  upon  whether  the  electron  momentum  k 
changes  during  the  transition. 


a  b 


Figure  11:  Direct  (a)  and  indirect  (b)  semiconductor  excitations.  Adapted  from  Biakemore.43 


When  the  valence  band  maximum  and  conduction  band  minimum  energies  occur  at  the 
same  momentum  as  shown  in  Figure  11a,  the  band  gap  transition  is  direct.  For  direct  transition 
semiconductors  the  spectral  absorbance  is  given  by  the  expression: 

a(hv)  =  A*(hv-Egap)1/2  (39) 

where  A*  is  a  constant  and  Egap  is  the  semiconductor  band  gap  energy  difference.  Examples  of 
direct  excitation  semiconductors  are  cadmium  sulfide  (CdS)  and  gallium  arsenide  (GaAs).43 

W  hen  the  valence  and  conduction  band  extrema  occur  at  different  momenta,  as  shown  in 
Figure  lib,  the  transition  is  indirect.  Indirect  semiconductor  spectral  absorbance  is  given  by  the 
expression: 
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(40) 


a(/zv)  =  B"{hv-Ega^ 

where  B*  is  also  a  constant.  Electron  momentum  is  conserved  during  indirect  transitions  by 
coupling  of  the  electron  excitation  absorbance  to  either  the  emission  or  absorption  of  a  lattice 
phonon  vibration.  Examples  of  indirect  excitation  semiconductors  include  the  elements  silicon 
and  germanium.43 

1.3 Titanium  Dioxide 
1.3.1  General  Properties 

The  element  titanium  is  the  9th  most  abundant  in  the  Earth's  crust.51  The  Rev.  William 
Gregor  first  noted  titanium  compounds  in  1791  and  in  1795  the  element  itself  was  recognized 
and  formally  named  by  M  artin  Klaproth.52  Titanium  metal  and  its  alloys  possess  a  high  strength- 
to-weight  ratio  and  resistance  to  corrosion  which  have  led  to  their  widespread  use  for  the 
fabrication  of  aircraft,  missile,  and  submarine  parts.  These  mechanical  properties  and  titanium's 
biocompatibility  make  it  useful  for  fabrication  of  medical  implant  devices  such  as  those  used  for 
orthopedic  repair  and  dental  reconstruction.53,54 

Titanium  dioxide  (or  titania)  has  been  observed  extra  terrestrially  in  the  spectra  of  M  -type 
red  stars  and  in  the  composition  of  lunar  soil  samples  returned  during  the  Apollo  missions.51 
Titanium  dioxide  is  a  widely  used  industrial  commodity.  Until  2004  the  U.S.  Department  of 
Commerce  prepared  an  annual  report  of  domestic  titanium  dioxide  production  and  stockpiles  and 
in  2003  the  total  annual  U  .S.  production  was  reported  as  approximately  1.6  million  tons.55  While 
titanium  dioxide  is  found  in  a  wide  variety  of  consumer  products  ranging  from  cosmetics  to 
sunscreen  formulations56,  its  dominant  application  presently  is  in  the  production  of  paint.  While 
macroscopic  crystals  of  titanium  dioxide  are  transparent,  particles  sized  below  the  wavelength 
range  of  visible  light  (approximately  300  to  700  nanometers)  disperse  light  to  an  extent  that 
those  particles  appear  to  be  white.  This  optical  property  of  sub-micron  titanium  dioxide  particles 
and  their  perceived  low  toxicity  made  this  material  an  attractive  alternative  to  lead  based 
pigments  that  were  banned  in  the  United  States  in  1978  owing  to  the  neurotoxicity  of  lead.57, 58 
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While  titanium  dioxide  is  not  known  to  possess  the  neurological  risks  associated  with  lead 
compounds,  there  is  an  accumulating  body  of  research  that  indicates  that  nanoparticulate  titania 
may  hold  other  risks.  Chronic  inhalation  of  nanoparticulate  titanium  dioxide  has  been  linked 
with  pulmonary  disease.59"67 

The  shift  from  lead-based  white  paint  pigments  to  titanium  dioxide  inadvertently  lead  to 
the  discovery  of  the  photocatalytic  properties  of  titanium  dioxide.  The  organic  binder 
compounds  in  early  titanium  dioxide  pigmented  paints  were  observed  to  decompose  following 
exposure  to  sunlight,  resulting  in  'chalking'  of  the  paint  surface.68  Carl  Renz  published  several 
papers  in  the  early  1920' s  from  studies  where  he  observed  that  titanium  dioxide  particles 
suspended  in  glycerine  would  change  color  from  white  to  blue-gray  or  black  following  exposure 
to  sunlight.69  The  titania  color  change  noted  by  Renz  was  an  indication  that  the  titanium  cations 
had  undergone  a  photoinitiated  reduction  from  Ti4+  to  Ti3+  accompanied  by  extraction  of 
titanium  dioxide  lattice  oxygen  to  oxidize  the  glycerine  suspension  agent. 

Titanium  dioxide  has  been  a  popular  subject  of  photocatalysis  research  for  decades. 
Henrich  and  Cox  described  titania  as  "almost  the  prototypical  transition- metal  oxide  in  surface 
studies...  ”70  in  1994,  and  the  quantity  of  titanium  dioxide  research  since  that  publication  has  not 
diminished.  The  popularity  of  this  oxide  is  likely  caused  by  its  low  cost,  ease  of  preparation  and 
handling,  and  wide  range  of  intriguing  physical  properties. 

Titanium  dioxide  is  commonly  found  in  nature  as  one  of  three  crystalline  polymorphs: 
anatase,  rutile,  and  brookite.  Under  high  pressure  (4—10  GPa)  titanium  dioxide  has  been 
reported  to  undergo  transformation  to  columbite,  baddeleyite,  and  fluorite  forms,  but  these 
polymorphs  are  not  known  to  play  a  role  in  nanoparticle  photo  catalysis.71'75  The  polymorphs 
encountered  in  nanoparticle  photocatalyst  research  are  almost  exclusively  anatase  and  rutile. 
When  examined  in  bulk  form  at  ambient  temperatures  and  pressures,  rutile  is  the  most 
thermodynamically  stable  form  of  titanium  dioxide,  although  anatase  is  thought  to  be  the  most 
photo  catalytically  active  polymorph.71, 76  78  Both  anatase  and  brookite  will  convert  to  rutile 
upon  heating,  although  the  conditions  for  this  conversion  are  thought  to  be  dependent  upon 
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particle  size79,80  with  the  thermodynamically  favored  polymorph  shifting  from  rutile  to  anatase  at 
particle  sizes  below  14  nm  (1.4  x  10'8  m). 81-83  The  unit  cells  for  anatase  and  rutile  are  given  in 
Figure  12. 
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Figure  12:  Unit  cell  structures  for  the  rutile  and  anatase  polymorphs  of  titanium  dioxide.  Reprinted 
from  Surface  Science  Reports  87,  UlrikeDiebold,  The  surface  science  of  titanium  dioxide,  53-229,  2003, 
with  permission  from  E  Isevier.84 


Titanium  dioxide  is  an  n-type  semiconductor.  Both  the  rutile  and  anatase  polymorphic 
forms  of  titanium  dioxide  are  indirect  band  gap  semiconductors.85  The  band  gap  energies  for  the 
three  natural  polymorphs  have  been  determined  by  several  investigators  through  spectroscopic 
measurements  of  single  crystal  or  nanoparticle  samples.  Selections  of  these  observations  of  band 
gap  energy  are  presented  in  Table  1: 


Ti02  Band  Cap  Energy,  Eg/  dJ 

Investigatui-s) 

Anatase 

Rutile 

Bnookite 

Kim,  etal. 8b 

3.22 

Shibata,  et  al. 87 

3.20 

3.26 

Koelsch,  etal. 88 

3.34,3.39 

3.40 

Kavan,  et  al. 85 

3.2 

3.0 

Tang,  etal. 89 

3.2 

Table  1:  Titanium  dioxide  single  crystal  and  nanoparticle  band  gap  energies. 
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Although  it  is  possible  to  estimate  the  band  gap  energy  for  a  sample  of  nanoparticles 
from  diffuse  reflectance  spectroscopic  measurements,  particle  size  itself  can  influence  this 
parameter  through  quantum  confinement  effects. 

Upon  excitation,  an  electron  (e)  and  an  electron  hole  (h+)  can  become  bound  in  a 
hydrogen  like  entity  called  an  exciton.  Excitons  have  quantized  energy  levels  given  by  the 
equation:48 


E 

ex,n 


M 

2  2 


me  n 


(41) 


where  m  is  the  electron  rest  mass  (9.109  x  10'31  kg)90,  er  is  the  dielectric  constant  of  the 
semiconductor  material,  and  n  is  an  integer.  Mr  is  the  reduced  mass  of  the  exciton,  given  by  the 
equation: 
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M  m  m, 
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(42) 


where  me*  is  the  effective  electron  mass  and  rrih*  is  the  effective  hole  mass.  Effective  mass  is  the 
apparent  mass  of  an  entity  exhibited  during  transport  in  a  crystal  lattice.  Excitons  can  be  weakly 
bound  (large  distance  of  separation  between  electron  and  hole),  or  strongly  bound  (small  distance 
of  separation  between  electron  and  hole).  The  energy  necessary  to  create  an  exciton  need  not  be 
equivalent  to  the  band  gap  energy  for  the  semiconductor.48 

Roduner  described  the  considerations  that  determine  whether  strong  or  weak  quantum 
confinement  of  electrons  occur  within  a  nanoparticle.91  The  Bohr  radii  for  an  electron,  ae,  and  a 
hole,  ah,  are  given  by  the  equations: 


.  h2£Q£2 

4nm  e1 


(43) 
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(44) 


where  h  is  Planck's  constant  (6.626  x  10'34  J  •  s),  eo  is  the  vacuum  permittivity  constant  (8.854  x 
10 12  J  1  •  C2  •  m"1),  and  e  is  the  elementary  charge  (1.602  x  10  19  C).  When  the  radius  of  a 
semiconductor  particle  approaches  ae,  the  bandgap  energy  for  the  particle  increases.  The  lowest 
energy  transition  is  given  by: 


AE  ~ 


h2 

8  Mr2 


(45) 


where  r  is  the  particle  radius.  The  term  M  is  either  ju,  the  reduced  exciton  mass  for  cases  where 
r  <<  ae  and  r  <<  a/,  (strong  confinement),  or  M  is  the  effective  electron  mass  me  for  cases  where 
r  <<  de  and  r  >>  dh  (electron  confinement).  The  energy  change  for  the  lowest  non-degenerate 
transition  Z3  can  be  written  as: 


A E(Z3)  =  Eg-Eb  + 


0.67  h2 
8  Mr2 


(46) 


where  £ b  is  the  exciton  binding  energy  in  bulk  (r  — ►  qo)  and  £g  is  the  band  gap  energy.  For  either 
the  case  of  strong  confinement  or  electron  confinement  where  r  <<  ae,  the  non-degenerate 
energy  transition  can  be  approximated  by  Equation  (45)  and  Equation  (46)  can  be  rewritten  as: 


E  = 

8 


033 /r 
8  Mr2 


+  Eb 


(47) 


If  Eb  is  taken  to  be  constant,  clearly  as  the  particle  radius  r  decreases  the  band  gap  energy  Eg 
increases.  The  converse  of  this  statement,  band  gap  energy  decreases  (corresponding  to  a  UV- 
Vis  absorbance  maxima  shifted  to  longer  wavelength)  with  increasing  particle  size,  is  observed 
experimentally.  Kormann92  demonstrated  an  absorbance  red-shift  (i.e.  absorbance  shift  to  larger 
wavelength)  for  titanium  dioxide  during  preparation  of  colloidal  nanoparticles  over  thirty 
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minutes  of  particle  growth.  Other  quantum  confinement  induced  absorbance  maxima  shifts 
include  those  demonstrated  by  Wang  and  Herron  for  CdS  nanoparticles93,  Lockwood  for 
dissolving  Si  particles94,95,  and  King  etal.  for  titanium  dioxide  films  over  the  course  of  multiple 
annealing  cycles.96 

1.3.2  Semiconductor  Doping  and  T rap  States 

A  semiconductor  that  relies  solely  upon  the  magnitude  of  its  band  gap  energy  for 
excitation  of  electrons  to  the  conduction  band  is  an  intrinsic  semiconductor.  Often  the  band  gap 
energy  for  a  semiconductor  is  too  large  to  permit  adequate  population  of  the  conduction  band  to 
occur  from  room  temperature  thermal  excitation  alone.  One  approach  to  enhancement  of  valence 
band  to  conduction  band  electron  excitation  is  to  deliberately  introduce  a  trace  level  contaminant, 
known  as  a  dopant,  to  the  lattice  structure.  Dopant  materials  are  selected  for  their  abilities  to 
either  contribute  electrons  to  the  conduction  band  or  accept  electrons  from  the  valence  band  at 
energies  significantly  lower  than  the  intrinsic  band  gap  energy.  Semiconductors  that  rely  upon 
dopants  to  increase  their  electrical  conductivity  are  known  as  extrinsic  semiconductors. 

Dopant  atoms  replace  atoms  within  the  crystal  lattice  structure,  but  their  Bloch  equation 
solutions  are  different  from  those  of  the  semiconductor  lattice  in  which  they  reside.  Because  the 
semiconductor  dopants  have  different  Bloch  wavefunctions,  dopants  are  able  to  reside  within  the 
forbidden  band  gap  zone  of  the  semiconductor.  Dopant  electron  states  that  reside  near  the 
conduction  band  are  electron  donors,  and  the  doping  that  occurs  is  n-type,  where  ' n'  represents 
negative  charge  carriers  (i.e.  electrons). 
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Figure  13:  Phosphorus  n-doping  of  a  germanium  semiconductor.97 
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Figure  13  is  an  example  of  n-doping  of  a  germanium  semiconductor  crystal  by  the 
addition  of  phosphorus.  Like  the  carbon  in  a  diamond,  germanium  atoms  crystallize  with  a 
tetrahedral  bonding  geometry  (Figure  13a).  While  germanium  itself  is  an  intrinsic 
semiconductor,  the  addition  of  trace  amounts  (less  than  one  part  in  100, 000)97  of  phosphorus  to 
the  molten  germanium  during  crystallization  leads  to  the  incorporation  of  the  dopant  phosphorus 
into  the  germanium  lattice  (Figure  13b).  Asa  group  V  element,  phosphorus  is  pentavalent  and 
once  incorporated  into  the  germanium  lattice  it  has  an  unused  bonding  electron.  With  a  binding 
energy  of  0.0046  eV  to  the  P+  cation,  this  'extra'  electron  is  easily  excited  into  the  germanium 
conduction  band  at  room  temperature  (298K),  where  the  thermal  excitation  Ar7" is  only  0.025  eV 
(Figure  13c).97  The  phosphorus  dopant  electron  states  are  not  bound  by  the  Bloch  equations  for 
germanium,  and  these  phosphorus  donor  states  lie  just  below  the  conduction  band. 

Dopants  can  also  enhance  semiconductor  conductivity  through  acceptance  of  electrons 
from  the  valence  band.  Figure  14  describes  positive  charge  carrier,  p-doping  of  germanium.  In 
this  example  the  tetrahedral  germanium  bonding  geometry  (Figure  14a)  incorporates  a  boron 
atom  into  its  lattice  structure  during  crystallization.  As  a  group  III  element,  boron  is  trivalent 
and  there  is  a  disruption  in  the  germanium  tetrahedral  bonding  network,  indicated  by  a  broken 
line  in  Figure  14b.  The  boron  dopant  atoms  are  able  to  accept  valence  band  electrons  from 
adjacent  lattice  germanium  atoms.  Once  germanium  atoms  lose  a  valence  electron,  a  positive 
charge  carrier  entity,  or  a  'hole',  is  created  in  the  valence  band  and  electron  conduction  is 
facilitated. 


Figure  14:  Boron  p-doping  of  a  germanium  semiconductor. 
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Semiconductor  doping  is  relevant  to  the  electron  trapping  states  of  titanium  dioxide 
nanoparticles.  While  the  examples  of  semiconductor  doping  given  in  Figure  13  and  Figure  14 
were  accomplished  through  introduction  of  trace-level  foreign  materials  to  the  germanium 
crystal  lattice,  doping  effects  can  also  be  accomplished  through  chemical  reduction  of  a  small 
fraction  of  the  semiconductor.  Titanium  dioxide  can  be  made  non-stoichiometric  through 
heating  induced  loss  of  oxygen  from  its  crystal  lattice.  Figure  15  shows  a  schematic  for  the 
process  by  which  electrons  are  trapped  at  surface  states. 
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Figure  15:  Titanium  dioxide  surfacetrapping  of  electrons. 


FI  eating  of  the  titanium  dioxide  in  Figure  15a  leads  to  the  loss  of  oxygen  from  the  titania  crystal 
lattice.  The  oxygen  departing  the  lattice  leaves  behind  electrons  which  populate  the  Ti4+  3d 
orbitals  to  form  Ti3+  or  T i2+  as  shown  in  Figure  15b.  Because  the  T i3+  and  Ti2+  ions  are  not  part 
of  the  Bloch  equation  solution  that  defines  the  T iO 2  band  structure,  these  ions  can  energetically 
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exist  within  the  band  gap  as  sites  of  trapped  electrons.  Gopel  et  al.  have  shown  with  electron 
spectroscopy  techniques  that  the  titanium  3d  orbital  is  energetically  close  to  the  titanium  dioxide 
conduction  band.98  This  evidence  strongly  indicates  that  reduced  Ti3+  and  Ti2+  cations  are  the 
sites  of  shallow  trapped  electrons.  Ghosh  et  al.  identified  eight  electron  trapping  Ti3+  defect 
states  within  the  band  gap  of  single  crystal  titanium  dioxide  from  thermoluminescence  and 
thermally  stimulated  current  measurements.99  Because  of  the  large  surface-to-area  ratio  of 
nanoparticles,  the  physical  location  of  these  trap  sites  is  assumed  to  be  on  the  particle  surface.100 
Following  their  formation,  the  shallow  trapped  electrons  are  available  for  thermal  excitation  to 
the  conduction  band,  as  shown  in  Figure  15c.  The  formation  of  Ti3+  and  T i2+  defect  states  in  the 
titania  band  gap  is  itself  a  type  of  intrinsic  doping.  Within  semiconductor  nanoparticles  the 
number  of  trapped  charges  may  greatly  exceed  the  free  electron  population.101 

While  the  repository  of  trapped  electrons  within  the  band  gap  of  undoped  titanium 
dioxide  is  likely  the  titanium  cation  itself,  removal  of  oxygen  is  not  the  sole  means  of  trapping 
electrons.  The  surface  of  nanoparticles  feature  step  edges  and  corners,  shown  in  Figure  16  on 
page  38,  and  the  coordinate  unsaturated  titanium  cations  that  exist  at  these  surface  features  are 
susceptible  to  reductive  electron  trapping.  Other  sources  of  the  electrons  trapped  at  titanium  ion 
lattice  defects  include  electrochemical  reduction102'104  and  surface  adsorbates105.  Injection  of 
electrons  into  titanium  dioxide  nanoparticles  from  organic  dyes  is  the  mechanism  through  which 
titania  nanoparticles  have  been  incorporated  into  photovoltaic  cel  Is.18, 19'100, 106,107 

In  addition  to  study  of  intrinsic  Ti3+  trapped  electron  defects,  the  addition  of  dopants  such 
as  nitrogen108"114,  sulfur115"119,  fluorine120'124  to  titanium  dioxide  nanoparticles  is  an  active  area  of 
research.  Research  reported  in  the  literature  indicates  that  the  addition  of  these  dopants  either 
reduces  the  titanium  dioxide  band  gap  to  energies  that  can  be  spanned  by  visible  light,  or  these 
dopants  form  trap  states  in  the  band  gap  from  which  visible  light  or  thermal  excitation  can 
promote  electrons  to  the  conduction  band. 
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1.3.3  The  Titanium  Dioxide  Crystal  Surface 

The  unit  cell  of  a  perfect  crystal  lattice  extends  to  infinity  in  all  directions.  Within  a 
perfect  lattice,  an  electron  excited  into  a  conduction  band  state  will  propagate  infinitely  with  its 
associated  Bloch  wave  function.  This  perfect  crystal  lattice  is  a  construct  useful  for  the 
derivation  of  theory  but  it  does  not  describe  a  real  crystal  or  a  crystalline  nanoparticle. 

In  contrast  to  the  ideal  cases,  real  crystals  are  bounded  by  surfaces.  The  surfaces  of  a 
crystal  constitute  a  potential  energy  barrier  that  confines  mobile  charge  carriers  within  the  lattice. 
Electrons  that  have  acquired  sufficient  energy  within  the  lattice  are  able  to  overcome  the 
potential  energy  barrier  at  the  surfaces  and  escape  the  lattice.  Thermionic  electron  emission 
from  metals  is  a  manifestation  of  electrons  overcoming  the  lattice  potential  energy  barrier. 

Another  aspect  of  crystalline  surface  asymmetry  is  reconstruction.  The  atoms  or  ions 
present  on  the  surface  of  a  crystal  will  rearrange  so  as  to  minimize  the  energy  of  the  surface. 
This  shifting  of  bond  lengths  and  geometric  positions  for  the  surface  layer  will,  in  turn,  alter  the 
geometry  of  the  lattice  plane  immediately  adjacent  to  the  surface  layer.  The  truncation  of  a 
crystal  lattice  at  the  surface  can  leave  atoms  or  ions  incompletely  coordinated,  effectively  able  to 
bond  to  other  species.  Insufficiently  coordinated  surface  ions  can  result  in  a  reactive  crystal 
surface.  Quite  often  these  terminal  surface  ions  will  minimize  energy  by  reacting  with 
atmospheric  water  to  form  hydroxylated  surface  species.  A  surface  covered  with  highly  polar 
hydroxy  groups  then  becomes  hydrophilic,  i.e.  wettable.  Prevention  or  removal  of  surface 
coverage  with  atmospheric  components  represents  perhaps  the  greatest  challenge  facing  surface 
science  researchers,  necessitating  the  use  of  ultra-high  vacuum  approaches. 

A  nother  feature  of  the  crystal  surface  is  lack  of  uniformity.  If  the  surface  of  a  crystal  was 
formed  by  removal  of  the  atom  or  ion  flux  that  lead  to  the  crystal's  growth,  any  irregularities  that 
were  present  in  the  interior  of  the  lattice  structure  will  be  represented  on  the  surface,  much  like 
imperfections  on  a  painted  surface  that  result  from  underlying  irregularities.  If  the  crystal  lattice 
was  formed  by  cleavage  of  a  larger  crystal  on  a  lattice  plane,  it  is  unlikely  that  the  newly  formed 
surfaces  will  be  equivalent.  M  echanical  or  thermal  events  may  partially  remove  lattice  ions  from 
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the  surface.  Partial  loss  of  lattice  atoms  en  masse  can  either  lead  to  the  formation  of  terraces  or 
islands  on  the  surface.  Discrete  lattice  component  losses  from  the  surface  of  a  lattice  can  lead  to 
isolated  vacancies  that  interrupt  the  regularity  of  the  surface  pattern.  Bombardment  of  the  crystal 
surface  by  gas  molecules  can  lead  to 'foreign'  adatom  features  on  the  surface.  The  final  result  of 
material  transport  to  and  from  a  crystal  surface  is  an  irregular  terrain  possessing  many  potential 
sites  for  surface  reactivity  and/or  electron  trapping.  A  very  generic  example  of  some  of  these 
surface  irregularities  is  given  in  Figure  16. 


adatom  terrace 


Figure  1&  Surface  irregularities  on  a  crystal  lattice125 


The  surface  of  titanium  dioxide  possesses  several  possibilities  for  charge  trapping  that  are 
apparent  on  one  of  the  surface  reconstructions  of  the  rutile  polymorph.  Electron  trapping  that  is 
less  than  1  eV  below  the  conduction  band  edge  is  defined  as  shallow,  while  trap  states  that 
exceed  1.5  eV  beneath  the  conduction  band  edge  are  defined  as  deep.126 

Figure  17  illustrates  the  (110)  surface  of  the  rutile  titanium  dioxide  polymorph.  Other 
possible  surface  configurations  for  rutile  exist,  but  as  a  low  energy  configuration  the  (110) 
surface  represents  the  most  likely  surface  orientation  to  be  encountered.  This  surface  geometry 
possesses  two  structural  features  that  play  a  role  in  electron  trapping. 
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6f-Ti  H  adatom 


[110] 

[001]  ^  O  O  atoms  O  H  adatoms  (OH**) 

Figure  17:  Aspects  of  the  rutile  Ti02  surface.  Reprinted  with  permission  from  Matthiesen,  et  al., 

Observation  of  All  the  Intermediate  Steps  of  a  Chemical  Reaction  on  an  Oxide  Surface  by  Scanning 

T  unneiing  M  icroscopy.  ACS  Nano  2009,  A-j  .Copyright  2009  American  Chemical  Society.127 

Titanium  cations  within  a  titanium  dioxide  lattice  are  coordinated  in  a  distorted 
octahedral  geometry  with  oxygen  anions.  The  rutile  lattice  is  formed  by  a  repeating  orthogonal 
arrangement  of  these  octahedral  complexes.  At  the  rutile  surface,  some  of  the  complexes  are 
truncated  leading  to  the  formation  of  coordinate  unsaturated  Ti4+  cations,  indicated  by  ' 5f-T i '  in 
Figure  17.  These  exposed,  unsaturated  Ti4+  sites  are  Lewis  acid  sites,  capable  of  receiving 
electron  density  at  their  exposed  face.  In  addition  to  forming  bonds  with  reacting  species  on  the 
surface  of  the  nanoparticle,  these  sites  are  also  susceptible  to  reduction,  i.e.  they  can  function  as 
electron  traps.  Thompson  has  referred  to  the  Ti3+  species  within  titanium  dioxide  nanoparticles 
as  an  electron  trapped  in  the  band  gap  region.126,127 

Other  octahedral  complexes  are  oriented  with  their  central  axes  normal  to  the  lattice 
surface,  leading  to  the  formation  of  parallel  rows  of  oxygen  ions  that  bridge  the  Ti4+  cations 
immediately  below.  These  rows  of  bridge  bonded  oxygen  cations  on  the  face  of  the  lattice  can 
be  thermally  excited  to  the  point  of  leaving  the  lattice.  Titanium  cations  left  behind  by  thermal 
oxygen  loss  from  the  crystal  surface  acquire  an  electron  and  are  reduced  from  the  Ti4+  to  Ti3+ 
state.  These  reduced  titanium  cations  constitute  surface  sites  of  trapped  electrons.  Examples  of 
bridge  oxygen  vacancies  are  indicated  in  Figure  17  by  broken  line  sites  where  oxygen  ions 
formerly  resided.  One  such  site  on  the  far  left  side  of  the  diagram  is  labeled  'Obr  vac.'.  In 
addition  to  the  possibility  of  being  thermally  driven  off  of  the  lattice  entirely,  the  bridge  oxygen 
ions  in  Figure  17  can  function  as  Lewis  bases  and  accept  protonation.  A  protonated  bridge 
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oxygen  is  indicated  on  the  far  right  side  of  the  Figure  17  by  the  'H  adatom'  label.  Hydroxyl 
groups  located  on  the  surfaces  of  titanium  dioxide  nanoparticles  can  function  as  hydrogen 
bonding  sites,  capable  of  promoting  surface  adsorption  of  complementary  species. 
Szczepankiewicz  has  also  identified  surface  hydroxyls  as  a  site  for  delocalized  shallow  electron 
trapping.128 

Diagrams  such  as  Figure  17  can  lead  to  a  mistaken  conclusion  that  a  crystal  lattice,  even 
a  flawed  one,  is  a  static  and  passive  terrain.  The  reality  is  that  the  surface  of  crystalline  titanium 
dioxide  is  a  highly  dynamic  and  reactive  locale.  It  is  even  possible  to  follow  the  movements  of 
oxygen  atoms  and  adsorbates  from  external  sources  in  real  time  by  electron  tunneling 
microscopy,  as  has  been  recently  demonstrated  by  M  atthiesen.127 
1.3.4  Titanium  Dioxide  Nanoparticles  and  Shallow  Electron  Trap  States 

Nanoparticles  have  attracted  intense  research  interest  in  part  due  to  their  phenomenally 
high  surface  areas  which  can  reach  hundreds  of  square  meters  per  gram  of  material.  This  high 
ratio  of  area  to  mass  produces  an  unusually  large  stage  upon  which  chemical  reactions  can  be 
orchestrated.  The  increased  surface  area  of  nanoparticles  magnifies  the  contribution  of  crystal 
lattice  defects  to  the  particles'  chemical  reactivity.  A  large  surface  area  for  nanoparticles 
translates  into  a  large  number  of  surface  electron  trapping  states.  Some  researchers  have  noted 
that,  on  the  surface  of  semiconductor  nanoparticles,  the  density  of  trap  states  may  greatly 
outnumber  the  population  of  mobile  electrons.101 

Ghosh  identified  six  definite  and  two  likely  shallow  trap  states  ranging  from  0.27  to  0.87 
eV  below  the  conduction  band  edge  in  single  crystal  titanium  dioxide  from  thermoluminescence 
and  thermally  stimulated  current  observations,  and  speculated  that  these  trap  states  were  likely 
Ti3+  reduced  from  Ti4+  in  the  crystal  lattice.99, 129  The  low  energy  gap  between  these  states  and 
the  conduction  band  makes  electrons  residing  in  these  traps  thermally  available  for  transport  to 
the  conduction  band  and  participation  in  chemical  reaction.  Because  shallow  trap  states  are  the 
result  of  either  unsaturated  Ti4+  or  bound  hydroxyl  species  attached  to  the  titanium  dioxide 
nanoparticle,  these  electron  traps  are  a  surface  phenomenon.  The  trap  states  can  function  as 
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intermediate  repositories  for  electron  storage  that  permit  surface  catalysis  to  proceed  at  energies 
far  below  those  required  for  direct  valence  to  conduction  band  excitation. 

1.4 Titanium  Dioxide  Nanoparticle  Preparation 
1.4.1  Liquid  Preparation  Techniques 

Preparation  of  metal  oxides  from  liquid  reagent  precursors  is  one  of  the  simplest 
approaches  to  forming  nanoparticles.  Synthesis  of  particles  in  a  liquid  environment  affords 
possibly  the  greatest  level  of  control  over  the  chemical  environment  under  which  the  particles  are 
formed.  These  procedures  are  performed  with  common  laboratory  equipment,  simple  techniques, 
and  relatively  inexpensive  reagents.  Because  of  this  simplicity  the  majority  of  titanium  dioxide 
nanoparticle  research  is  conducted  with  materials  prepared  from  liquid  synthesis  techniques. 

To  produce  a  nanoparticle  is  to  grow  a  crystal.  The  challenge  and  principle  difficulty  is 
to  control  the  size,  size  distribution,  and  geometry  of  the  nanoparticles  produced.  M  onodisperse 
particle  growth,  or  the  formation  of  discrete  particles  of  identical  size  and  shape,  is  the  objective 
of  most  nanoparticle  syntheses.  Two  approaches  to  nanoparticle  synthesis  have  evolved:  forced 
hydrolysis  and  colloidal  particle,  i.e.  sol-gel  techniques.  Matijevic  and  Sapieszko  described  the 
general  conditions  necessary  to  monodisperse  nanoparticle  growth  from  solutions  by  forced 
hydrolysis.130 

The  first  stage  of  nanoparticle  growth  within  a  solution  is  the  formation  of  clusters  from 
solvated  ions.  Following  cluster  formation,  the  formation  of  nanoparticle  nuclei  from  the 
clusters  must  be  accomplished.  If  the  conditions  under  which  these  initial  stages  occur  are  not 
optimal,  uncontrolled  coagulation  occurs. 

Nanoparticle  growth  relies  upon  minimization  of  the  electrostatic  repulsion  between  the 
precursors  in  solution.  This  requirement  is  generally  accomplished  by  control  of  the  reaction 
solution  pH  to  adjust  the  solvent  environment  to  a  point  isoelectric  with  the  surface  precursor 
charge.  Another  approach  is  to  increase  the  solution  ion  strength  to  a  level  where  the  precursor 
surface  charges  are  screened  from  one  another.  Typical  approaches  to  ion  strength  adjustment 
include  addition  of  chlorides  and  sulfates. 
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Reactant  concentration  and  time  for  particle  growth  also  play  important  roles  in  the 
quality  of  monodisperse  particle  growth..  If  precursor  concentration  is  too  high,  uncontrolled 
agglomeration  of  crystallites  will  occur  with  rapid  precipitation.  Once  optimal  reactant 
concentrations  have  been  determined,  excessive  prolongation  of  crystal  growth  will  lead  to  an 
increase  of  variation  in  crystal  size. 

M  atijevic  et  al.  are  believed  to  have  reported  the  first  instance  of  monodisperse  titanium 
dioxide  particle  synthesis  from  a  titanium  tetrachloride  precursor  with  sulfate  ion  strength 
adjustment.131  This  work  was  a  continuation  of  earlier  studies  on  metal  oxide  particle  formation 
from  hydrous  sols.132  The  M  atijevic  titania  particles  ranged  from  1  to  4  pm  (1  to  4  x  10'6  m) 
which  exceeds  most  definitions  of  a  nanoparticle,  but  his  techniques  for  particle  growth  were 
used  in  more  recent  work. 

Kim  et  al.  reported  upon  the  preparation  of  a  stable TiOCh  aqueous  solution  from  which 
nanoparticles  were  latter  prepared.133  Their  procedure  began  with  slow  addition  of  water  to 
titanium  tetrachloride  cooled  to  0  QC  followed  by  dilution  in  distilled  water.  The  resulting 
precursor  reagent  was  reported  stable  for  up  to  one  year  when  stored  at  room  temperature. 
Nanoparticle  formation  was  accomplished  simply  by  heating  theTiOCh  solution  to  100  QC.  The 
particles  produced  were  the  anatase  titania  polymorph,  exhibiting  moderate  crystallinity  without 
further  treatment.  By  altering  the  rate  of  solution  heating,  rutile  nanoparticles  were  also  prepared 
from  the  same  TiOC l2  starting  material.  The  resulting  nanoparticles  were  reported  to  be 
spherical  with  a  mean  diameter  of  200  nm  (2.00  x  10'7  m). 

Other  researchers134"142  have  continued  forced  hydrolysis  of  titanium  tetrachloride  with 
variations  on  the  M  atijevic  techniques.  Addition  of  acid  catalysts,  use  of  nonaqueous  solvents, 
changes  in  precursor  concentrations,  different  crystal  growth  temperature  regimes,  and  variation 
in  nanoparticle  growth  time  have  all  been  applied  to  custom  nanoparticle  synthesis.  The  results 
of  these  studies  have  been  nanoparticles  with  sizes  from  10  to  50  nm  that  exhibit  high 
photoactivity.  In  all  of  the  cited  works,  the  nanoparticles  that  were  produced  were  agglomerated 
into  fused  growths  of  nanoparticles.  The  causes  were  likely  particle  collisions  while  growing  in 
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solution  and/or  post  synthesis  heat  treatments  intended  to  increase  crystallinity  and  alter  the 
polymorph  composition  of  the  nanoparticles. 

Sol  gelation  techniques  for  nanoparticle  synthesis  have  become  common  in  the  literature. 
A  sol  is  a  suspension  of  solvated  colloidal  particles  which  do  not  precipitate.  If  these  particles 
are  allowed  to  grow  until  their  structures  merge,  a  gel  is  formed  where  solvent  is  captured  within 
void  spaces  of  the  gel  network.  The  structure  remaining  after  removal  of  the  solvent  is  a  high 
surface  area  solid. 

The  sol-gel  approach  to  nanoparticle  synthesis  is  based  upon  synthesis  of  the  colloidal 
particles  followed  by  cessation  of  particle  growth  before  a  macromolecular  gel  network  is 
formed.  The  precursor  reagent  is  an  organometallic  compound,  usually  an  alkoxide,  that  is 
easily  hydrolyzed  to  form  a  hydroxide  and  an  alcohol143: 

M(OR)n  +xH20  -»  M(OR)nx(OH)x  +  xROH  (48) 

Formation  of  the  particles  proceeds  by  a  condensation  reaction  where  either  water  or  alcohol  is 
eliminated  with  the  addition  of  each  metal  hydroxide  unit: 

M(OR),x(OH)x  +  M(OR)nx (OH)x  ->  M 2 (OR) 2(n _x) (OH) 2x  jO  +  H20  (49) 

M(OR)n.x(OH)x  +  M(OR)nx(OH)x  -»  M2(OH)2(n_x)(OR)2x  ,0  +  ROH  (50) 

Typically  these  reactions  are  conducted  with  the  metal  alkoxide  and  the  water  reactant 
dispersed  into  anhydrous  organic  solvents,  with  alcohols  being  the  most  frequently  encountered 
solvent  system  in  use.  The  reagent  concentrations  and  molar  ratio  must  be  optimized  to  promote 
homogeneous  sol  growth  rather  than  uncontrolled  precipitation.  Ogihara143  identified  the 
optimal  range  of  alkoxide  precursor  concentration  as  0.01  ~  0.2  mol/dm3  with  a  corresponding 
water  reactant  concentration  range  of  0.05  ~  0.5  mol/dm3  for  monodisperse  titanium  dioxide 
nanoparticle  preparation.  The  reactants  must  be  completely  mixed  before  the  onset  of  particle 
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nucleation.  Insoluble  impurities  must  be  removed  from  the  reactant  solutions  prior  to  use,  as 
extraneous  particulates  serve  as  unintentional  nucleation  sites.  Finally,  the  stability  of  the  colloid 
suspension  must  be  maintained  either  through  the  use  of  steric  or  electrostatic  repulsion. 

Like  forced  hydrolysis,  much  work  has  been  done  to  find  optimal  sol-gel  reaction 
conditions  for  control  of  particle  size,  geometry,  and  composition.  The  alkoxide  precursors  used 
have  generally  been  titanium  tetraethoxide  (TTE)144,  titanium  tetra-n-butoxide  (TTB)145"150,  or 
titanium  tetraisopropoxide  (TTIP)148, 151'157.  The  alkoxide  is  usually  prepared  in  alcohol  and 
added  slowly  to  either  concentrated  acid  catalyst  (hydrochloric,  acetic,  or  nitric  acids)  or 
combined  with  water  dissolved  in  alcohol  solution.  Salt  additives  are  sometimes  used  to  direct 
the  phase  composition  of  the  nanoparticles.  The  sol  formed  following  the  initial  mixing  of 
reactants  is  allowed  to  'age'  to  promote  nanoparticle  growth  for  intervals  ranging  from  hours  to 
weeks.  Deshpande  et  al.149  even  examined  the  effects  of  particle  growth  over  one  year. 

The  product  of  many  sol-gel  preparations  are  amorphous  anatase  nanoparticles. 
Crystallization  of  the  nanoparticles  and  partial  conversion  to  the  more  thermodynamically  stable 
rutile  polymorph  is  usually  accomplished  through  heat  treatment,  referred  to  as  calcining.  The 
dried  amorphous  particles  are  heated  to  temperatures  ranging  from  100  to  900  QC  for  several 
hours,  with  samples  of  the  particles  taken  at  regular  intervals  for  x-ray  diffraction  analysis  of 
crystallinity  and  composition. 

One  negative  outcome  of  calcination  is  that  at  higher  temperatures,  smaller  nanoparticle 
coalesce  into  larger  particles  with  a  loss  of  surface  area.  M  any  researchers  have  attempted  to 
avoid  this  alteration  of  nanoparticle  size  by  inducing  crystallinity  with  milder  thermal  treatment 
of  the  nanoparticles  while  they  are  still  in  solution.150, 153,156 

Solvent  based  nanoparticle  synthesis  presents  a  relatively  simple  approach  to  production 
of  particles  of  high  purity  and  controlled  composition.  Unfortunately,  preparation  of 
monodisperse  (i.e.  uniform  sized)  non-aggregated  nanoparticles  through  this  technique  remains 
difficult. 
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142  Vapor  Pyrolysis 

The  commercially  prepared  titanium  dioxide  nanoparticle  catalyst  frequently  used  as  a 
reference  for  comparison  is  Degussa  (now  EvonikA.G.)  P25.  This  nanoparticle  preparation  has 
a  mean  particle  size  of  25  nm  (2.5  x  10"8  m)  and  is  estimated  to  be  80%  anatase  and  20%  rutile. 
P25  titanium  dioxide  is  prepared  by  flame  hydrolysis  of  titanium  tetrachloride.  Flame  hydrolysis 
operates  by  injection  of  a  precursor  reagent  or  reagents  into  a  flame.  Oxidation  of  the  precursors 
takes  place  within  the  flame  and  particles  of  oxide  accrete,  sinter,  and  are  ejected.  Similar 
processes  are  used  industrially  for  titanium  dioxide  paint  pigment  manufacture158  and  fabrication 
of  telecommunications  optical  fiber.159, 160  The  two  methods  frequently  utilized  by  researchers 
for  pyrolytic  synthesis  of  titanium  dioxide  nanoparticles  are  the  tube  furnace  reactor  method  and 
the  laminar  diffusion  flame  torch. 

The  tube  furnace  reactor158  operates  by  accepting  vaporized  titanium  precursor  (typically 
titanium  tetraisopropoxide)  from  a  bubbler-nebulizer  apparatus  into  a  quartz  tube  that  passes 
through  an  electric  furnace.  W  ithin  the  furnace  heated  zone,  the  vaporized  precursor  reacts  with 
oxygen  and  titanium  dioxide  particles  are  formed.  During  their  passage  through  the  furnace,  the 
oxide  particles  grow  and  sinter.  The  nanoparticles  exiting  the  tube  furnace  are  captured  in  an 
electrostatic  dust  precipitator  that  uses  a  high  voltage  electric  field  to  force  the  electrically 
charged  nanoparticles  to  stick  to  an  electrode.  A  schematic  of  the  tube  furnace  reactor  is  shown 
in  Figure  18.  A  similar  process  for  continuous  generation  of  titanium  dioxide  nanoparticles  from 
a  plasma  induced  reaction  of  titanium  tetraisopropoxide  and  water  within  a  tube  reactor  was 
recently  described  by  Chen  etal.161 
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Figure  18:  Titanium  dioxide  tube  furnace  pyrolysis  reactor.158 


Pyrolysis  of  titania  precursor  within  a  torch  flame,  followed  by  nanoparticle  collection  in 
an  exhaust  hood  or  upon  a  water-cooled  metal  plate,  is  a  more  commonly  found  approach. 
Pratim  Biswas162'176  has  extensively  researched  the  operational  parameters  of  laminar  flow 
torches  used  to  hydrolyze  titanium  dioxide  alkoxide  precursors  (usually  titanium 
tetraisopropoxide).  An  example  schematic  of  Biswas’  titanium  dioxide  pyrolysis  apparatus  is 
given  in  Figure  19: 


Figure  19:  Flame  pyrolysis  nanoparticle  generation.168 
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The  laminar  flow  torch  used  for  vapor  pyrolysis  consists  of  concentric  quartz  tubes,  into  which 
individual  gaseous  components  of  the  hydrolysis  reaction  are  introduced.  Mixing  of  fuel, 
oxidizer,  and  precursor  vapor  occur  immediately  above  the  torch  head.  One  innovation  of  the 
Biswas  apparatus  is  the  air-cooled  quench  ring.  Thisfeature  is  used  to  direct  jets  of  air  onto  the 
hydrolysis  flame  for  the  purpose  of  reducing  nanoparticle  aggregation  through  sintering.  Particle 
collection  from  the  flame  is  accomplished  by  drawing  air  into  a  conical  metal  hood  affixed  to  the 
top  of  the  quench  ring  and  passing  the  air/nanoparticle  stream  through  an  electrostatic  particle 
precipitator.  In  addition  to  the  precipitator,  instrumentation  for  characterizing  nanoparticle 
morphology  through  light  scattering  behavior  can  also  be  attached  to  the  downstream  particle 
collection  train.  An  alternative  to  the  particle  quench  ring  /  collection  hood  assembly  is  a  water 
cooled  metal  plate.  The  metal  collection  plate  is  attached  to  glass  or  silicon  plates  and  this 
apparatus  is  used  to  study  nanoparticle  films,  such  as  would  be  used  for  electrodes  in 
photovoltaic  cell  fabrication. 

The  high  temperatures  associated  with  flame  pyrolysis  have  several  consequences.  The 
high  temperature  environment  in  which  nanoparticles  are  formed  within  the  flame  leads  to 
greater  crystallinity  with  less  need  for  post  synthesis  calcination.  Particle  geometry  from  flame 
synthesis  is  almost  exclusively  spherical.  The  distribution  of  temperatures  that  exist  within  the 
torch  flame  can  cause  a  polydisperse  nanoparticle  size  distribution  and  cause  extensive  particle 
sintering  into  agglomerates.  Sintering  and  agglomeration  can  be  somewhat  mitigated  through 
alteration  of  operating  parameters  for  the  pyrolysis  apparatus.  While  removal  of  particle  size 
dispersion  cannot  be  eliminated  completely,  the  continuous  operation  mode  of  flame  pyrolysis 
nanoparticle  generation  lends  itself  to  application  of  particle  size  selection  apparatus  such  as  the 
multiple-stage  differential  mobility  analyzer  described  by  Chen  et  al.177  and  applied  to  flame 
pyrolysis  by  Jiang,  Chen,  and  Biswas.168 
L4.3  Laser  Techniques 

Synthesis  of  nanoparticles  with  laser  induced  vaporization  carries  several  advantages 
over  vapor  pyrolysis  conducted  at  or  near  atmospheric  pressure.178  Unlike  a  torch  flame,  a  laser 
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beam  does  not  require  or  generate  near-atmospheric  gas  pressures.  A  laser  can  be  set  up  outside 
of  a  closed  chamber  and  its  energy  output  can  be  directed  by  optics  into  that  chamber  to  heat  a 
sample  contained  under  vacuum  conditions  or  some  non-atmospheric  gas  regime.  This  ability  to 
remove  the  source  of  material  heating  from  the  experiment  locale  reduces  the  likelihood  of 
contamination  of  both  the  precursor  materials  and  the  final  nanoparticle  product  by  making 
greater  control  of  the  experimental  conditions  possible. 

Lasers  are  capable  of  generating  greater  temperatures  than  conventional  flames.  Direct 
vaporization  of  refractory  metal  oxides  and  high  melting  point  metals  becomes  possible  with 
laser  heating.  Laser  heating  of  nanoparticle  precursor  materials  in  an  enclosed  chamber  opens 
the  possibility  of  reactive  synthesis  by  simply  introducing  small  quantities  of  gas.  A  pure  metal 
sample  can  be  vaporized  to  form  nanoparticles  and  either  during  the  vaporization  or  immediately 
afterwards,  a  small  quantity  of  oxygen  can  be  introduced  to  form  metal  oxide  nanoparticles. 
Changing  the  gas  to  nitrogen  can  produce  the  corresponding  metal  nitride.  Altering  the  pressure 
of  a  gas  or  gases  admitted  to  the  experiment  chamber  for  reactive  synthesis  can  influence  the 
stoichiometry  of  the  resulting  compounds  and  profoundly  change  their  characteristics. 

The  configuration  of  a  chamber  for  synthesis  of  nanoparticles  is  not  complex.  A  generic 
system  is  diagrammed  in  Figure  20.  A  pellet  or  rod  of  the  precursor  material  is  mounted  in  the 
synthesis  chamber  on  an  apparatus  that  allows  for  rotation  positioning  during  operation  of  the 
chamber.  Unless  the  chamber  is  designed  in  a  manner  to  place  the  enclosed  precursor  target  on 
the  laser  beam  axis,  front-surface  mirrors  are  used  to  direct  the  beam  into  the  chamber  through  a 
window  selected  for  its  transparency  to  the  laser  operating  wavelength.  The  nanoparticles  that 
condense  from  the  laser  heated  vapor  can  be  collected  on  a  stationary  surface  inside  the  chamber 
(Figure  20a)  or  can  be  withdrawn  outside  the  chamber  through  a  sampling  port  (Figure  20b) 
operated  at  a  lower  pressure  than  the  interior  of  the  chamber.  Chamber  pressure  is  controlled  by 
an  attached  vacuum  system  and  the  controlled  introduction  of  reagent  gases. 
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Figure  20:  A  typical  setup  for  synthesis  of  nanoparticles  with  laser  heating.  Nanoparticles  are 
captured  on  a  stationary  collector  surface  (a)  or  are  removed  from  the  chamber  during  synthesis  (b). 


Nd  -  YAG  (neodymium  -  yttrium-  aluminum  -  garnet)  solid-state  lasers  are  a  common 
choice  for  nanoparticle  synthesis.178"182  Although  these  lasers  operate  on  a  fundamental 
wavelength  of  1.06  qm,  generating  higher  order  harmonics  can  be  used  to  create  wavelengths  of 
0.213,  0.265,  0.355,  and  0.530  pm.183  Nd-YAG  lasers  can  be  operated  in  continuous  wave 
mode,  although  greater  power  output  is  achieved  through  pulsed  operation. 

Another  laser  used  for  nanoparticle  synthesis  is  the  carbon  dioxide  gas  laser.  These 
lasers  operate  on  a  fixed  wavelength  10.6  pm  in  the  infrared  region  of  the  electromagnetic 
spectrum.  Carbon  dioxide  lasers  are  capable  of  higher  power  output  than  their  Nd-YAG  solid- 
state  counterparts,  with  power  levels  of  1000  watts  possible  with  continuous  mode  operation.184 
This  high  power  capability  makes  C02  lasers  a  frequent  choice  for  industrial  applications  where 
cutting  or  melting  of  material  is  required. 

There  are  two  approaches  to  laser  operation  for  nanoparticle  synthesis.  Ablation  is  the 
result  of  pulsed  laser  operation.  In  this  mode,  a  portion  of  the  target  material  is  rapidly  heated, 
liquefied,  and  vaporized  during  the  'on'  portion  of  the  laser  pulse  cycle.  With  the  short  duration 
of  a  laser  pulse,  usually  between  10"6  ~  10"9  seconds,  the  removal  of  material  from  the  target 
surface  occurs  almost  explosively.  With  each  pulse,  the  vaporized  material  that  departs  the 
target  imparts  momentum  onto  the  bulk  material  it  leaves  behind.  During  the  'off  portion  of  the 
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laser  pulse  cycle,  the  target  material  has  some  time  to  dissipate  energy  from  the  immediately 
preceding  pulse  burst. 

Continuous  operation  of  a  laser  directed  onto  a  target  leads  to  vaporization  operation.  In 
this  mode,  the  precursor  target  material  remains  molten  in  the  region  where  the  laser  is  focused, 
and  the  molten  material  vaporizes  continuously.  With  no  interruption  in  material  vaporization, 
the  environment  of  the  synthesis  chamber  achieves  saturation  more  rapidly  and  maintains 
conditions  for  continuous  condensation  of  nanoparticles. 

The  dynamics  of  particle  condensation  for  laser  synthesis  of  nanoparticles  can  be 
complex.  It  is  not  unusual  for  solid  fragments  of  the  target  material  to  be  ejected  during  heating 
and  become  intermixed  with  the  nanoparticles  that  condense  from  the  gas  phase.  The 
wavelength  of  the  laser  beam  will  determine  the  nature  of  the  gas  environment.  For  titanium 
dioxide,  M  atsubara185  and  Kawasaki186  note  that  if  the  wavelength  of  the  laser  is  too  great  to 
impart  energy  greater  than  the  precursor  target  band  gap  (i.e.  the  titanium  dioxide  3.0  eV  band 
gap),  then  the  vapor  environment  that  forms  the  nanoparticles  will  be  molecular  species  such  as 
Ti02  and  Ti203.  A  Nd-YAG  laser  operated  at  266  nm  wavelength  will  have  a  photon  energy 
4.66  eV,  which  exceeds  the  titania  precursor  target's  band  gap  and  makes  possible  complete 
dissociation  of  the Ti02  solid  into  gaseous  atoms  and  ions. 

Laser  nanoparticle  synthesis  does  have  some  disadvantages.  The  cost  of  the  equipment  is 
high  compared  to  liquid  reagent  nanoparticle  synthesis  or  vapor  pyrolysis.  The  nanoparticle 
yields  are  generally  smaller,  and  particle  aggregation  still  occurs. 

L5 1  nfrared  Spectroscopy  of  Electrons 

L  5.1  Evidence  for  the  Infrared  Detection  of  Conduction  Band  Electrons 

The  elements  silicon  and  germanium  are  semiconductors,  with  band  gap  energies  of  1.12 
and  0.67  eV,  respectively.187  The  transparency  of  silicon  and  opacity  of  germanium  to  infrared 
radiation  were  utilized  by  Becker  and  Fan188  to  demonstrate  the  absorption  of  infrared  radiation 
by  mobile  charge  carriers  (i.e.  delocalized  electrons)  in  solid  state  matter.  Their  technique  used 
photographic  plates  that  were  narrowly  sensitized  to  radiation  in  the  infrared  region  at  1.3  pm 


50 


(1.3  x  10"6  m)  wavelength.  The  infrared  sensitized  plates  were  either  covered  with  thin  sheets 
(0.3  mm  thickness)  of  elemental  silicon  or  germanium  and  exposed.  The  silicon  covered  plate 
showed  an  image  from  its  exposure  whereas  the  germanium  covered  plate  did  not.  The  infrared 
transmittance  of  bulk  germanium  begins  at  1.6  pm  wavelength  and  above,  so  the  cutoff  at  1.3  pm 
is  expected.  The  smaller  bandgap  of  elemental  germanium  makes  it  more  likely  to  absorb 
energy,  either  electromagnetic  or  thermal,  with  the  result  being  valence  to  conduction  band 
electron  transitions.  The  larger  band  gap  energy  of  silicon  requires  more  energetic  (shorter 
wavelength)  radiation  to  initiate  electron  transitions  and,  as  a  result,  the  absorbance  cutoff  for 
silicon  is  lower  at  wavelengths  of  1.0  pm  and  below. 

Becker  and  Fan  recorded  infrared  spectra  for  the  thin  sheet  of  silicon  over  the 
wavelengths  of  0.8  to  2.3  pm  (12500  to  4350  cm'1)  at  temperatures  from  77  K  to  663  K .  As  the 
temperature  of  the  silicon  was  increased,  the  baseline  of  its  infrared  spectrum  shifted  upward  as  a 
result  of  greater  infrared  absorbance. 

This  infrared  absorbance  trend  as  a  function  of  temperature  for  bulk  silicon  can  be 
explained  by  thermal  excitation  of  its  valence  band  electrons  to  the  conduction  band.  At  higher 
temperatures,  the  silicon  became  more  infrared  opaque  as  a  result  of  more  abundant  thermal 
energy  to  initiate  electron  transitions.  The  infrared  absorbance  of  silicon  was  shown  to  be 
directly  proportional  to  the  population  of  mobile  electrons.  The  mobile  electron  population  can 
only  be  increased  by  energy  inputs  sufficient  to  initiate  valence  band  to  conduction  band 
transitions.  This  absorbance  trend  is  analogous  to  the  increasing  spectral  absorbance  of  a 
compound  as  a  function  of  its  concentration. 

Shortly  after  Becker  and  Fan's  work,  H arrick189  conducted  a  series  of  observations  that 
confirmed  infrared  absorbance  by  mobile  electrons.  FI  is  experiments  were  conducted  by  noting 
the  infrared  absorbance  of  a  germanium  bar  with  one  end  heavily  p-doped.  The  bar  was  inserted 
into  a  simple  electric  circuit,  and  infrared  absorbance  was  recorded  for  varying  populations  of 
injected  'holes',  which  increase  the  mobility  of  valence  band  electrons.  H arrick  found  that 
increasing  the  valence  band  electron  mobility  through  hole  injection  registers  a  corresponding 
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increase  in  infrared  absorptivity,  and  conversely  an  extraction  of  positive  charge  carrier  holes 
diminishes  electron  mobility  and  reduces  infrared  absorbance. 

1.5.2  Characteristics  of  the  Ti02  IR  Baseline 

The  behavior  of  a  titanium  dioxide  nanoparticle  infrared  spectrum  baseline  is  the  result  of 
several  interacting  processes.  The  foremost  contributor  to  the  spectrum  is  the  absorbance  exerted 
by  free  charge  carriers,  i.e.  the  conduction  band  electrons.  Secondary  effects  on  the  infrared 
spectrum  result  from  lattice  vibrations  coupled  to  the  conduction  band  electron  absorbances. 

The  atoms  or  ions  that  constitute  a  crystal  lattice  are  not  static  entities.  Each  atom  or  ion 
is  capable  of  limited  movement  and  this  freedom  results  in  lattice  element  oscillations  about 
fixed  positions.  When  an  atom  or  ion  displacement  propagates  through  the  lattice,  the  result  is  a 
quantized  displacement  wave  known  as  a  phonon. 

Two  types  of  distinct  vibrational  wave  motion  are  possible  within  a  crystal  lattice. 
Longitudinal  waves  occur  when  the  motions  of  atoms  or  ions  are  parallel  or  antiparallel  to  the 
direction  of  wave  propagation.  A  transverse  wave  occurs  when  the  direction  of  atom  or  ion 
displacement  is  orthogonal  to  the  direction  of  wave  propagation.  Figure  21  shows  examples  of 
both  kinds  of  wave  displacement. 
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Figure  21:  Vibrational  wave  propagation  in  crystal  lattices.  A  longitudinal  or  compression  wave  is 
shown  in  (a)  and  a  transverse  wave  is  shown  in  (b). 
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Lattice  phonons  are  classified  on  the  basis  of  whether  neighboring  atoms  or  ions  are 
moving  in  phase  (concerted  motion  in  the  same  direction)  or  out-of-phase  (motion  in  opposing 
directions).  When  adjacent  atom  vibrations  are  in  phase,  the  phonon  is  classified  as  acoustic,  and 
when  the  vibrations  are  out-of-phase  the  phonon  is  classified  as  optical.  Acoustic  phonon 
motion  couples  to  conduction  band  electron  infrared  absorbance  and  contributes  to  the  titanium 
dioxide  infrared  absorbance  spectrum  behavior. 

The  solid  line  in  Figure  22  is  an  absorbance  spectrum  recorded  from  titanium  dioxide 
nanoparticles.  This  spectrum  is  representative  of  the  baseline  absorbance  trends  for  titania 
nanoparticles  commonly  observed.  The  monotonic  increase  in  the  baseline  absorbance  from 
3500  to  1200  cm  1  is  caused  by  conduction  band  electron  absorbance  of  infrared  radiation.  The 
conduction  band  is  comprised  of  a  continuum  of  states.  Excitation  of  electrons  into  the  band  and 
transitions  of  electrons  within  the  conduction  band  energy  states  leads  to  a  Boltzmann-like 
electron  energy  distribution,  as  is  illustrated  in  the  Figure  22  inset. 


Figure 22:  Phonon  coupling  to  free  carrier  absorbance  in  titanium  dioxide.  The  inset  diagram190  shows 
the  origin  of  the  conduction  band  (CB)  electron  energy  states  distribution  from  thermal  excitation  out 
of  shallow  trap  stateT  i 3+  ion. 
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The  large  peak  at  1200  cm  1  is  due  to  phonon  absorbance.  At  the  low  energies  associated 
with  radiation  from  750  to  1200  cm"1,  the  titanium  dioxide  lattice  absorbance  blocks  infrared 
radiation  from  the  spectrometer  radiation  source.  From  1200  cm"1  to  higher  frequencies,  two 
trends  shape  the  baseline  absorbance. 

The  infrared  absorbance  of  free  charge  carriers  (i.e.  conduction  band  electrons)  is 
described  by  the  polynomial  equation: 

af  =  AAL5  +  5A2-s  +  CA3-5  (51) 

where  af  is  the  free  charge  carrier  absorbance.50  The  radiation  wavelength  is  given  by  X,  and  the 
variables  A,  B,  and  C  are  constants.  When  the  expression  is  converted  from  wavelength  to 
frequency,  v,  the  equation  becomes: 

af  =  Av'1'5  +  Bv~25  +  C'V3-5  (52) 

A',  B',  and  C'  are  also  constants.  The  first  term  in  Equation  (51)  with  the  1.5  exponential  gives 
the  infrared  absorbance  of  electrons  coupled  to  lattice  acoustic  phonons.  The  second  term  with 
the  2.5  exponential  is  the  conduction  band  electron  optical  phonon  coupling  absorbance.  The 
last  term  with  the  3.5  exponential  is  the  absorbance  for  phonon  coupling  to  ionized  impurities  in 
the  crystal  lattice.  These  term  definitions  for  Equation  (51)  carry  over  to  Equation  (52)  where 
the  exponentials  are  inverted  from  the  wavelength-to-frequency  conversion.  Between 
approximately  3500  to  1600  cm"1  and  1600  to  1200  cm"1  in  the  Figure  22  spectrum,  acoustic 
phonon  coupling  to  the  conduction  band  electrons  dominates  the  absorbance  as  is  indicated  by 
the  slightly  offset  Ao  15  and  Bd  1,5  trend  lines  shown  for  comparison.  The  change  in  phonon 
coupling  trends  at  1650  cm"1  might  be  attributed  to  superposition  of  two  absorbance  spectra.  The 
additional  spectral  trend  may  originate  from  a  thin  layer  of  tungsten  oxide  (WO3)  formed  on  the 
surface  of  the  tungsten  mesh  support  atthetungsten-titanium  dioxide  interface. 
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The  experimental  apparatus  used  to  obtain  the  spectrum  shown  in  Figure  22  uses  an 
electrochemical ly  etched  tungsten  mesh  for  both  mechanical  support  and  resistive  heating  of  the 
titanium  dioxide  nanoparticle  samples.  Details  of  this  apparatus  and  sample  preparation  are 
presented  in  section  2.5.2  .  The  titanium  dioxide  nanoparticles  that  are  observed  with  an  infrared 
spectrometer  are  pressed  into  the  openings  of  the  tungsten  support  mesh,  covering  a  circular 
region  of  approximately  one  centimeter  diameter.  The  nanoparticle-mesh  assembly  is  then 
oriented  normal  to  the  infrared  spectrometer  source  beam  within  the  sample  compartment.  This 
geometry  causes  infrared  radiation  from  the  instrument  to  enter  numerous  pores  in  the  mesh 
where  nanoparticles  are  packed.  Oxidation  of  the  tungsten  support  might  occur  either  from 
gaseous  oxygen  introduced  into  the  system  during  oxidative  sample  cleaning  procedures,  or  from 
abstraction  of  lattice  oxygen  directly  from  the  crystalline  titanium  dioxide  at  the  nanoparticle- 
tungsten  interface  within  the  mesh  openings  as  shown  in  Figure  23.  Tungsten  trioxide  is  an 
indirect  semiconductor191  with  a  band  gap  energy  similar  to  titanium  dioxide  (3.10  eV)192.  The 
tungsten  trioxide  present  as  a  thin  coating  on  the  support  mesh  within  the  sample  orifices  might 
contribute  a  weak  absorbance  signal  as  shown  in  Figure  23a  This  additional  signal,  combined 
with  the  titanium  dioxide  absorbance  spectrum  in  Figure  23b  could  produce  a  composite 
spectrum  as  shown  in  Figure  23b  This  combined  spectrum  would  have  two  acoustic  phonon 
couplings,  and  thus  two  separate  0  15  absorbance  trends. 
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Figure  23:  Superposition  of  interfacial  W03  infrared  spectrum  (a)  with  Ti02  sample  spectrum  (b)  to 
produce  observed  composite  infrared  spectrum  (c). 


Perturbation  of  the  infrared  spectral  baseline  signal  of  titanium  dioxide  nanoparticles 
presents  a  sensitive  measure  of  the  electron  density  within  that  material's  conduction  band.  The 
extent  of  thermal  excitation  of  electrons  to  the  titania  conduction  band  for  a  given  sample 
temperature  can  be  a  relative  indication  of  the  lattice  defect  density  on  the  nanoparticles  where 
electron  trapping  occurs.  Electron  trapping  sites  at  the  edge  of  the  conduction  band  energy  on 
the  nanoparticles  are  the  means  for  thermal  catalysis  of  adsorbate  reactions.  Although  multiple 
methods  are  available  to  characterize  the  energy  difference  between  these  shallow  electron  traps 
and  the  conduction  band,  infrared  spectroscopy  offers  possibly  the  simplest  approach.  M  y  goal 
is  to  demonstrate  a  method  for  accurately  characterizing  that  electron  shallow  trap  state  energy 
difference  for  titanium  dioxide  nanoparticles  from  infrared  spectroscopic  data  recorded  over 
multiple  sample  temperatures. 
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1.5.3  E  lectron  Density  Measurement  from  Infrared  Absorption 

Electrons  belong  to  the  class  of  particles  known  as  fermions.  This  class  of  particles  is 
characterized  by  half-integral  'spin'  (i.e.  spin  =  ±  1/2,  ±  3/2,  ±  5/2,  ... )  and  includes  electrons, 
protons,  3FI  nuclei,  and  13C  nuclei.  Particles  that  have  integral  spin  (i.e.  spin  =  ±  1,  ±  2,  ±  3,  ... ) 
are  bosons,  and  include  photons,  2H  nuclei,  and  12C  nuclei.  Use  of  the  word  'spin'  should  not  be 
taken  as  literal,  rotational  motion  of  the  electron.  The  angular  momentum  characteristic  of  the 
electron,  discovered  by  Goudsmidt  and  Uhlenbeck193  is  a  parameter  introduced  by  Dirac  in  his 
modification  of  the  Schrodi nger  equation  for  an  electron  to  account  for  the  relativistic  effects  of 
electron  motion.40 

Asa  fermion,  the  probability  of  electron  occupancy  of  an  energy  state  is  described  by  the 
Fermi  probability,  F(E)\ 


F{E)  = 


1  +  Exp 


E-E 


Fermi 


KT 


(53) 


where  kB  is  Boltzmann's  constant  (1.381  x  10"23  J  •  K  1),  T  is  the  thermodynamic  temperature, 
and  Efermi  is  the  Fermi  energy  for  that  material.194  A  derivation  of  the  Fermi  probability 
expression  is  given  in  Appendix  1.  Figure  24a  shows  the  Fermi  probability  for  a  semiconductor 
at  a  nonzero  temperature.  The  top  of  the  valence  band  in  Figure  24  is  set  as  the  zero  energy,  and 
the  bottom  edge  of  the  conduction  band  is  separated  from  the  valence  band  by  band  gap  energy, 
Egap.  Note  from  the  Figure  24a  magnified  inset,  that  at  T  >0K  the  Fermi  probability  of  electron 
occupancy  of  the  conduction  band  is  not  zero. 
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Figure  24:  Development  of  conduction  band  electron  density  expression.  The  dark  shaded  regions  are 
integrated  areas.™ 


The  conduction  band  of  a  semiconductor  crystal  lattice  is  composed  of  a  continuum  of 
discrete  energy  states.  The  number  of  states  with  energy  E  that  can  be  occupied  within  a  crystal 
lattice  is  given  by  the  density  of  states,  p(E)\ 


P(E)  = 


(54) 


where  me  is  the  electron  rest  mass  (9.1094  x  10"31  kg)  and  h  is  Planck's  constant  (6.626  x  10'34  J 
•  s).194  A  derivation  for  the  density  of  states  expression  is  given  in  Appendix  2.  Figure  24b 
shows  the  number  of  available  conduction  band  energy  states  as  a  dark  shaded  region. 

The  density  of  electrons  N(E)  between  an  energy  state  E  and  E  +  dE  in  a  crystal  lattice  is 
the  product  of  the  density  p(E)  of  those  energy  states  in  the  lattice  and  the  Fermi  probability  that 
those  states  will  be  occupied: 

N(E)dE  =p(E)F(E)clE  (55) 


where  dE  is  an  infinitesimal  increment  of  energy  E.  Substituting  for p(E)  and  N(E)\ 
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N(E)clE  = 


1  +  Exp 


E-E 


-clE 


Fermi 


KT 


(56) 


The  thermal  energy  available  at  a  temperature  T  is  given  by  kBT.  For  the  studies 
presented  here  the  nanoparticle  annealing  experiments  were  run  from  an  ambient  temperature  of 
298  K  to  a  final  temperature  of  675  K.  The  energy  available  at  675  K  is  approximately  0.058 
eV.  With  £  at  the  bottom  of  the  conduction  band  for  titanium  dioxide,  Egap  -  EFermi  is 
approximately  1.5  eV,  which  is  30  times  greater  than  kBT.  At  the  lower  temperatures,  this  energy 
disparity  is  even  greater,  so  an  approximation  can  be  made  for  the  Fermi  probability: 
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The  density  of  electrons  contained  within  the  conduction  band  can  be  determined  through 
integration  of  Equation  (58)  from  bottom  energy  of  the  conduction  band  (Egap)  to  top  of  the 
conduction  band.  The  top  of  the  band  can  be  taken  as  £  =  oo.  The  density  of  conduction  band 
electrons,  NCb,  is  calculated  from: 


Ncb(E)=  J  p(E)F(E)dE  (59) 

Er 

Gap 
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Ncb(E)  =  AT3'2  Exp 
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(60) 
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The  dark  shaded  region  in  Figure  24c  is  proportional  to  the  density  of  conduction  band  electrons 
at  temperature  7. 

The  focus  of  my  thesis  is  to  determine  the  energy  of  transition  from  a  shallow  trap  state, 
just  below  the  conduction  band,  into  the  conduction  band.  The  thermal  excitation  is  from  Estate  to 
the  conduction  band,  and  not  across  EGap .  A  n  appropriate  substitution  is: 


N(E)=BTmExp 


-CE 


State 


>5/2 


where  B  ■ 


(meKkB) 

h3 


and  C  : 


2  kr 


(61) 


The  infrared  baseline  absorbance  for  titanium  dioxide  nanoparticles  at  temperature  T  is 
proportional  to  the  number  of  mobile  electrons  populating  the  conduction  band.  The  relationship 
between  the  spectral  absorbance  and  concentration  of  an  absorbing  species  is  given  by  the  Beer- 
LambertLaw: 


A  =ebc  (62) 

where  A  is  the  absorbance,  s  is  the  molar  absorptivity,  bis  the  absorbing  media  path  length,  and  c 
is  the  absorbing  species  concentration.  The  molar  absorptivity  term  of  the  Beer-Lambert  Law  is 
a  constant  that  is  wavelength  and  material  specific.  The  absorption  path  length  for  a  solid 
sample  can  also  be  considered  constant.  Integration  of  the  absorbances  over  a  range  of 
wavelengths  (or  wavenumber  frequencies)  for  a  given  concentration  of  absorbing  species  will 
give  an  area  that  is  proportional  to  the  concentration  of  those  absorbing  species.  The 
proportionality  constant  for  this  area-concentration  relation  is  the  integrated  trend  of  the 
absorptivity  constants  as  they  vary  with  wavelength. 

The  population  of  electrons  in  the  conduction  band  of  a  semiconductor  absorbs  infrared 
radiation  over  a  broad  range  of  wavelengths  through  intra-band  electron  transitions.190  Once 
excited  into  the  conduction  band,  a  population  of  electrons  will  energy  partition  across  the  range 
of  available  energies  (0.45  eV  for  the  4000  to  400  cm"1  infrared  spectrum)  in  a  Boltzmann-like 
distribution.  It  is  this  intra-band  transition  over  the  continuum  of  available  electron  states  in  the 
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conduction  band  that  produces  the  characteristic  baseline  absorbance  profile  for  the 
semiconductor.  As  more  electrons  enter  the  conduction  band  from  increased  thermal  excitation, 
more  electrons  are  present  to  absorb  available  energy  and  the  baseline  absorbance  increases.  The 
population  of  electrons  observed  under  constant  conditions  within  a  fixed  quantity  of 
semiconductor  material  is  an  electron  concentration.  It  follows  from  the  Beer-Lambert  Law  that 
an  integrated  area  Z  under  an  infrared  spectrum  baseline  is  also  proportional  to  the  conduction 
band  electron  population: 


Z  °c  BTmExp 


-CE 


State 


->5/2 
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2  kr 


(63) 


If  a  proportionality  constant  D  is  introduced  to  Equation  (63)  a  new  equation  is  created  that  can 
be  solved  for  the  semiconductor  shallow  trap  state  energy,  Estate'. 
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A  set  of  transmission  infrared  spectra  can  be  recorded  for  a  sample  of  titanium  dioxide 
nanoparticles  held  at  varying  constant  temperature  setpoints.  For  each  spectrum,  a  part  of  the 
baseline  is  integrated.  The  baseline  area  vs.  nanoparticle  temperature  data  set  are  fit  to  Equation 
(64)  using  nonlinear  regression  techniques  where  D  and  Estate  are  parameters  of  the  regression. 
In  this  manner,  the  shallow  trap  state  energy  Estate  can  be  experimentally  measured  with  infrared 
spectroscopy. 
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Chapter  2:  Nanopartide Preparation,  Characterization,  and  Spectroscopy 


Z1  Introduction 

My  work  with  titanium  dioxide  nanoparticles  began  with  synthesis  of  nanoparticles  by 
laser  evaporation  techniques.  This  method  of  producing  nanoparticles  was  chosen  on  the  basis 
of  an  existing  body  of  experience  by  previous  researchers  in  our  group,  the  greater  likelihood  of 
synthesizing  nanoparticles  free  of  contamination,  and  the  simplicity  of  the  preparation  technique. 
Laser  evaporation  produces  lower  yields  of  titanium  dioxide  nanoparticles  than  solution  based  or 
flame  pyrolysis  techniques,  but  the  resulting  particles  require  no  post  synthesis  calcination  to 
induce  crystallinity 

Characterization  of  the  nanoparticles  principally  consisted  of  examination  by 
transmission  electron  microscopy.  X-ray  diffraction  studies  and  high-resolution  transmission 
electron  microscopy  was  performed  on  a  limited  set  of  samples  to  confirm  crystallinity  and  to 
determine  which  polymorphs  (rutile  and  anatase)  were  present  within  the  nanoparticles. 

Nanoparticle  samples  were  pressed  onto  a  tungsten  mesh  and  mounted  for  infrared 
spectroscopic  examination  under  high  vacuum  conditions.  Spectra  were  acquired  under 
temperature  controlled  conditions,  where  annealing  was  accomplished  by  resistive  heating  of  the 
tungsten  support  mesh.  The  resulting  spectra  were  saved  as  text  files,  which  were  latter 
processed  with  custom  FORTRAN  77  programs  and  M  icrosoft  Excel®  spreadsheets. 

22  Nanopartide  Preparation  by  Laser  Vaporization 
ZZ1  Description  of  Apparatus 

Nanoparticle  preparation  by  laser  vaporization  was  accomplished  by  use  of  the  apparatus 
designed  and  constructed  by  Professor  Brian  Tissue195,  196  and  later  used  by  Dr.  Wesley 
Gordon.197  199  The  chamber  itself  is  cylindrical,  externally  measuring  approximately  54  cm  in 
length  and  20  cm  in  diameter,  and  constructed  from  stainless  steel.  The  gas  pressure  within  the 
chamber  is  controlled  by  pumping  with  an  Alcatel  10  liter  rotary  vane  pump  and  pressurization 
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from  a  regulated,  ultra-high  purity  grade  oxygen  cylinder.  Both  vacuum  and  gas  ingress  to  the 
chamber  are  controlled  by  user  adjustable  valves  located  at  the  chamber.  The  pressure  within  the 
chamber  was  monitored  with  a  Granville-Phillips  Model  275  Convectron  vacuum  pressure 
sensor.  The  chamber  pressure  signal  from  the  Convectron  sensor  was  read  on  an  MKS 
Instruments  Baratron  type  107A  analog  display.  The  oxygen  pressure  within  the  chamber  was 
estimated  from  calibration  data  supplied  with  the  Convectron  sensor.  Figure  25  shows  a 
schematic  of  the  vaporization  chamber.  The  vaporization  laser  beam  was  focused  with  a  zinc 
selenide  lens  external  and  adjacent  to  the  chamber.  The  chamber  side  entry  window  for  the 
vaporization  beam  was  also  formed  from  zinc  selenide  for  infrared  radiation  transparency. 


Figure  25:  Schematic  of  laser  evaporation  chamber  used  for  nanoparticie  synthesis  (end  view).  Adapted  from  Tissue.195 


The  laser  used  for  vaporization  was  a  Syn  Rad,  M  odel  648-5-28W-1344  carbon  dioxide 
laser.  The  maximum  power  rating  of  the  laser  was  100  Watts.  Control  of  the  laser  was  through 
a  Syn  Rad  UC-1000  module,  which  featured  pulsed  and  continuous  operating  modes  and  a 
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continuously  variable  power  adjustment.  The  laser  was  water  cooled  with  a  recirculating  water 
chiller.  The  laser  beam  was  directed  to  the  entrance  optics  of  the  vaporization  chamber  by  an 
optical  system  diagrammed  in  Figure  26.  The  carbon  dioxide  laser  beam  passed  through  a  1:100 
ratio  beam  splitter  with  the  1%  reflected  beam  directed  to  a  M  elles  Griot  optical  power  sensor 
head,  rated  at  30  Watts  maximum  power.  The  digital  readout  display  that  accompanied  the 
sensor  head  was  set  to  read  units  of  Watts.  The  effective  vaporization  laser  power  level  was 
estimated  by  multiplication  of  the  displayed  power  reading  by  a  factor  of  100X.  Vaporization 
experiments  were  typically  run  at  100%  power  settings,  which  resulted  in  approximately  60 
Watts  laser  power. 
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Figure  26:  Vaporization  chamber  laser  optics  (top  view). 


Working  in  the  presence  of  a  60  Watt  laser  beam,  capable  of  heating  ceramics  to 
incandescence,  makes  heightened  observance  of  safety  precautions  necessary.  During  operation 
of  the  vaporization  laser,  access  to  areas  of  the  optical  bench  where  the  beam  passes  are  blocked 
by  two  aluminum  stop  plates  mounted  on  the  front  of  the  laser  optical  table.  Prior  to  any 
operation  of  the  carbon  dioxide  laser,  a  secondary  low-power  helium-neon  laser,  aligned  coaxial 
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to  the  primary  vaporization  beam,  is  used  to  check  the  laser  optics  alignment.  The  greatest 
likelihood  of  beam  misalignment  is  at  the  vertical  beam  steering  mirror  assembly  adjacent  to  the 
vaporization  chamber.  Here,  adjustment  of  the  mirrors  during  laser  operation  to  move  the 
vaporization  locale  on  the  target  pellet  can  lead  to  stray  reflections  from  the  mirrors,  zinc 
selenide  focusing  lens,  and/or  the  zinc  selenide  entrance  window  of  the  chamber.  Prior  to 
operation  of  the  primary  vaporization  laser,  checks  are  made  in  the  vicinity  of  the  chamber 
entrance  port  for  unexpected  He-Ne  reflections  from  the  optics  system  and  alignment 
adjustments  are  made  as  needed  to  ensure  that  all  laser  beams  are  coaxial  and  directed  only  into 
the  vaporization  chamber. 

A  secondary  purpose  for  the  helium-neon  laser  is  to  check  the  path  of  the  vaporization 
laser  beam  inside  of  the  chamber.  Prior  to  sealing  and  evacuating  the  chamber,  the  beam 
placement  onto  the  appropriate  region  of  the  starting  material  pellet  must  be  verified.  The 
uninterrupted  passage  of  the  laser  beam  through  the  entrance  orifice  cut  into  the  nanoparticle 
collector  target  surface  also  must  be  confirmed.  Any  irregularities  in  laser  beam  or  collector 
surface  alignment  are  corrected  before  operation  of  the  vaporization  laser. 

Nanoparticle  condensates  are  collected  within  the  vacuum  chamber  on  a  stainless  steel 
surface.  The  preferred  collector  geometry  was  a  hemispherical  dome  arrange  symmetrically 
above  the  starting  material  target  pellet.  Both  the  interior  surfaces  of  the  chamber  and  the 
collector  hemisphere  assembly  were  cleaned  prior  to  use  with  95%  ethyl  alcohol  and  laboratory 
wipe  tissues.  The  collector  hemisphere  was  supported  upon  an  adjustable  metal  bracket  and 
placed  so  that  its  lower  edge  was  even  with  the  top  of  the  target  pellet  surface.  Proper  alignment 
of  the  collector  dome  allowed  observation  of  the  target  pellet  surface  from  a  forward  view  port 
and  unencumbered  passage  of  the  vaporization  laser  beam  through  another  port,  90  degrees  from 
the  view  port.  Details  of  the  hemisphere  placement  are  given  in  Figure  27. 
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collector  support  bracket 


a:  collector  hemisphere  top  view  b:  vaporization  chamber  front  view 

Figure  27:Details of  hemispherical  collector  dome  placement. 


2.2.2  Nanoparticle  Synthesis 

The  starting  material  for  all  nanoparticle  syntheses  experiments  was  a  partially  reduced 
titanium  dioxide  pellet,  part  number  T -2051  purchased  from  CERAC,  Inc.  These  pellets  were 
assayed  at  99.9%  titanium  oxide,  rutile  form.  The  partially  reduced  form  of  the  titania  (non- 
stoichiometric  TiOx,  x  <  2)  imparted  a  dark  blue-black  coloration  on  the  pellet  surfaces.  Before 
and  after  each  vaporization  synthesis,  the  target  pellets  were  weighed  to  estimate  the  amount  of 
material  lost  from  the  pellet  during  vaporization. 

Following  setup  of  the  vaporization  chamber  and  satisfactory  laser  optics  alignment,  the 
vaporization  chamber  was  sealed  and  evacuated  with  the  rotary  vane  pump  to  a  pressure  of 
approximately  10  T orr  on  the  Convectron  gauge.  The  vacuum  valve  was  closed  and  the  chamber 
was  backfilled  with  ultra  high  purity  (UHP)  oxygen  to  a  pressure  of  approximately  100  Torr. 
The  chamber  was  then  re-evacuated  and  backfilled  an  additional  two  times  to  remove  residues  of 
air  from  the  apparatus.  Following  the  third  U  FI  P  oxygen  flush,  the  chamber  was  evacuated  and 
backfilled  with  oxygen  to  the  required  operating  pressure  for  the  nanoparticle  synthesis 
experiment.  Upon  stabilization  of  the  vaporization  chamber  to  its  operating  pressure,  the 
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vaporization  laser  was  switched  on  and  set  to  a  power  of  approximately  30%.  The  target  pellet 
was  slowly  rotated  with  the  30%  beam  impinging  its  surface  for  approximately  two  minutes  to 
gradually  anneal  the  surface  of  the  pellet.  Following  the  annealing  operation,  the  laser  power 
was  increased  to  100%  and  the  start  time  for  the  vaporization  run,  in  addition  to  the  chamber 
pressure  and  laser  power,  were  recorded. 

During  the  laser  vaporization  interval,  the  chamber  pressure  was  maintained  by  manual 
adjustment  with  introduction  of  additional  oxygen  and  evacuation  through  the  appropriate  gas 
inlet  and  vacuum  outlet  control  valves.  The  target  pellet  was  manually  rotated  to  ensure  that  no 
single  region  of  its  surface  was  over  eroded  by  the  laser  beam.  When  all  regions  for  a  particular 
radius  of  the  target  pellet  were  eroded,  the  vaporization  beam  was  relocated  to  an  uneroded 
radius  by  adjustment  of  the  laser  beam  steering  mirrors  external  to  the  chamber. 

Vaporization  synthesis  of  titanium  dioxide  nanoparticles  typically  was  conducted  for  one 
hour  to  ninety  minutes  of  laser  operation,  dependant  on  the  condition  of  the  target  pellet  as  the 
vaporization  progressed.  At  the  end  of  target  vaporization,  the  carbon  dioxide  laser  power  was 
gradually  reduced  to  zero  percent  and  the  laser  was  shutdown.  The  oxygen  input  to  the  chamber 
for  pressure  maintenance  was  shut  off  and  the  vacuum  metering  valve  was  opened  full  to 
evacuate  the  chamber  of  any  residual  titanium  dioxide  vapor  or  suspended  particles.  Following 
approximately  five  minutes  of  evacuation,  the  vacuum  metering  valve  was  closed  and  the  pump 
system  powered  down.  The  U  H P  oxygen  was  then  opened  to  increase  the  chamber  pressure  to 
atmospheric  so  that  the  chamber  front  access  door  could  be  opened. 

Following  a  sufficient  time  to  allow  the  collector  hemisphere  dome  to  cool,  the 
vaporization  chamber  front  access  door  was  fully  opened  and  the  collector  removed  to  a  clean 
work  surface.  The  nanoparticle  film  was  scraped  from  the  interior  surface  of  the  collector  dome 
with  a  plastic  laboratory  spatula  and  transferred  to  a  tared  shell  vial.  The  masses  of  the  vial  + 
particles,  in  addition  to  the  post  vaporization  mass  of  the  target  pellet  were  measured  on  a 
laboratory  balance  and  recorded.  Post  synthesis  tasks  included  cleanup  of  the  vaporization 
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chamber  with  ethanol  and  laboratory  tissue  wipes,  shutdown  of  the  UH  P  oxygen  supply  cylinder, 
and  replacement  of  protective  covers  over  the  vaporization  laser  optics. 

An  initial  concern  of  the  nanoparticle  laser  syntheses  was  effective  control  of  mean 
particle  diameter  and  the  range  of  particle  diameters  produced  for  a  set  of  operating  conditions. 
Two  approaches  to  affecting  particle  size  control  that  were  investigated.  One  technique  was 
variation  of  the  target-to-pel let  distance  r,  (see  Figure  27b)  by  using  collector  hemispheres  of 
varying  radii.  The  second  and  more  successful  approach  investigated  was  nanoparticle  synthesis 
at  different  oxygen  pressures.  The  results  of  this  work  are  summarized  in  the  next  chapter. 

2.3 Transmission  Electron  Microscopy 
2.3.1 1  mage  Resolution 

In  1873  Ernest  Abbe  described  the  limit  of  optical  microscopy  image  resolution  as 
approximately  200  nm  (2  x  10"7  m).200  Two  structures  that  are  physically  apart  by  greater  than 
approximately  one  half  of  the  lower  wavelength  of  visible  light,  400  nm,  cannot  be  distinguished 
with  conventional  optical  microscopy.  The  resolving  power  gis  defined  as: 

8  =0.61 — ^ —  (67) 

n  sin  a 

where  2  is  the  wavelength  of  the  light,  n  is  the  refractive  index  of  the  medium  (n=l  for  vacuum), 
and  a  is  the  lens  acceptance  half  angle.201  The  0.61  factor  originates  from  the  observation  that 
only  part  of  the  first  diffraction  maximum  is  necessary  to  form  a  resolved  image.202 

M  agnification  that  exceeds  the  Abbe  theoretical  limit  became  possible  with  the  discovery 
of  the  electron  and  its  wave-particle  duality.  Davisson  and  Germer  had  observed  Bragg 
diffraction  of  an  electron  beam  by  a  crystal  of  nickel  in  1927  which  strongly  suggested  that 
electrons  could  be  dealt  with  as  electromagnetic  radiation.203, 204  Over  1925  to  1927  FI  ans  B  usch 
discovered  that  a  magnetic  field  could  be  used  as  an  electron  lens  in  a  similar  manner  to  the 
function  of  glass  lenses  for  light.  By  the  close  of  the  1930' s  the  light  based  optical  microscopy 
resolution  limit  had  been  exceeded  by  researchers  conducting  magnification  with  electrons 
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directed  in  either  magnetic  or  electrostatic  electron  'optics'.202  Electron  microscopy  had  become 
a  viable  research  tool. 

The  de  Broglie  wavelength  /l  for  an  electron  is  calculated  from  the  equation: 

h  h 

A=-=—  (68) 

P  mev 

where  h  is  the  Planck  constant,  p  is  the  electron  momentum,  me  is  the  electron  rest  mass,  and  v  is 
the  electron  velocity.90  The  velocity  of  an  electron  is  determined  from  a  classical  mechanics 
treatment  of  its  acceleration  within  an  electric  field.  At  the  instant  of  an  electron's  departure 
from  a  cathode  in  a  vacuum,  all  of  the  particle's  energy  is  potential.  This  potential  energy  is  the 
direct  result  of  the  electron's  repulsion  from  the  negatively  charged  cathode  and  its  attraction  to 
the  opposing  positively  charged  anode.  As  the  electron  travels  to  the  anode,  it  trades  potential 
energy  for  kinetic  energy.  By  the  time  that  the  electron  reaches  the  anode  its  energy  is 
completely  kinetic  and  given  by: 

Ue=—mv 2  (69) 

2  e 

where  U  is  the  electric  field  potential  in  volts  and  e  is  the  elementary  charge,  previously  given. 
Equation  (69)  can  be  rearranged  to  give  the  velocity  of  an  electron  as  a  function  of  the  voltage 
potential  through  which  it  was  accelerated: 


(70) 


Substitution  of  Equation  (70)  into  Equation  (68)  makes  the  calculation  of  electron  wavelength  as 
a  function  of  acceleration  field  voltage  possible.  Examples  of  accelerating  voltages,  their 
resulting  electron  wavelengths,  and  the  corresponding  image  resolutions  for  a  typical 
transmission  electron  microscope  are  given  in  T able  2: 
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Voltage  (kV) 

Wavelength  (nm) 

Resolution  (nm) 

20 

0.0087 

0.44 

40 

0.0061 

0.31 

60 

0.0050 

0.25 

80 

0.0043 

0.21 

100 

0.0039 

0.19 

1000 

0.00087 

0.10 

Table  2:  Comparison  of  accelerating  voltage,  electron  wavelength,  and  resolving  power  for  a 
transmission  electron  microscope.201 


2.3.2  E  lectron  M  icroscope  0  peration  M  odes 

Conventional  transmission  electron  microscopes  can  operate  in  two  basic  modes.  If  the 
electrons  impinge  upon  the  lattice  planes  of  the  specimen  meet  the  Bragg  criteria,  diffraction 
occurs.  Alternately,  the  electrons  can  pass  through  the  sample  undeflected  from  their  original 
paths.  Through  selection  of  apertures  and  magnetic  lens  settings,  either  the  transmitted  electrons 
can  be  focused  to  form  a  magnified  image  (bright  field  mode)  or  the  diffracted  electrons  can  be 
focused  to  form  a  projected  diffraction  pattern  (dark  field  mode).205  The  diffraction  pattern  of 
the  specimen  forms  at  the  back  focal  plane  of  the  microscope  objective  lens  and  is  intrinsic  to  the 
operation  of  the  instrument,  whether  the  application  accesses  the  diffraction  information  or  the 
bright  field  magnified  image.  The  ray-trace  diagrams  of  both  modes  are  shown  in  Figure  28  for 
a  single  objective  lens  and  schematic  representations  of  the  electron  optics  for  a  transmission 
electron  microscope  operating  in  both  image  magnification  and  electron  diffraction  modes  are 
given  in  Figure 29. 
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Figure  28:  Electron  diffraction  in  a  transmission  electron  microscope.205  Broken  lines  indicate  the 
diffracted  electron  beam  and  solid  lines  correspond  to  the  transmitted  electron  beam.  The  'lens’  in  the 
figure  is  actually  a  magnetic  field. 


Figure  28  shows  the  possibilities  of  either  focusing  transmitted  electrons  on  the  optical 
axis  for  a  bright  field  image  or  off-axis  imaging  for  a  dark-field  diffraction  image.  In  both  cases, 
a  single  focus  point  on  the  objective  lens  back  focal  plane  is  used  to  form  an  image.  Another 
possibility  exists  where  a  much  larger  part  of  the  back  focal  plane  information  is  utilized. 
Through  manipulation  of  the  magnetic  lens  settings  and  Fourier  transform  processing  of  back 
focal  plane  information  it  is  possible  to  greatly  increase  the  resolution  of  the  magnified  image. 
High  resolution  transmission  electron  microscopy  images  extend  magnification  of  the  specimen 
to  the  extent  that  the  patterns  of  order  on  the  atomic  scale  can  be  discerned.  These  patterns  are 
actually  interference  patterns  between  the  microscope  electron  beam  and  the  electron  density  of 
the  specimen's  atomic  structure.  With  high-resolution  imaging  techniques,  detail  such  as  lattice 
defects  and  crystal  plane  orientations  are  observable. 
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Figure  29:  Diagram  for  electron  optics  of  a  transmission  electron  microscope  operated  in  (a) 
magnification,  and  (b)  electron  diffraction  modes.  Note  that  the  magnetic  lenses  are  indicated  by 
cross-hatched  boxes.202 


2.3.3  Nanoparticle  Preparation  for  Electron  M  icroscopy 

My  transmission  electron  microscopy  work  on  titanium  dioxide  nanoparticles 
encompassed  all  of  the  imaging  techniques  described  above.  I  utilized  the  Zeiss  Model  10CA 
located  at  the  Virginia-M  aryland  Regional  College  of  Veterinary  Medicine  for  diameter 
measurement  of  nanoparticles  and  observation  of  aggregation  and  particle  morphology.  The 
Philips  EM  420  microscope  presently  located  at  the  Institute  for  Critical  Technology  and 
Applied  Science  (ICTAS)  was  briefly  used  to  observe  electron  diffraction  patterns  for  some 
selected  samples.  Also  at  ICTAS,  the  FEI  Titan  300  microscope  was  used  for  high  resolution 
imaging  of  some  selected  nanoparticle  preparations.  The  data  collected  from  both  the  Philips 
EM  420  and  the  FEI  Titan  300  were  used  to  observe  the  crystallinity  of  the  nanoparticle 
preparations  examined. 

The  most  common  electron  microscopy  performed  was  size  measurement  imaging  on  the 
Zeiss  10CA  microscope.  Samples  that  were  examined  on  this  instrument  were  placed  on 
electron  microscopy  grids.  The  grids  that  I  used  were  a  #150  copper  mesh  coated  with  a  varnish 
support  film  layer,  the  Electron  M  icroscopy  Science's  model  FCF150-Cu  grid. 
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Sample  preparation  for  transmission  electron  microscopy  was  accomplished  by 
dispersion  of  a  small  quantity  of  nanoparticles  (<0.001  g)  into  reagent  grade  100%  ethyl  alcohol. 
The  dispersions  were  prepared  in  unused  7  mL  screw-cap  shell  vials  that  were  rinsed  with  two 
aliquots  of  the  dispersion  solvent  immediately  prior  to  use. 

Following  their  solvent  rinses,  the  sample  dispersion  vials  were  completely  filled  with 
additional  dispersion  solvent  and  immediately  capped  pending  sample  introduction. 
Nanoparticle  samples  were  transferred  to  the  alcohol  filled  shell  vials  with  unused,  disposable 
glass  Pasteur  pipettes  that  had  been  alcohol  rinsed  immediately  prior  to  use.  The  rinsed  pipette 
was  immediately  dried  with  a  stream  of  ultra  high  purity  nitrogen  and  the  pipette  tip  was  then 
used  to  pick  up  a  small  quantity  of  nanoparticles.  The  particle  loaded  pipette  tip  loaded  was  then 
immediately  immersed  into  a  dispersion  solvent  filled  shell  vial,  at  which  time  the  nanoparticles 
were  rinsed  from  the  pipette  into  the  solvent.  Following  withdrawal  of  the  glass  pipette,  the  shell 
vial  was  capped  and  agitated  in  an  ultrasonic  bath  for  approximately  thirty  minutes  to  disperse 
the  transferred  nanoparticles. 

Following  ultrasonic  dispersal  of  the  nanoparticles  into  the  ethanol,  a  miniature 
disposable  polyethylene  transfer  pipette  was  used  to  place  one  drop  of  the  alcohol  nanoparticle 
dispersion  onto  an  unused  microscopy  sample  grid.  The  dispersion  alcohol  was  evaporated  from 
the  sample  grid  within  a  covered  Petri  dish  at  room  temperature.  Following  complete 
evaporation  of  dispersal  solvent,  the  prepared  microscopy  grids  were  transferred  to  a  storage  box 
pending  examination  at  one  of  the  aforementioned  electron  microscopes. 

24  X-ray  Diffraction 
24>lThe  Bragg  Condition 

The  diffraction  of  x-rays  by  crystalline  material  was  discovered  by  Max  von  Laue  in 
1912  and  applied  to  structure  determination  by  W.H .  and  W.L.  Bragg  (father  and  son)  in  1913. 206 
The  relationship  between  radiation  wavelength,  angle  of  incidence  to  the  surface  of  a  crystal,  and 
the  spacing  of  atomic  planes  within  the  lattice  is  given  by  Bragg's  law: 
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n  A  =2  d  sin0 


(71) 


The  parameter  n  is  the  order  of  reflection,  or  the  path  difference  in  number  of  wavelengths  of 
radiation  scattered  by  adjacent  planes  of  lattice  atoms.  The  scattered  radiation  wavelength  is  X, 
and  the  distance  between  adjacent  planes  of  radiation  scattering  lattice  atoms  is  the  parameter  d. 
The  angle  6,  or  the  Bragg  angle,  describes  the  angle  of  refraction  where  radiation  interacting 
with  adjacent  planes  of  lattice  atoms,  constructively  interferes.  Collectively,  these  parameters 
are  referred  to  as  the  'Bragg  condition'.  The  interrelation  of  these  parameters  is  illustrated  in 
Figure  30.  The  superficial  appearance  of  the  relationship  between  incident  and  diffracted 
radiation  in  Figure  30  strongly  suggests  that  this  phenomenon  is  similar  to  the  reflection  of  light, 
although  that  depiction  is  an  oversimplification.  Flowever,  the  similarity  of  Bragg  diffraction  to 
radiation  reflection  has  become  part  of  the  terminology  of  this  discipline  and  diffraction  maxima 
are  still  referred  to  as  'reflections'  in  the  literature. 


2.4.2  Production  of  X-rays 

The  radiation  used  to  perform  x-ray  diffraction  is  generated  by  high  energy  electron 
bombardment  of  a  metal  target  under  vacuum  conditions.  Thermionic  emission  from  a  filament 
within  the  x-ray  tube  generates  free  electrons.  Biasing  of  the  filament  potential  with  respect  to 
the  target  anode  to  a  voltage  potential  of  approximately  104  volts  accelerates  the  electrons  to  high 
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kinetic  energies.  These  electron  energies  are  sufficient  to  remove  core  (inner  shell)  electrons 
from  within  the  atoms  of  the  target  metal.  T ransitions  of  electrons  falling  into  the  newly  formed 
electron  vacancies  define  the  radiation  wavelength.  Figure  31  shows  some  of  the  commonly 
encountered  transitions  used  in  the  description  of  x-ray  radiation  wavelengths. 


Figure  31:  Electron  transitions  associated  with  some  x-ray  radiation  sources.206 


The  presence  of  subshells  in  atomic  structure  leads  to  multiplicities  of  the  electron 
transitions  shown  in  Figure  31.  For  example,  there  might  simultaneously  exist  a  Kai,  Ka2,  and 
Kp  transition  for  a  given  metal.  The  wavelengths  given  for  most  x-ray  sources  used  in 
diffraction  work  are  usually  weighted  averages. 

2.4.3  Powder  Diffraction 

The  technique  used  for  x-ray  characterization  of  nanoparticles  is  powder  diffraction.  The 
x-ray  beam  is  diffracted  from  multiple  crystals  at  random  orientations  within  the  sample  mount. 
The  diffracted  beams  strike  a  detector  plane,  either  photographic  film  or  the  path  of  a  moving 
electronic  detector.  The  random  orientation  of  diffracting  crystals  causes  not  one  diffraction  spot 
for  a  characteristic  set  of  Bragg  conditions,  but  many.  The  multiplicity  of  diffraction  spots  for  a 
unique  set  of  conditions  causes  the  formation  of  a  Debeye-Scherrer  diffraction  ring,  as  shown  in 
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Figure  32.  A  plot  of  the  diffraction  ring  intensities  against  their  associated  20  values  produces  a 
powder  x-ray  diffraction  pattern. 


Figure 32:  Powder  x-ray  diffraction. 


2.4.4  Sample  Preparation  and  Diffraction  Pattern  Collection 

My  x-ray  diffraction  characterization  of  titanium  dioxide  nanoparticles  was  conducted 
upon  theOxford  Diffraction,  Gemini  diffractometer  located  in  the  Virginia  Tech  crystallography 
laboratory.  Powder  diffraction  are  not  commonly  performed  on  this  instrument,  as  its  main  work 
is  performed  by  single  crystal  diffraction  techniques.  With  the  small  quantities  of  nanoparticle 
material  available  for  x-ray  analysis,  my  work  had  to  be  performed  on  a  micro-scale. 

The  nanoparticle  samples  were  packed  into  x-ray  diffraction  capillary  tubes  sold  by  the 
Charles  Supper  Company.  These  tubes  are  manufactured  from  silica  glass  in  an  assortment  of 
sized  ranging  from  0.3  to  0.7  mm  inside  diameter.  The  procedure  for  loading  a  sample  tube  was 
to  crush  a  small  (~5  mg)  portion  of  the  nanoparticle  sample  with  a  blunt  plastic  probe  inside  a 
clean  shell  vial.  This  treatment  was  necessary  to  break  up  nanoparticle  agglomerates  and 
facilitate  easier  packing  of  the  nanoparticles  into  the  glass  capillary  tube.  Following  this 
treatment,  the  nanoparticle  sample  was  scooped  into  the  broad,  open  end  of  the  capillary  tube 
and  packed  into  the  opposing  closed  end  of  the  tube  by  ultrasonic  vibration.  Once  a  column  of 
approximately  1.5  cm  of  nanoparticles  had  been  packed  into  the  closed  ent  of  the  tube,  the  excess 
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portion  of  the  sample  tube  was  removed  with  a  glass  cutting  tool.  The  open  end  of  the  tube  was 
sealed  with  a  small  portion  of  bee's  wax.  The  closed  end  of  the  packed  sample  tube  was  then 
fixed  into  a  brass  mounting  collet  sized  for  mounting  on  the  Gemini  diffractometer  goniometer 
sample  stage. 


Figure  33:  Gemini  x-ray  diffractometer  powder  sample  tube  preparation  steps. 


The  nanoparticle  diffraction  patterns  were  collected  for  a  range  of  20  angles  from  five  to 
sixty  degrees.  The  radiation  wavelength  was  set  to  the  copper  Ka  emission  line.  Diffractometer 
control  and  data  acquisition  was  managed  through  Oxford  Diffraction's  CrysAlis  Professional 
software.  The  diffraction  data  files  generated  were  written  in  an  ASCII  text  format,  and  the 
diffraction  data  they  contained  were  examined  with  data  plotting  software. 

2.5  Infrared  Observation  of  Electron  T ransitions 
2.5.1  Fourier  T ransform  IR  Spectroscopy 

M  id-infrared  spectroscopy  over  the  region  of  wavelengths  spanning  2.5  to  25  p,m  (2.5  ~ 
25  x  10"6  m)  is  most  commonly  used  for  the  identification  of  chemical  bonds.  For  a  bond  to 
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absorb  energy  in  this  region  of  the  electromagnetic  spectrum,  it  must  have  a  dipole  moment.  The 
absorbance  of  infrared  radiation  results  in  an  oscillatory  change  in  the  bond  geometry  as  its 
dipole  moment  is  changed.  These  bond  geometry  changes  are  observed  as  bond  stretching, 
bending,  or  scissoring  movements. 

Infrared  spectra  are  most  commonly  presented  as  energy  absorbance  over  a  frequency 
domain.  Because  of  the  large  frequencies  associated  with  micron  wavelength  radiation 
(approaching  hundreds  of  GHz),  infrared  spectra  have  historically  been  presented  with  frequency 
units  of  wavenumbers,  whose  units  are  cm"1. 

Infrared  spectroscopy  is  an  energy  dispersive  technique.  Prior  to  the  availability  of 
inexpensive  microcomputers,  infrared  energy  dispersion  for  spectroscopy  was  most  commonly 
accomplished  through  sequential  scanning  of  wavelengths  with  a  ruled  diffraction  grating  or 
prism  based  monochromator.  This  approach  generally  required  the  spectroscopist  to  make 
compromises  between  scanning  speed  and  quality  of  spectral  resolution.  Faster  scanning  of 
wavelengths  was  accomplished  by  loss  of  fine  spectral  detail.  The  time  requirements  of  this 
sequential  scanning  approach  also  made  it  impractical  for  applications  such  as  reaction 
intermediate  observation  or  kinetic  studies. 

The  replacement  of  a  grating  based  monochromator  in  the  spectrometer  with  an 
interferometer  accomplishes  a  significant  improvement  in  spectrum  acquisition  speed.  An 
example  of  one  common  design,  the  1891  M  ichelson  interferometer,  is  shown  in  Figure  34. 
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Figure  34:  The  M  ichelson  interferometer.  Note  that  the  compensating  plate  has  been  omitted  from  the 
stationary  mirror  path  Lx  for  clarity.  The  location  of  a  sample  for  spectrometry  would  be  in  the 
recombined  beam  prior  to  the  detector.184 


The  interferometer  operates  by  splitting  light  from  a  common  source  into  two  beams  of 
approximately  equal  intensity.  This  optical  division  is  accomplished  through  the  use  of  a 
partially  reflective  mirror,  known  as  a  beam  splitter.  Each  of  the  two  beams  is  directed  over  a 
separate  optical  path  and  then  recombined  immediately  prior  to  detection.  In  Figure  34  those 
separate  paths  are  designated  Li  and  L2.  By  alteration  of  one  optical  path's  length  (L 2)  through 
the  use  of  a  moving  mirror,  constructive  and  destructive  interferences  are  caused  within  the 
recombined  Li  and  L2  beams.  For  a  monochromatic,  i.e.  single  wavelength,  radiation  source  the 
output  of  the  interferometer  detector  is  modulated  by  the  varying  extent  and  type  of  Li  and  L2 
interference  in  the  recombined  beam,  as  the  mobile  mirror  is  linearly  translated.  The  resultant 
detector  output  is  an  oscillating  cosine  waveform  of  single  frequency. 

Spectroscopy  usually  requires  a  continuous  broadband  radiation  source  and  the  result  of 
scanning  the  interferometer  mirror  displacement  in  this  case  is  a  superposition  of  many  cosine 
waves,  each  waveform  having  a  different  frequency.  This  superposition  results  in  secondary 
interferences  and  the  creation  of  beat  frequency  overtones.  Consequently,  the  interferometer 
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detector  output  for  a  broadband  source  is  a  complex  signal  whose  intensity  varies  as  a  function  of 
mirror  spatial  displacement.  This  pattern  is  known  as  an  interferogram.  The  general  equation  for 
the  spatial  domain  interferogram  is: 

1(8 )  =  B(v)cos2nv8dv  (72) 

where  1(3)  is  the  signal  intensity  as  a  function  of  mirror  displacement  3,  and  B(v)  is  the  spectral 
absorbance  as  a  function  of  wavenumber  frequency  v.184  The  1(3)  term  is  the  interferogram  and 
th eB(v)  term  is  the  frequency  domain  spectrum  contained  within  it. 

In  order  to  extract  a  frequency  domain  spectrum  from  an  interferogram,  Fourier 
transform  techniques  are  applied  to  Equation  (72).  The  result  of  this  transformation  is  the 
equation: 


B(v)=2r  I(8)cos2Kv8d8  (73) 

JO 

This  equation  represents  the  frequency  domain  spectrum  that  was  acquired  through 
interferometry.  Obviously,  a  physical  instrument  cannot  make  infinite  mirror  translations  so  a 
compromise  is  necessary.  Scanning  of  the  mobile  mirror  in  a  spectrometer  is  typically  limited  to 
several  hundred  to  a  few  thousand  translations,  and  a  numerical  approximation  (i.e.  fast  Fourier 
transform  or  similar  techniques)  of  Equation  (73)  is  applied  through  software.  Execution  of 
these  actions  result  in  an  absorbance  infrared  spectrum  of  comparable  quality  to  a  sequential 
scan,  but  only  requiring  a  fraction  of  the  time  required  for  the  latter. 

Irradiation  of  a  sample  with  a  broadband  infrared  radiation  signal  and  application  of 
interferometry  techniques  gives  Fourier  transform  infrared  (FTIR)  spectroscopy  the  multiplex,  or 
Fellgett  advantage  over  monochromator  scanning.  The  entire  spectrum  is  acquired  in  a  parallel 
fashion  where  all  radiation  frequencies  and  their  associated  sample  absorbances  are  acquired 
simultaneously.  This  efficiency  gain  is  the  basis  for  the  FTIR  scan  speed  advantage  over 
monochromator  based  designs. 
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The  benefits  of  high  speed  in  spectroscopy  include  avoidance  of  changing  conditions  in 
the  instrument  or  sample  altering  the  spectrum  during  its  acquisition.  Rapid  acquisition  of 
spectra  also  open  the  opportunity  to  accumulate  multiple  spectra  which  can  be  averaged  or  added 
to  increase  the  signal-to-noise  ratio  of  the  data.  Fast  spectroscopy  also  makes  possible  the 
observation  of  transient  chemical  species  over  time  in  a  changing  reaction  environment,  making 
transition  state  intermediate  species  and  kinetic  studies  possible. 

2.5.2  The  Vacuum  FTIR  Cell 

The  high  surface  area  of  titanium  dioxide  nanoparticles,  ranging  from  10  to  100  m2  •  g  \ 
add  significant  difficulty  to  spectroscopic  studies.  Nanoparticles  rapidly  adsorb  materials  from 
the  air,  even  at  low  concentration.  This  capacity  for  gas  adsorption  is  one  of  the  strongest 
incentives  for  interest  in  their  application  to  heterogeneous  catalysis.  Hydrocarbon  adsorbates 
will  oxidize  on  the  surface  of  the  particles,  possibly  removing  oxygen  from  the  nanoparticle 
crystalline  lattice  and  reducing  the  titanium  cations  remaining.  The  electronegativity  of  oxygen 
in  water  adsorbates  will  depress  electron  population  in  the  conduction  band,  as  will  molecular 
oxygen  adsorbates.  Clearly,  any  attempt  at  spectroscopic  observation  of  electron  behavior  within 
titanium  dioxide  nanoparticles  requires  a  controlled,  adsorbate  free  nanoparticle  environment. 
The  vacuum  FTIR  cell  provides  that  control. 

The  vacuum  cell  apparatus  we  use  for  nanoparticle  FTIR  studies  was  designed  and  built 
at  the  University  of  Pittsburgh  for  the  research  group  of  Professor  John  Y  ates.  A  schematic  of 
the  cell,  in  cut  away  profile,  is  given  in  Figure  35.  The  FTIR  high  vacuum  cell  is  built  upon  a 
cubical  stainless  steel  box,  equipped  with  vacuum  flange  openings.  On  opposing  sides  of  the  cell 
are  potassium  bromide  IR  windows  and  on  the  front  face  of  the  cell  is  a  third  optical  port  (not 
shown  in  Figure  35)  equipped  with  a  sapphire  window.  This  third  port  is  used  as  a  light  injection 
window  for  conducting  ultraviolet  irradiation  experiments.  The  top  of  the  cell  is  extended 
approximately  25  cm  with  a  pair  of  stainless  steel  flange/pipe  assemblies.  Within  this  cylindrical 
extension  are  electrical  vacuum  pass-through  seals  for  supplying  power  and  thermocouple 
connections  to  the  nanoparticle  stage  inside  the  cell. 


81 


thermocouple  leads 


Figure  35:  FTIR  high  vacuum  cell.  Front  port  for  ultraviolet  irradiation  omitted  for  clarity. 


The  FTIR  vacuum  cell  can  accept  a  Dewar  sleeve  on  the  upper  cylindrical  extension  for 
liquid  nitrogen  cryogenic  cooling.  This  procedure  is  performed  for  trapping  of  trace  water  onto 
the  upper  walls  of  the  cell  during  nanoparticle  cleaning  and  conditioning  procedures,  in  advance 
of  making  experimental  observations  under  dehydroxyl ated  conditions. 

The  vacuum  cell  is  attached  with  a  flexible  metal  line  to  a  gas  manifold  that  is  used  for 
evacuating  the  cell,  monitoring  its  pressure,  and  allowing  the  controlled  introduction  of  gases  at 
low  partial  pressure  for  surface  adsorption  studies.  Vacuum  pumping  is  accomplished  with  a 
Pfeiffer  Vacuum  60  L's1  turbomolecular  pump  and  a  Varian  20  L*s_1  ion  pump  in  tandem. 
Normal  operating  pressure  for  the  FTIR  vacuum  cell  is  1  x  10"8Torr. 

The  FTIR  vacuum  cell  is  mounted  in  the  beam  path  of  a  M  attson  Instruments,  Inc.  M  odel 
RS-10000  Fourier  transform  infrared  spectrometer,  as  shown  in  Figure  36.  This  spectrometer  is 
equipped  with  a  cryogenically  cooled  mercury-cadmium-tellurium  (MCT)  infrared  detector. 
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FTIR  control  and  data  acquisition  is  managed  through  a  serial  interface  connection  to  a  desktop 
computer  running  M  attson  WinFirst®  spectroscopy  software.  The  vacuum  cell  is  attached  to  the 
spectrometer  sample  compartment  with  two  electromechanical,  stepping  motor  actuators  for 
vertical  and  horizontal  positioning.  The  spectrometer  optics  windows  within  its  sample 
compartment  are  adapted  with  polymer  rings  that  are  used  to  secure  a  nitrile  sleeve,  gas-tight 
connection  between  the  FTIR  and  the  vacuum  cell  windows  (indicated  as  a  and  b  in  Figure  36). 
Both  polymer  rings  are  equipped  with  side  tubing  connections  for  a  dry  nitrogen  gas  purge. 
These  flexible  nitrile  purged  tunnels  minimize  or  eliminate  atmospheric  moisture  contributions 
to  the  IR  spectrum  background  without  interfering  with  cell  translation. 


mesh  thermocouple  leads  1 1  high  current  power  leads 


Dewar  sleeve  (optional) 


to  vacuum  manifold 


z-axis  vacuum  cell  translation 


horozontal  axis  (in/out)  X 
vacuum  cell  translation  * 


Figure 36:  Vacuum  cell  installed  onto  Mattson  RS-10000  FTIR  spectrometer. 


The  fixture  for  mounting  of  nanoparticle  samples,  positioned  within  the  vacuum  cell, 
consists  of  two  nickel  clamps  that  are  attached  by  large  diameter  wiring  to  an  external  high- 
current  power  supply.  These  clamps  are  used  to  hold  a  piece  of  electrochemical ly  etched 
tungsten  mesh  in  the  optical  path  of  the  FTIR  and  provide  electrical  connection  of  opposing 
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edges  of  the  mesh  to  the  power  supply  circuit.  This  tungsten  mesh,  purchased  from  the  Bukbee- 
M  eyers  Co.,  is  formed  by  etching  a  repeating  pattern  of  square  holes  (0.3  mm  x  0.3  mm)  into  a 
thin  sheet  of  rolled  tungsten.  Under  magnification  the  resulting  etched  mesh  resembles  window 
screen  material  and  the  mesh  is  approximately  80%  transparent  to  light.  Details  of  the  clamp 
assembly  are  shown  in  Figure  37. 


to  vacuum  flange 
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power  leads  (2) 

nickel  mesh  clamps  (2) 

thermocouple  weld 

reference  position  (empty) 
sample  no.  1 

sample  no.  2 
tungsten  mesh 


Figure  37:  Tungsten  mesh  nanoparticle  stage  assembly. 


In  addition  to  its  optical  transparency,  tungsten  mesh  is  utilized  as  a  mounting  media  for 
nanoparticles  in  the  IR  cell  because  of  its  electrical  qualities.  Tungsten  is  more  resistive  than 
copper  (5.39  x  10"8  versus  1.71  x  10'8  Q  •  m  at  298  K  )3  which  causes  the  mesh  to  resistively  heat 
during  the  passage  of  electrical  current.  By  varying  the  current  applied  to  the  mesh,  the 
temperature  of  the  tungsten  mesh  and  the  nanoparticle  samples  pressed  into  its  orifices  can  be 
controlled  to  within  ±  IK. 

The  temperature  of  the  heated  tungsten  mesh  and  the  nanoparticles  pressed  into  it  is 
monitored  and  controlled  through  the  use  of  a  type  K  thermocouple  welded  onto  the  mesh.  The 
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thermocouple  junction  is  prepared  by  welding  an  approximately  six  inch  length  of  0.002  inch 
O.D.  alumel  wire  to  a  six  inch  length  of  0.002  inch  O.D.  chromel  wire  at  a  single  location.  The 
junction  is  then  spot  welded  directly  to  a  position  on  the  tungsten  mesh  above  the  SI  reference 
position  as  indicated  in  Figure  37.  The  thermocouple  signal  is  passed  out  of  the  cell  through 
vacuum  tight  electrical  feed-through  connections.  The  temperature  signal  is  read  on  a 
proportion-integral-derivative  (PID)  controller  which  regulates  the  supply  of  heating  current  to 
the  tungsten  mesh  and  maintains  the  mesh  at  a  setpoint  temperature.  The  temperature  setpoint  is 
ramped  by  either  a  manual  control  setting  or  use  of  a  programmable  voltage  ramp  generator. 

The  tungsten  mesh  fixture  within  vacuum  FTIR  cell  has  three  sample  positions.  The  top 
location,  indicated  in  Figure  37  by  a  broken  line  circle,  is  typically  left  empty.  An  infrared 
spectrum  of  this  bare  tungsten  mesh  location  is  recorded  at  the  beginning  of  an  experiment  and 
used  to  background  correct  the  nanoparticle  spectra.  When  a  new  mesh-nanoparticle  assembly  is 
loaded  into  the  vacuum  cell,  the  electromechanical  actuators  for  the  cell  mount  are  slewed 
through  a  range  of  horizontal  and  vertical  positions  for  each  nanoparticle  sample  location  on  the 
tungsten  mesh  while  the  peak-to-peak  voltage  of  the  FTIR  interferogram  is  monitored  for  a 
maximum  reading.  The  nanoparticle  sample  location  on  the  mesh  where  the  interferogram 
voltage  difference  maximizes  is  the  location  of  maximum  optical  density  for  that  sample.  These 
optical  density  maxima  for  both  samples  SI  and  S2  in  Figure  37  are  the  locations  where  infrared 
spectra  are  taken  during  experimental  observations,  as  these  are  the  sites  were  sample  path  length 
is  assumed  to  be  greatest.  Through  entry  of  the  numerical  coordinates  for  both  of  these 
absorbance  maxima  into  the  actuator  control  computer,  the  vacuum  cell  can  be  quickly 
repositioned  to  either  sample  location  during  an  experiment. 

2.5.3  Sample  Preparation  and  Conditioning 

A  new  piece  of  tungsten  mesh  is  cut  to  size  with  fine  scissors  using  either  a  paper 
template  or  a  previously  used  mesh  as  a  guide.  H  oles  are  located  and  punched  into  the  mesh  for 
mounting  on  the  vacuum  cell  nickel  clamps.  Once  the  mesh  is  cut  to  size,  it  is  cleaned  by 
ultrasonic  agitation  in  a  beaker  of  reagent  grade  methyl  alcohol  for  thirty  minutes,  followed  by 
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another  thirty  minutes  in  a  beaker  of  reagent  grade  acetone.  Following  the  second  solvent 
cleaning,  the  mesh  is  allowed  to  air  dry  while  loosely  wrapped  in  laboratory  tissue  wipes. 

Prior  to  loading  nanoparticle  samples  onto  the  mesh,  a  temperature  monitoring 
thermocouple  must  be  fabricated  and  attached  directly  to  the  mesh.  Two  pieces,  each 
approximately  six  inches  in  length,  of  0.002  inch  O.D.  chromel  and  alumel  wire  arejoined  with  a 
resistance  spot  weld,  leaving  a  one-half  inch  length  of  each  wire  extending  beyond  the  weld 
joint.  These  short  extensions  of  thermocouple  wires  are  threaded  into  two  adjacent  holes  in  the 
tungsten  mesh  at  a  location  just  below  the  point  of  attachment  for  the  thermocouple  to  the  mesh 
surface.  The  electrodes  of  the  resistance  welding  tool  are  then  used  to  press  the  thermocouple 
junction  against  the  tungsten  mesh  surface  and  the  junction  is  spot  welded  onto  the  mesh. 
Following  attachment  of  the  thermocouple,  excess  wire  threaded  through  the  mesh  is  trimmed 
and  the  response  of  the  thermocouple  is  checked  with  a  digital  multimeter.  If  the  thermocouple 
shows  a  voltage  response  with  temperature  change,  the  mesh  is  ready  for  nanoparticle  sample 
loading. 

The  tungsten  mesh  and  attached  thermocouple  are  placed  within  a  custom  press  mold 
assembly,  shown  in  Figure  38.  With  the  mesh  correctly  positioned,  the  top  half  of  the  mold  is 
secured  with  machine  screws.  A  5  mg  portion  of  each  nanoparticle  sample  is  loaded  into  its 
position  (S2  or  S3  in  Figure  38)  and  the  piston  is  seated  in  the  mold  above  the  loose 
nanoparticles.  The  mold  and  piston  are  then  placed  on  a  hydraulic  press  and  the  piston  is 
compressed  at  13  kpsi  (1.3  x  104  lb*in"2)  for  two  minutes.  Following  nanoparticle  loading  of  the 
tungsten  mesh  sample  positions,  the  top  half  of  the  mold  is  removed  and  the  mesh  is  mounted  on 
the  vacuum  cell  nickel  clamps. 
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Figure  38:  Nanoparticle  press  mold  assembly. 


Following  placement  of  the  tungsten  mesh  and  nanoparticles  into  the  vacuum  cell,  the 
cell  is  temporarily  wrapped  with  laboratory  resistive  heating  tapes  and  attached  to  the 
aforementioned  FTIR  vacuum  system.  The  cell  is  heated  to  approximately  100  °C  and 
evacuated  with  the  turbomolecular  pump  for  approximately  twelve  to  twenty-four  hours  for 
water  removal. 

After  the  initial  evacuation  of  the  vacuum  cell,  the  nanoparticles  are  subjected  to  oxygen 
treatments  for  the  oxidative  removal  of  adsorbed  organic  materials.  The  heat  tape  is  removed  and 
the  vacuum  cell  is  reinstalled  into  the  FTIR  sample  compartment.  Vacuum  pumping  of  the  cell 
is  continued  and  the  mesh  is  resistively  heated  to  675  K  over  four  hours.  The  vacuum  pumping 
is  discontinued  and  the  cell  is  pressurized  to  20  Torr  with  oxygen  while  resistive  mesh  heating  is 
continued  at  675  K  for  an  additional  hour.  Following  the  first  oxygen  treatment,  the  vacuum  cell 
is  pumped  to  1  x  10'6  Torr.  The  cell  is  once  again  pressurized  to  20  Torr  with  oxygen  and 
allowed  to  continue  heating  at  675  K  for  an  additional  30  minutes.  Following  this  second 
oxygen  treatment  the  cell  is  pumped  back  to  1  x  10'6  Torr  and  allowed  to  cool  to  room 
temperature.  The  effect  of  these  treatments  is  to  remove  water  and  hydrocarbon  adsorbates  from 
the  nanoparticles. 
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2.5.4 Temperature  Ramp  Infrared  Spectra  Collection 

Thermal  ramping  of  the  titanium  dioxide  nanoparticles  within  the  FUR  vacuum  chamber 
is  ordinarily  performed  to  either  observe  the  condition  of  the  nanoparticle  samples  following 
adsorption  studies  or  to  remove  water  residues  from  the  particles.  The  vacuum  cell  is  not 
continually  pumped  when  it  is  not  in  use,  but  rather  an  isolation  valve  on  the  line  from  the  cell  to 
the  vacuum  manifold  is  closed.  This  mode  of  idling  the  apparatus  avoids  the  necessity  of 
continual  operation  of  the  turbomolecular  and  ion  pumps.  However,  even  when  the  cell  is  under 
active  vacuum  pumping,  a  slow  ingress  of  water  into  the  system  from  the  outside  environment  is 
unavoidable.  The  large  surface  area  of  the  nanoparticle  samples  within  the  vacuum  cell  capture 
and  accumulate  water  from  this  continual  flux  of  low  level  water  leakage.  A  technique  has  been 
developed  to  dehydroxylate  the  nanoparticle  surfaces  for  sufficient  intervals  that  permit 
spectroscopic  experiments  over  several  hours  under  nearly  water-free  conditions. 

A  Dewar  vessel  was  designed  (see  Figure  35  and  Figure  36)  to  fit  over  the  cylindrical 
extension  of  the  vacuum  cell  body.  This  insulated  vessel  seals  against  the  stainless  steel  body  of 
the  cell.  During  its  use  the  Dewar  sleeve  is  filled  with  a  coolant,  usually  liquid  nitrogen.  Once 
cooled  to  ~  77  K  by  the  liquid  nitrogen,  the  upper  regions  of  the  vacuum  cell  become  a  water 
trapping  surface.  The  nanoparticle  samples  temperature  is  then  resistively  heated  to 
approximately  400  K  while  the  FTIR  cell  is  evacuated  and  the  cryogenic  water  trapping  Dewar 
sleeve  is  in  place.  The  trace  level  of  adsorbed  water  on  the  nanoparticle  surfaces  is  thermally 
desorbed  by  resistive  heating  and  upon  contact  with  the  cryogenically  cooled  vacuum  cell  wall, 
the  vapor  condenses  and  is  trapped  as  ice.  Following  this  treatment,  the  nanoparticles  exist  in  a 
highly  dehydroxyl ated  state. 

Temperature  ramping  experiments  of  the  FTIR  vacuum  chamber  are  generally  conducted 
under  automated  control.  Initially  a  spectrum  is  taken  of  the  'bare'  tungsten  mesh  SI  region  (see 
Figure  37  and  Figure  38)  at  ambient  temperature  (298  K).  This  spectrum  is  used  for  background 
subtraction  from  the  subsequent  spectra.  This  background  subtraction  is  intended  to  remove 
spectral  contributions  of  the  tungsten  and  tungsten  oxide  from  the  titanium  dioxide  nanoparticle 
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spectra.  The  high-current  power  supply  PID  controller  setpoint  is  then  switched  over  to 
automatic  ramping  from  an  external  voltage  ramp  circuit.  As  the  setpoint  voltage  is  increased  by 
the  external  voltage  ramp  generator,  the  electric  current  passing  through  the  tungsten  mesh  is 
increased  by  the  PID  controller  and  resistive  heating  of  nanoparticle  samples  increases. 
Temperature  read-back  from  the  tungsten  mesh  thermocouple  is  displayed  on  the  PID  controller 
unit  and  is  recorded  manually,  concurrent  with  spectra  acquisition  by  the  FTIR  data  system. 
Spectral  data  files  are  recorded  in  comma  separated  value  (CSV)  ASCII  text  files. 
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Chapter  3:  Data  Analysis 


3.1  E  lectron  M  icroscopy  Particle  Diameter  M  easurements 

The  chemical  reactivity  and  physical  properties  of  nanoparticles  are  size  dependent. 
T hese  effects  make  accurate  measurement  of  nanoparticle  size  and  morphology  an  important  part 
of  characterization.  I  chose  transmission  electron  microscopy  (TEM )  for  these  measurements 
because  it  offered  a  relatively  simple,  inexpensive,  and  direct  method  of  obtaining  this  data. 

All  nanoparticle  diameter  data  was  obtained  from  the  Zeiss  IOC  A  transmission  electron 
microscope  at  the  Virginia-M  aryland  Regional  Veterinary  School.  This  instrument  produces  a 
computer  graphics  file  that  the  user  copies  and  examines  offline.  The  file  format  is  the  tagged 
image  file  format  (TIFF)  which  can  be  opened  and  manipulated  by  most  graphics  editing 
programs.  Each  file  that  this  microscope  generates  includes  an  image  scale  bar  which  can  be 
used  to  calibrate  image  processing  software. 

Nanoparticle  diameter  measurements  are  made  from  the  TIFF  files  with  the  ImageJ 
software  package207,  distributed  by  the  National  Institute  of  Health.  Prior  to  making  particle 
diameter  measurements,  the  software  is  calibrated  by  drawing  a  measurement  line  along  the 
length  of  the  TEM  image  calibration  bar.  The  length  of  the  calibration  bar  and  its  measurement 
units  are  then  entered  into  the  program.  Following  this  image  scale  calibration  of  the  software, 
all  measurements  made  for  that  TIFF  file  produce  calibrated  results. 

Particle  diameter  measurements  are  accomplished  in  ImageJ  by  drawing  a  line  with  a 
computer  mouse  cursor  across  the  particle  followed  by  execution  of  the  'Measure'  command. 
This  command  causes  the  software  to  calculate  the  particle  diameter  in  the  units  previously 
calibrated  from  the  TIFF  image  calibration  bar,  and  append  the  result  to  an  ANSII  text  file  of 
measurement  results.  Because  the  TEM  images  commonly  contain  hundreds  to  thousands  of 
nanoparticle  features,  I  wrote  a  macro  in  ImageJ  that  executed  the  'measure'  command  from  a 
keystroke,  and  then  uniquely  labeled  the  measured  particle  diameter  to  avoid  double 
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measurements  of  particles  within  an  image.  At  the  conclusion  of  diameter  measurements  for  a 
TIFF  file,  the  file  was  saved  with  an  appended  file  name  for  archival  purposes  and  the  diameter 
measurements  text  file  was  saved  for  subsequent  processing. 

WhileTEM  was  likely  the  most  accurate  method  of  obtaining  nanoparticle  diameters,  the 
images  were  usually  not  completely  accessible  for  diameter  measurements.  Particle 
agglomeration,  particularly  for  smaller  diameter  nanoparticles,  lead  to  opaque  aggregates  of 
material  with  largely  obscured  particle  boundaries.  In  these  cases,  diameter  measurements  were 
made  for  particles  on  the  periphery  of  the  aggregates,  where  particle  boundaries  were  more 
clearly  defined.  In  spite  of  this  limitation  of  some  images,  it  was  still  common  to  obtain  50-200 
measurements  from  a  single  TIFF  file.  A  summary  table  of  particle  diameter  measurements  and 
particle  synthesis  conditions  is  included  in  Appendix  5,  T able  6  on  page  127. 

Control  of  nanoparticle  size  and  size  distribution  is  a  concern  common  to  most 
researchers  working  with  these  materials.  The  laser  vaporization  technique  that  I  used  offered 
two  possible  approaches  to  particle  diameter  control. 

The  first  method  of  diameter  control  attempted  was  control  of  the  linear  distance  from  the 
precursor  material  pellet  undergoing  vaporization  to  the  nanoparticle  collector  surface. 
Following  some  initial  experimentation  with  sheet  stainless  steel  collector  surfaces  that  were 
arranged  in  a  semi-cylindrical  geometry,  a  hemispherical  geometry  was  selected  for  collector 
surfaces.  A  hemisphere  provides  a  symmetrical  geometry  and  an  equal  path  length  from  the 
precursor  pellet  located  at  its  center  to  the  collector  surface. 

Two  experiments  were  performed  with  stainless  steel  collector  hemispheres  of  different 
radii  (2.3  and  6.3  cm).  Nanoparticles  were  prepared  with  both  hemispheres  under  identical 
vaporization  conditions  (37  Torr  02, 100%  laser  power,  1  hour  vaporization  time).  The  results  of 
these  experiments  are  summarized  in  Appendix  5,  Table  6  for  the  SPB014  and  SPB015 
nanoparticle  batches.  Although  the  larger  collector  radius  (6.3  cm)  produced  a  1  nm  larger  mean 
particle  diameter,  the  ±3  -  4  nm  error  in  the  diameter  measurements  for  these  batches  makes  any 
conclusion  that  path  distance  significantly  influenced  their  size  questionable. 
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The  second  approach  to  influencing  nanoparticle  size  was  vaporization  chamber  pressure 
control.  For  these  experiments,  also  summarized  in  Table  6,  a  hemispherical  stainless  steel 
collector  surface  of  4.4  cm  radius  was  used.  Following  three  cycles  of  evacuation  and  ultra  high- 
purity  (U H P)  oxygen  flushing,  the  vaporization  chamber  was  backfilled  with  UFIP  oxygen  and 
held  at  pressures  varying  from  1  to  500  Torr  for  each  experiment.  In  all  experiments,  the  laser 
was  operated  at  100%  power  output  (60  Watts)  and  vaporization  was  performed  for  one  hour. 
The  data  from  these  experiments  are  plotted  in  Figure  39. 

While  no  strongly  correlated  statistical  trend  is  demonstrated  from  the  plot  in  Figure  39, 
the  data  does  appear  to  support  the  generalization  that  higher  chamber  pressures  lead  to  larger 
nanoparticle  diameters.  Another  effect  that  appears  to  increase  with  chamber  pressure  is  the 
distribution  of  particle  diameters  produced.  H  igher  pressure  (500  Torr)  growth  appears  to  nearly 
approximate  a  bimodal  particle  size  distribution,  as  is  demonstrated  by  the  error  bars  and  TEM 
micrograph  insets  of  Figure  39.  Nanoparticles  produced  at  low  pressure  (5  Torr)  appear,  in 
comparison,  to  approach  monodisperse  characteristics. 

These  pressure-size  effects  can  be  explained  by  a  decrease  in  mean  free  path  for  particles 
of  material  vaporized  from  the  precursor  pellet,  as  chamber  pressure  is  increased.  As  the 
chamber  pressure  is  increased  the  frequency  of  particle-particle  and  particle-oxygen  molecule 
collisions  are  also  increased.  This  higher  collision  frequency  causes  an  increasingly  Boltzmann- 
like  distribution  of  nanoparticle  velocities.  The  slower  particles  experience  more  collisions  and 
grow  to  larger  sizes  before  they  impact  the  collector  surface.  Also,  a  stationary  large  diameter 
particle  on  the  collector  surface  continues  to  offer  a  larger  target  for  collisions  with  gaseous 
material,  so  some  extent  of  size-skewed  particle  growth  might  continue. 
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Figure  39 :  Nanoparticle  mean  diameter  asa  function  of  laser  vaporization  chamber  02  pressure.  Data 
given  in  Appendix  5,  Table  6.  The  scale  bars  on  the  inset  micrographs  all  indicatea  100  nm  dimension. 


My  calculations  of  titanium  dioxide  nanoparticle  shallow  trap  state  energies  were 
obtained  from  two  nanoparticle  synthesis  batches,  SPB018  (4.0  ±  0.6  nm  mean  diameter)  and 
SPB020  (30  ±  10  nm  mean  diameter).  The  diameter  measurements  data  for  these  batches  are 
included  in  Appendix  5,  Table  6.  TEM  micrographs  and  particle  size  histograms  are  presented 
in  Appendices  7  and  8  on  pages  128  and  129.  A  high  resolution  TEM  image  from  the  Titan  300 
microscope  at  the  Virginia  Tech  Institute  for  Critical  Technology  and  Applied  Science  (ICTAS) 
is  presented  in  A  ppendix  8  on  page  130. 
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3,2  Powder  X-ray  Diffraction 

Powder  x-ray  diffraction  (XRD)  patterns  of  nanoparticles  are  used  for  identification  of 
phase  composition  within  the  particles,  i.e.  confirmation  of  the  presence  of  the  rutile  and/or 
anatase  polymorphs  and  their  relative  amounts.  Nanoparticle  mean  diameter  can  be  indirectly 
obtained  from  XRD  through  the  use  of  the  Scherrer  equation: 


PW,  = 

2 


0.941 

LcosO 


(74) 


where  PW1/2\s  the  reflection  peak  width  at  half  height,  X  is  the  x-ray  radiation  wavelength,  and  0 
is  the  Bragg  angle  of  incidence.208,209  The  particle  diameter  is  the  variable 

Another  motive  for  collection  of  XRD  patterns  for  nanoparticles  is  confirmation  of 
nanoparticle  crystallinity.  Crystallinity  is  less  of  a  concern  for  laser  vaporization  particle 
synthesis  as  it  is  for  sol-gel  precipitation  techniques,  as  the  former  technique  provides  its  own 
heat  treatment  to  crystallize  condensed  material. 

M  y  powder  XRD  analyses  were  performed  on  particles  held  within  glass  capillary  tubes. 
Although  glass  is  an  amorphous  material,  it  still  contributes  a  broad  x-ray  diffraction  peak  and  all 
of  my  XRD  patterns  are  superimposed  upon  this  background  contribution.  Figure  40  is  an  XRD 
pattern  from  an  empty  0.3  mm  I.D.  sample  capillary. 

The  use  of  pattern  subtraction  was  attempted  for  removal  of  the  sample  tube  background 
contribution  to  the  XRD  pattern.  Differences  in  tube  orientation  with  respect  to  the  x-ray  beam 
and  variations  in  tube  diameter  made  subtraction  of  a  'blank'  from  an  empty  capillary  from  a 
diffraction  pattern  lead  to  negative  peak  deflections.  The  chance  of  inducing  artificial  features  to 
the  XRD  pattern  or  removing  useful  information  led  to  my  decision  to  refrain  from  using  this 
type  of  background  subtraction. 
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Figure  40:  X-ray  diffraction  pattern  for  an  empty,  0.3  mm  ID  sample  capillary. 


Another  intriguing  approach  to  removal  of  the  silica  background  was  a  data  processing 
treatment  based  on  the  work  of  Sonneveld  and  Visser.210  I  wrote  a  computer  program  that 
applied  their  baseline  isolation  algorithm  for  subtraction  of  silica  contributions  to  the  XRD 
patterns.  The  algorithm  (and  my  program)  required  a  user  input  for  the  number  of  baseline 
smoothing  cycles  applied  to  the  extracted  background.  M  y  work  with  this  program  showed  that 
the  excessive  use  of  this  baseline  smoothing  could  potentially  remove  some  useful  XRD  pattern 
information,  but  an  example  of  the  results  of  Sonneveld-Visser  baseline  correction  is  included  in 
Figure  41. 
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Figure  41:  Degussa  P25  titanium  dioxide  nanoparticle  XRD  pattern  with  Sonneveld  -  Visser  baseline 
correction.  Pattern  (a)  is  the  uncorrected  data,  (c)  is  the  extracted  and  smoothed  baseline,  and  (c)  is 
the  baseline  corrected  XRD  pattern. 


XRD  characterization  of  nanoparticles  was  limited  to  the  Degussa  P25  titanium  dioxide, 
SPB018  (4.0  ±  0.6  nm  mean  diameter)  and  the  SPB041  nanoparticle  synthesis  batches.  In  the 
case  of  SPB041  electron  microscopy  characterization  has  not  yet  been  performed,  but  its 
synthesis  conditions  (100  Torr  O2,  4.4  cm  collector  radius)  are  consistent  with  nanoparticles  in 
the  10  ~  20  nm  diameter  range.  The  Figure  42b  XRD  pattern  for  SPB041  is  consistent  with  a 
mean  particle  diameter  at  the  upper  end  of  the  10  ~  20  nm  range,  as  is  apparent  when  comparing 
the  pattern  to  Figure  42a  for  Degussa  P25  which  has  a  mean  particle  diameter  of  25  nm. 
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Figure  42:  X-ray  diffraction  nanoparticle  characterization.  Patterns  are  for  (a)  Deguss  P25,  (b) 
SPB041,  and  (c)  SPB018  titanium  dioxide.  Reflections  identified  with  '?'  are  tentative. 


The  XRD  pattern  reflection  assignments  are  based  on  the  work  of  Zhou  et  al.211,  where 
Degussa  P25  was  characterized.  The  patterns  in  Figure  42  were  offset  from  each  other  for  clarity 
in  the  pattern  plot,  but  no  other  modification  of  the  data  has  been  performed.  It  is  interesting  to 
note  that  for  Figure  42c,  from  the  4.0  ±  0.6  nm  diameter  particles,  the  dominant  titanium  dioxide 
polymorph  appears  to  be  rutile,  with  very  little  or  no  anatase  present  in  the  pattern.  This 
outcome  would  seem  to  be  in  conflict  with  the  earlier  observation  that  anatase  is  the 
thermodynamically  favored  polymorph  for  particle  diameters  below  14  nm. 81-83  This  phenomena 
might  be  a  characteristic  of  particles  made  by  laser  vaporization  synthesis,  and  it  merits  further 
investigation. 
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3.3  FTI R  Data  Processing  and  Shallow T rap  Energy  Calculation 

M  y  calculation  of  titanium  dioxide  nanoparticle  shallow  trap  state  energy  is  based  upon  a 
mathematical  treatment  of  infrared  spectra.  Conduction  band  electron  population  effects  on  the 
infrared  spectrum  baseline  are  observed  over  all  infrared  frequencies  between  4000  to  1000  cm"1. 
Although  the  temperature  ramp  spectra  were  recorded  under  vacuum,  spectral  absorbances  from 
residual  water  on  the  nanoparticles  in  the  vacuum  cell  is  still  observed.  For  this  reason  I  decided 
to  use  a  region  of  the  infrared  spectrum  from  2600  to  1900  cm"1  for  my  calculations.  This  region 
is  isolated  from  hydrogen  bonded  OH  absorbance  over  3700  to  2600  cm"1  and  water  H-O-H 
bending  modes  at  1650  cm"1. 

The  Mattson  WinFirst  spectroscopy  software  saves  infrared  spectra  as  space-delimited 
ASCII  text  files.  Each  file  contains  several  hundred  data  points,  and  collectively  there  were 
approximately  one  hundred  or  more  files  to  examine.  My  approach  was  to  write  a  computer 
program  in  G77  FORTRAN  (WaveShave.f,  Appendix  3)  which  used  a  text  file  list  of  data  files 
to  direct  sequential  processing  of  spectra.  Each  data  file  on  the  list  was  opened  and  read  by 
WaveShave.f  and  rewritten  as  a  second  text  file  which  only  contained  spectral  data  points 
between  my  predefined  range  of  infrared  frequencies,  2600  to  1900  cm"1.  When  the  file 
processing  run  was  complete,  a  new  set  of  truncated  infrared  spectrum  data  files  had  been 
created. 

The  association  of  infrared  spectrum  baseline  area  given  in  Equations  (63)  to  (66)  with 
the  Fermi-Dirac  behavior  of  conduction  band  electrons  requires  an  integrated  area  under  the 
spectrum  baseline  for  each  truncated  spectrum  file  produced  by  data  processing  with 
WaveShave.f.  I  wrote  a  second  G77  FORTRAN  program  (T rapl nt2.f.  Appendix  4)  which  was 
adopted  from  a  program  given  in  Nyhoff  and  Leestma's  FORTRAN  77  for  Engineers  and 
Scientists. 212  My  program  applied  trapezoidal  numerical  integration  to  the  truncated  spectra 
baselines  previously  generated  by  WaveShave.f.  Trapezoidal  integration  is  based  upon  the 
approximation: 
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(75) 


Ja  f(x)dx  «  ^[/(x0)  +  2f(xl)  +  2/U2)  +  ••  •  +  2f(xn  l)  +  (xn)] 

where  zlx  =  (b  -  a)  / n  and  Xj  =  a  +  i/lx. 213  The  variables  a  and  b  were  the  starting  and  ending 
frequencies  of  the  truncated  spectrum,  n  was  the  number  of  data  points  contained  within  the  data 
file,  and  Jxwas  the  spacing  between  absorbance  data  points  within  the  file.  Like  WaveShave.f, 
the  program  T rapl nt2.f  uses  an  ASCII  text  file  containing  a  list  of  data  files  (those  previously 
generated  by  WaveShave.f)  to  direct  processing.  The  integration  results  from  the  batch 
processing  are  written  to  a  single  ASCII  text  file  as  a  summary  of  data  file  names  and  their 
associated  integration  result. 

Following  integration  of  the  truncated  FTIR  spectra,  a  curve  fit  to  the  Fermi-Dirac 
relation  previously  given  in  Equation  (64)  was  necessary  (see  page  61).  The  baseline  integration 
areas  were  matched  to  their  associated  sample  temperatures  and  these  data  were  fit  to  the  Fermi- 
Dirac  relation  using  a  non-linear  least  squares  technique  described  by  Brown.214, 215  This  curve 
fit  technique  is  performed  with  a  spreadsheet  constructed  in  Microsoft  Excel®.  Reiterative 
processing  of  the  data  set  within  the  spreadsheet  is  accomplished  with  the  Solver®  application 
that  is  included  with  the  spreadsheet  program.  The  variables  of  the  regression  are  the  integrated 
baseline  area  Z  and  the  nanoparticle  temperature  T.  The  regression  parameters  are  the 
proportionality  constant  D  and  the  shallow  trap  energy  state  Estate  An  example  of  my 
spreadsheet  constructions  is  shown  in  Figure  43. 
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Figure  43:  Nonlinear  regression  curve  fit  spreadsheet  for  theFermi-Dirac  relationship. 

Brown’s  approach  is  to  use  Solver®  to  reiteratively  adjust  the  fit  parameters  specified  by 
the  user  while  simultaneously  seeking  to  maximize  coefficient  of  determination,  RSQ  in  Figure 

43.  For  perfect  correlation  of  the  independent  and  dependent  variables  (T  and  Z  in  my  work) 
RSQ  would  be  one.  A  summary  of  the  regression  results  for  the  SPB018  (4.0  ±  0.6  nm) 
nanoparticles  is  given  in  Table  3  and  for  the  SPB020  (30  ±  10  nm)  nanoparticles  in  Table  4. 
Plots  of  the  regression  data  sets  and  their  fitted  Fermi  -  Dirac  trend  lines  are  provided  in  Figure 

44. 
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Regression  Parameters 

Data  Set 

D 

( A.U. •  cm2*  ne) 

Estate 

(eV) 

RSQ 

11730/07 

3.48  xlO18 

0.290 

0.999 

12/03/07 

3.10X1018 

0.277 

0.996 

12/04/07 

1.22  x  1017 

0.457 

0.993 

12/10/07 

5.01  x  1018 

0.302 

0.992 

12/20/07 

7.14  x1a18 

0.377 

0.997 

Mean  (eV): 

0.34 

Standard  Deviation  (eV):  0.08 


Table  3:  Fermi  -  Dirac  regression  results  for  SPB018  (4.0  ±0.6  nm)Ti02  nanoparticles. 


Regression  Parameters 

Data  Set 

D 

( A.U. •  cm2*  ne) 

Estate 

(eV) 

RSQ 

12/05/07 

9.84  x  10" 18 

0.398 

0.992 

12/06/07 

2.98  x  1018 

0.259 

0.976 

12/18/07 

2.65  x  1018 

0.207 

0.998 

Mean  (eV): 

0.3 

Standard  Deviation  (eV): 

0.1 

Table 4:  Fermi  -  Dirac  regression  results  for  SPB020  (30  ±10  nm)Ti02  nanoparticles. 
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Temperature  /  K 


Figure  44:  Integrated  area  vs.  temperature  plots  for  4  nm  (circle  symbols)  and  30  nm  (triangle 
symbols)  Ti02.  Trend  lines  are  fitted  Fermi  -  Dirac  relation.  Data  sets  designated  with  (a)  were 
collected  with  cryogenic  water  trapping. 
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Chapter  4:  Discussion  of  Results  and  Conclusions 


41  Review  of  Research  Objectives 

Titanium  dioxide  nanoparticles  catalyze  oxidation-reduction  chemistry  by  using  external 
energy  inputs  (heat  or  electromagnetic  radiation)  to  promote  electron  transport  between  the 
particle  surface  and  adsorbates.  This  catalytic  cycle  is  initiated  by  excitation  of  electrons  to  the 
titania  conduction  band,  causing  the  migration  of  both  electrons  and  electron  vacancies  (i.e. 
'holes')  to  the  particle  surface.  Thermal  excitation  of  electrons  provides  an  alternative  to  high- 
intensity  ultraviolet  irradiation  of  nanoparticles  for  promotion  of  adsorbate  reaction  catalysis. 
The  electrons  entering  the  titania  conduction  band  as  the  result  of  thermal  excitation  originate 
from  structural  defects  in  the  titanium  dioxide  crystal  lattice  structure  that  are  collectively  named 
shallow  trap  states.  Characterization  of  the  energy  requirements  for  thermal  excitation  of 
electrons  from  these  shallow  trap  states  is  important  to  understanding  and  optimizing  titanium 
dioxide  nanoparticle  catalysis.  M  y  work  sought  to  determine  the  shallow  trap  state  energies  for 
titanium  dioxide  nanoparticles  that  were  formed  through  laser  vaporization-condensation 
synthesis  techniques.  M  y  measurement  of  shallow  trap  state  energies  was  accomplished  through 
collection  of  transmission  infrared  spectra  of  the  nanoparticles  over  temperatures  spanning  298 
to  675  K,  followed  by  evaluation  of  the  spectroscopic  data  with  an  expression  derived  from 
Fermi-Dirac  electron  theory. 

42  Nanopartide  Sizes  and  Size  Distributions 

One  early  idea  for  directing  the  size  of  titanium  dioxide  nanoparticles  produced  from 
laser  vaporization  was  to  vary  the  radial  distance  from  the  precursor  material  pellet  to  the 
collector  surface.  By  increasing  the  volume  in  which  accreting  particles  could  undergo  collisions 
with  gaseous  titanium  dioxide  species,  it  was  hoped  that  larger  particles  could  be  generated. 

Two  hemispherical  steel  collector  surfaces  were  prepared  with  radii  of  2.3  and  6.3  cm. 
The  4.0  cm  difference  in  radii  correspond  to  a  volume  difference  of  approximately  500  cm3. 
Nanoparticle  batches  were  prepared  with  each  collector  at  a  chamber  pressure  of  37  Torr.  The 
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particle  diameter  data  for  these  batches,  SPB014  and  SPB015  in  Table  6  on  page  127,  suggests 
that  reducing  the  collector  volume  by  approximately  one-half  liter  had  negligible  effect  on  the 
nanoparticle  mean  diameter.  The  difference  in  diameters  between  these  batches  was 
approximately  8%.  Electron  micrographs  for  both  batches  are  presented  in  Figure  45. 


a  b 

Figure  45:  T ransmission  electron  micrographs  for  SPB014  and  SPB015  titanium  dioxide  nanoparticle 
batches.  The  precursor  pellet  to  collector  surface  distance  was  (a)  6.3  cm  and  (b)  2.3  cm.  Both 
syntheses  conducted  at  37  Torr  0  2  chamber  pressure. 


The  effect  of  pressure  manipulation  on  nanoparticle  mean  diameter  is  much  easier  to 
discern.  While  the  data  in  T able  6  for  nanoparticle  syntheses  with  the  4.4  cm  radius  collector  at 
different  pressures  is  not  well  correlated  (R2  =  0.86  in  Figure  39),  there  is  a  monotonic  increase 
in  particle  diameter  with  increasing  vaporization  chamber  pressure.  There  also  is  a  trend  of 
increasing  particle  size  distribution  (note  the  increasing  standard  deviations  with  pressure  in 
T  able  6). 

As  the  chamber  pressure  increases,  the  mean  free  path  of  the  gaseous  species  within  the 
chamber  reduces.  Atkins  and  de  Paula  note  that  mean  free  path  reduces  by  one-half  for  each 
doubling  of  pressure.90  For  the  Table  6  data  vaporization  experiments  were  performed  at 
pressures  between  1  and  500  Torr  O2.  The  difference  in  these  pressures  represents 
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approximately  nine  doublings  of  pressure.  The  corresponding  nine  reductions  of  the  mean  free 
path  by  factors  of  one-half  would  have  a  cumulative  effect  of  a  2  x  10"3  reduction.  This 
magnitude  of  mean  free  path  reduction  translates  into  greatly  increased  collision  frequency 
between  accreting  nanoparticles  and  gaseous  material  from  the  liquefied  precursor  pellet. 

Increased  collisions  in  the  gas  phase  from  increased  chamber  pressure  also  provide  a 
greater  opportunity  for  M  axwell-Boltzmann  statistics  to  broaden  the  distribution  of  particle 
diameters.  The  gaseous  species  within  the  vaporization  chamber  partitions  in  a  Boltzmann 
distribution  at  any  of  the  pressures  used  for  nanoparticle  synthesis.  At  higher  pressures,  with 
increased  collision  frequency,  that  distribution  of  energies  is  imparted  upon  the  accreting 
particles.  The  result  is  that  the  particles  acquire  a  Boltzmann  distribution  of  energies  which 
affects  their  number  of  collisions.  For  example,  a  particle  impacted  by  more  high  kinetic  energy 
particles  will  itself  acquire  a  higher  kinetic  energy  than  an  identical  particle  that  experienced 
fewer  high  energy  impacts.  The  distribution  of  nanoparticle  speeds  leads  to  a  distribution  of 
growth  rates  that  broaden  the  nanoparticle  size  distribution.  This  effect  appears  to  be  directly 
proportional  to  the  vaporization  chamber  pressure. 

4.3  Shallow  T rap  State  E  nergy  M  easurements 

The  Fermi-D irac  regression  results  for  4  nm  and  30  nm  particles  were  considered 
separately  in  Table  3  and  Table  4  respectively,  on  page  101.  Chestnoy  et  al.  216  observed  that 
when  the  diameter  of  3  nm  cadmium  sulfide  nanoparticles  approached  the  diameter  of  their 
shallow  trap  states,  strong  mixing  of  the  trap  state  and  the  conduction  band  wavefunctions  occur. 
This  mixing  suggests  that  there  may  be  a  size  effect  on  the  shallow  trap  state  energies  observed 
for  particles  of  greatly  different  size.  M  y  observations  do  not  indicate  a  size  effect  on  shallow 
trap  state  energies.  The  mean  Estate  values  in  T able  3  and  T able  4  are  essentially  the  same  (0.34 
±  0.08  and  0.3  ±  0.1  eV  respectively).  Flowever,  it  may  be  inappropriate  to  consider  this 
comparison  as  a  test  of  size  effects  on  shallow  trap  state  energies,  as  the  larger  particles 
represented  in  T able  4  had  a  much  larger  distribution  of  diameters  than  the  particles  in  T able  3. 
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Shallow  trap  state  energies  reported  by  different  researchers  for  titanium  dioxide  span  the 
range  of  0.1  to  1.1  eV.99, 100, 217  221  The  techniques  used  to  make  these  measurements  included 
thermoluminescence  and  diffuse  reflectance  spectroscopy,  as  well  as  electrochemistry  and  purely 
theoretical,  computer  simulation.  M  y  0.3  ~  0.34  eV  range  of  values  from  infrared  spectroscopy 
compare  favorably  with  some  of  these  literature  values  that  are  presented  in  Table  5. 


Researcher  Shal lowTrap  State  Gap-s  (&J) 

Ghosh  et  al." 

Frindell  etal.217 

Lagemaatetal.100 

Szczepankiewicz  et  al.218 

Valentin  etal.219 

Warren  et  al.220 

Peter  etal.221 

Table  5:  Survey  of  reported  values  for  titanium  dioxide  shallow  trap  state  energy  gap,  relative  to 
bottom  edge  of  conduction  band,  (a)  indicates  a  calculated  value. 

The  placement  of  the  shallow  trap  state  energy  within  the  context  of  other  features  of  the 
titanium  dioxide  energy  band  diagram  is  shown  in  Figure  46.  As  previously  noted  in  T able  1  on 
page  30,  the  gap  separating  the  bottom  of  the  conduction  band  and  the  top  of  the  valence  band  is 
approximately  3.0  eV.  The  Fermi-D irac  expression  presented  in  Equation  (64)  that  relates  the 
integrated  infrared  baseline  area  to  conduction  band  electron  density  was  derived  with  an 
assumption  that  the  E  variable  in  the  exponential  term  represents  the  electron  trap  energy  closest 
to  the  conduction  band,  i.e.  a  shallow  trap  state.  The  energy  difference  between  the  conduction 
band  and  this  shallow  trap  state  is  also  indicated  in  Figure  46.  Other  features  noted  in  this  figure 


0.27,  0.28,  0.32,  0.39,  0.48,  0.56, 
0.60,  0.62,  0.76,  0.87 

0.1  ~0.6 

0.4  a 

0.42 

0.9,  1.1 

0.1 

0.1-0.25 


106 


are  the  vacuum  energy  level  energy  ( Evacuum )  for  titania,  i.e.  the  energy  at  which  an  electron  is 
completely  removed  from  the  semiconductor,  and  the  Fermi  energy  ( E Fermi )■ 


Egap  (~  3  eV) 


0.5  E, 


■gap 


vacuum  level 

CB 

conduction  band  bottom 

^  Epermi 

shallow  trap  state  — 

valence  band  top 

T 


0.3  ~  0.34  eV 


Figure  46:  The  shallow  trap  state  energy  difference  of  titanium  dioxide  relative  to  the  band  structure. 


Accurate  placement  of  the  Fermi  energy  for  a  semiconductor  nanoparticle  such  as 
titanium  dioxide  is  difficult.  Experimental  approaches  to  measuring  Fermi  energy  have  included 
electron  spectroscopy  techniques98  and  observation  of  electric  current- voltage  behavior  at  a 
junction  between  semiconductors  of  known  and  unknown  Fermi  energies.222  FI  ayes  and 
Stoneham  noted  that  for  a  pure,  undoped  semiconductor  the  Fermi  energy  is  located  in  the 
middle  of  the  band  gap,  but  when  dopants  are  introduced  the  Fermi  level  is  increased  above  the 
midpoint  of  the  gap.223  Gopel  et  al.  placed  the  Fermi  level  less  than  1  eV  below  the  conduction 
band.98  Given  that  electron  trap  state  populations  within  the  band  gap  are  variable,  location  of 
the  Fermi  energy  for  a  nanoparticle  without  direct  measurement  is  limited  to  estimates. 

4.4  E  stimation  of  E  lectron  Shallow  T raps  Population 

Titanium  cations  are  the  most  likely  locale  of  titanium  dioxide  nanoparticle  shallow 
electron  traps.  While  electron  micrographs  such  as  Figure  45  on  page  104,  Appendix  6  Figure 
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50  on  page  128,  and  Appendix  7  Figure  51  on  page  129  give  the  appearance  that  nanoparticles 
are  smooth  spheres,  at  the  atomic  scale  the  surface  of  the  particle  contains  numerous  steps  and 
corners  where  titanium  cations  are  incompletely  coordinated  with  oxygen  anions.  This  particle 
surface  atomic  'roughness'  is  discernable  in  the  high  resolution  transmission  electron  micrograph 
shown  in  Appendix  8  Figure  52.  Larger  particles  containing  multiple  titania  polymorph  phases 
may  have  grain  boundaries  that  also  can  contribute  irregularities  in  cation  coordination  leading  to 
electron  trapping.224  Coordinate  unsaturated  cations  at  the  particle  surface  are  particularly 
susceptible  to  capture  of  electrons  from  adsorbates,  with  a  reduction  of  oxidation  state  from  Ti4+ 
to  Ti3+.  Loss  of  oxygen  from  titanium  dioxide  with  an  accompanying  formation  of  T i3+  as  shown 
in  Figure  15  on  page  35  is  also  more  likely  to  occur  at  the  particle  surface.  Also,  Lagemaat  et  al. 
noted  that  lattice  defect  electron  trap  sites  are  presumed  to  be  located  on  the  surfaces  of 
nanoparticles  because  of  their  high  surf  ace- to- volume  ratio.100 

An  attempt  to  estimate  the  number  of  shallow  trap  sites  could  begin  with  the  Equation 
(64)  Fermi-Dirac  expression.  While  this  expression  gave  estimates  of  the  shallow  trap  energy 
level  that  were  in  reasonable  agreement  with  the  work  of  other  researchers,  as  shown  in  T able  5, 
closer  examination  of  Equation  (64)  does  not  offer  a  means  to  mathematically  estimate  the  total 
number  of  trap  sites  present  on  a  titania  nanoparticle  particle. 

Consider  the  example  of  a  4  nm  diameter  titania  particle  that  is  composed  entirely  of  the 
rutile  polymorph  (density  =  0.235  g  •  cm"3)84.  Evaluation  of  Equation  (61)  for  Estate  =  0.3  eV  at  a 
temperature  of  298  K  indicates  that  approximately  1.69  x  1016  electrons  are  excited  into  the 
conduction  band  for  one  gram  of  nanoparticles  (See  Appendix  9  on  page  131).  Implicit  in  this 
example  is  the  assumption  that  the  number  of  electrons  entering  the  conduction  band  equals  the 
number  of  electron-donor  shallow  trap  sites.  The  Fermi-Dirac  trends  in  Figure  44  demonstrate  a 
correlation  between  nanoparticle  temperature  and  FTIR  integrated  baseline  area.  This  increase 
in  integrated  FTIR  baseline  area  is  an  indication  of  increased  conduction  band  electron 
population.  If  evaluation  of  Equation  (61)  is  repeated  under  similar  constraints,  with  the 
nanoparticle  temperature  changed  to  675  K,  the  new  conduction  band  electron  population  is 
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approximately  1.51  x  1018  electrons  per  gram  of  rutile  nanoparticles..  If  the  number  of  electrons 
in  the  conduction  band  is  equal  or  proportional  to  the  number  of  shallow  trap  sites,  there  should 
exist  a  temperature  for  Equation  (61)  where  a  maximum  conduction  band  electron  population 
occurs.  The  number  of  conduction  band  electrons  at  this  temperature  could  serve  as  an  estimate 
of  the  number  of  electron  trap  states.  This  maxima  should  be  accessible  through  evaluation  of 
the  roots  of  a  first-order  derivative  of  Equation  (61). 

The  first-order  derivative  of  Equation  (61)  with  respect  to  the  independent  temperature 
variable  is: 
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where  B  and  C  are  aggregate  constants,  Estate  is  the  shallow  trap  state  energy,  and  T  is  the 
thermodynamic  temperature.  The  first  multiplied  term  of  Equation  (76)  and  the  right-most 
exponential  term  can  only  reach  a  zero  value  at  an  infinite  temperature,  as  shown  in  Equations 
(77)  and  (78).  The  remaining  middle  term  of  Equation  (76)  achieves  a  zero  value  at  a  negative 
thermodynamic  temperature  value,  which  has  no  physical  meaning,  as  shown  in  Equation  (80). 
These  results  exclude  a  purely  analytical  solution  to  estimation  of  shallow  trap  population  from 
the  Equation  (61)  Fermi-Dirac  expression. 
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(80) 


It  is  intuitive  that  a  nanoparticle  has  a  finite  number  of  shallow  state  trapped  electrons. 
An  experimental  determination  of  the  number  of  trap  states  from  recording  additional  infrared 
spectra  at  temperatures  above  675  K  may  not  succeed  at  determination  of  an  electron  thermal 
excitation  maximum.  Excessive  heating  of  titania  nanoparticles  will  lead  to  the  creation  of 
additional  oxygen  vacancies  with  creation  of  new  trap  sites.  A  nnealing  of  titanium  dioxide  will 
alter  the  polymorph  composition  of  mixed  phase  nanoparticles,  obscuring  the  trap  state 
population  of  the  starting  material. 

4.5  FT  I R  Surface  Chemistry  Observation  and  Electron  Monitoring 

Infrared  spectroscopy  is  an  indispensable  tool  for  the  study  of  nanoparticle  surfaces. 
FTIR  techniques  provide  a  fast,  nondestructive  method  for  identification  of  surface  adsorbates. 
Shifts  of  characteristic  infrared  absorbance  frequencies  also  provide  information  about  the 
surface  of  adsorption  in  addition  to  qualitative  identification  of  the  adsorbates.  For  example,  a 
red-shifted  infrared  absorbance  frequency  (i.e.  shifted  to  lower  value,  corresponding  to  a  longer 
wavelength)  is  an  indication  of  weakened  bonding  in  the  adsorbate.  This  condition  could  arise 
from  withdrawal  of  adsorbate  electron  density  through  bonding  to  a  surface  Lewis  acid  site. 
Alternately,  population  of  adsorbate  antibonding  molecular  orbitals  with  electrons  from  the 
adsorption  surface  could  weaken  the  adsorbate  bond  and  lead  to  a  red-shifted  absorbance 
frequency.  The  complementary  utility  of  using  the  same  infrared  spectroscopy  to  simultaneously 
monitor  electron  population  transitions  is  that  it  provides  information  to  discern  which  of  these 
cases  might  be  applicable. 

Infrared  monitoring  of  electron  population  offers  further  information,  beyond  correlation 
of  adsorbate  IR  frequency  shifts.  A  common  reaction  mechanism  step  for  catalytic  adsorbate 
oxidation  on  reducible  oxide  surfaces  (in  the  absence  of  gaseous  oxygen)  is  lattice  oxygen 
extraction  through  the  Mars  and  van  Krevelen  mechanism.225"229  For  titanium  dioxide,  this 
lattice  oxygen  extraction  is  accompanied  by  the  formation  of  additional,  populated  shallow 
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trapped  electrons  (see  Figure  15  on  page  35)  which  are  thermally  excited  to  the  conduction  band. 
Through  IR  monitoring  of  electron  populations,  a  lattice  oxygen  extraction  mechanism  could  be 
identified  for  adsorbate  oxidation.  Similar  scenarios  can  be  described  where  nanoparticle 
electron  population  can  be  used  to  monitor  surface  coverage  by  electron  withdrawing  adsorbates 
(e.g.  molecular  O2  or  H20)  or  their  displacement  from  the  particle  surface  by  reactive  adsorbates. 
46  Future  Work 

The  FTIR  vacuum  cell  offers  a  highly  controlled  environment  to  observe  nanoparticle 
surface  chemistry.  This  control  of  experimental  variables  comes  at  a  cost  of  time.  The 
procedures  necessary  to  place  new  samples  into  the  cell  and  condition  them  can  require  days  to 
weeks  before  the  cell  and  the  samples  within  it  are  ready  for  experimental  work.  In  addition  to 
slow  throughput,  observation  of  nanoparticle  catalysts  under  high  vacuum  is  not  representative 
of  typical  catalyst  applications  . 

Dr.  Dimitar  Panayotov  has  designed  and  had  fabricated  a  quartz  flow  cell  for  FTIR 
nanoparticle  observation  at  nearly  atmospheric  pressure.  This  spectrometer  cell  is  designed  to 
contain  a  square  tungsten  mesh,  identical  to  the  material  used  as  a  nanoparticle  carrier  in  the 
FTIR  vacuum  cell.  With  resistive  heating,  the  FTIR  flow  cell  is  capable  of  conducting 
nanoparticle  catalysis  studies  that  parallel  the  vacuum  cell  studies  with  a  continually  replenished 
supply  of  reactant. 

The  geometry  of  the  flow  cell  is  arranged  in  a  diamond  orientation,  with  gas  flow  into 
and  out  of  the  cell  occurring  at  opposite  corners  of  the  gas  cavity,  as  shown  in  Figure  47.  The 
gas  expansion  zone,  located  at  the  widest  part  of  the  gas  cavity,  falls  directly  above  the  location 
of  the  nanoparticle  sample  mounted  within  the  cell.  At  this  region  of  gas  expansion,  the  linear 
velocity  of  the  gas  decreases,  affording  greater  time  for  gas  -  nanoparticle  interaction. 


Ill 


Figure  47:  Nanopartide  FTIR  flow  cell.  The  photographic  inset  in  the  upper  left  corner  is  the 
prototype  flow  cell  and  flow  cell  mounting  bracket.  The  circular  region  in  the  photograph  is  a  sample 
of  Degussa  P25  titanium  dioxide  nanoparticles. 


The  FTIR  flow  cell  is  attached  to  an  external  gas  routing  system  that  allows  for  insertion 
of  gaseous  compounds  or  vaporized  liquids  into  the  carrier  gas  stream  flowing  through  the  cell. 
We  hope  to  use  this  apparatus  to  examine  chemical  warfare  agents  interaction  with  nanoparticle 
system  under  conditions  more  accurately  representative  of  applied  usage  of  the  catalysts. 

Another  advantage  of  the  gas  flow  cell  approach  is  quick  sample  through  put.  Because 
the  flow  cell  operates  at  approximately  atmospheric  pressure,  its  gas  seals  can  be  less  robust  than 
those  of  a  high  vacuum  chamber.  With  less  stringent  gas  seal  requirements,  the  flow  cell  can  be 
rapidly  set  up  or  disassembled.  It  is  hoped  than  when  the  flow  cell  apparatus  is  fully  operational, 
several  nanoparticle  samples  per  week  can  be  examined  as  opposed  to  one  or  two  sets  of 
nanoparticles  per  month  for  the  FTIR  vacuum  cell.  With  this  level  of  sample  through  put,  the 
flow  cell  can  serve  as  screening  tool  to  identify  nanoparticle  samples  that  show  promising 
characteristics  which  merit  more  detailed  examination  under  high  vacuum  conditions. 
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Other  interests  include  improvement  of  the  laser  vaporization  procedure  for  nanoparticle 
synthesis.  I  would  like  to  explore  application  of  filter  trapping  of  nanoparticles ,  electrostatic 
nanoparticle  capture,  and  size  segregation  of  particles  to  obtain  more  monodisperse  samples  of 
particles  for  particle  size  effect  studies. 
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Appendix  A:  Density  of  States  Derivation 


Blakemore230  presented  this  derivation  for  calculating  the  density  of  electron  states  within 
a  crystal  lattice.  Several  assumptions  are  implicit  in  this  derivation: 

•  M  otions  of  electrons  are  restricted  to  the  volume  of  the  crystal. 

•  There  are  a  finite  number  of  electron  states  for  the  specified  range  of  energies  and  momenta. 

•  The  crystal  lattice  is  much  larger  than  the  sizes  if  its  component  atoms. 

•  The  electron  'sees'  a  homogeneous  environment  (uniform  potential  field)  within  the  lattice. 

•  The  electron  potential  energy  \/ is  zero  within  the  lattice  and  infinite  outside  the  lattice. 

•  The  effects  of  the  crystal  surfaces  are  negligible. 

•  The  electron  wave  function ¥  is  periodic  along  each  Cartesian  axis  of  the  lattice. 

An  electron  delocalized  from  its  original  atom  within  a  crystal  I  attic  is  free  to  move 
within  the  confines  of  that  lattice  space.  Such  an  electron  has  momentum,  represented  by  the 
vector  p.  Quantum  mechanics  govern  the  actions  of  delocalized  electrons  and  electrons  possess 
wave-particle  duality.  Thede  Broglie  wavelength  of  an  electron  is: 

h  h 

A=-=—  (81) 

P  m,,v 

where  p  is  the  magnitude  of  the  vector  p,  me  is  the  effective  electron  mass,  and  h  is  Planck's 
constant.  Each  electron  can  also  be  represented  by  a  wave-vector: 

k  ——  where  fi=—  (82) 

h  In 

_  — ^ 

The  wave-vector  k  with  scalar  components  kx,  ky,  and  kz,  exists  in  a  Cartesian  coordinate  vector 

— ^ 

space  that  is  constructed  by  drawing  a  vector  k  from  the  origin  to  a  point  k.  The  spherical 

— ^ 

surface  formed  by  all  wave-vectors  k  for  a  given  electron  represents  the  degenerate  energy 

— ^ 

states  for  that  electron.  Figure  48  shows  /c-vector  space  for  an  electron: 
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Appendix  A:  Density  of  States  Derivation  (Continued) 


Figure 4&  k  vector  ^Dacefor  an  electron. 


For  an  electron  represented  by  a  k  wave-vector  the  following  conditions  must  be  true  for  the 
vector  components: 


p  =fik  p  =hk  p  =Tik 

1  x  x  •  y  y  •  z  z 


(83) 


The  kinetic  energy  of  the  electron  is: 


E  =  —(k2  +  k2  +  k2)=trk2 
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The  time  independent  Schrodinger  equation  for  the  electron  within  the  crystal  lattice  is: 


2m 


V2vF  =Ex¥ 


(85) 


An  acceptable  solution  for  the  wave  function  VF  is 


T  =  Cexp(z'A  •  r) 


(86) 
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Appendix  A:  Density  of  States  Derivation  (Continued) 


where  r  is  the  radial  displacement  vector  for  k.  The  wave  function  *F  is  defined  such  that  the 
probability  of  finding  an  electron  within  a  volume  element!  is  given  by  'l'  dr.  The  periodicity 
of  the  wave  function  'F  is  L.  Acceptable  solutions  for  the  wave  function  'F  correspond  to 
electron  states  within  the  lattice,  if  the  following  conditions  are  true: 


(87) 


where  nx,  ny,  and  nz  are  integers  and  can  assume  the  values  0,  ±1,  ±2,  ±3, ...  ,±n. 

The  Heisenberg  uncertainty  principle  requires  that  if  an  electron  position  is  determinate 
in  L  the  electron  must  be  indeterminate  in  momentum  and  wave-vector: 


h  'T.k 

A p  =—  and  A k  —  (88) 

L  L 

— ^ 

Only  specific  positions  in  Ar-space  correspond  to  acceptable  wave  functions  and  each  acceptable 
position  is  within  a  cell  of  volume  (27i/L)3.  Figure  49  shows  a  range  of  such  positions  in  the  x-y 
plane  and  an  example  of  a  cell: 


Figured  Acceptable  wave  function  positions  in  thex-y  plane 


The  number  of  electron  states  within  a  sphere  of  radius  A: is  found  by  dividing  the  volume  of  the 
sphere  by  the  volume  of  a  single  wave  function  cell: 
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Appendix  A:  Density  of  States  Derivation  (Continued) 


f  4  A 

r  l1) 

N  = 

13  J 

&3L3 
6n 2 


(89) 


The  change  in  the  number  of  states  for  a  change  in  sphere  radius  Aris: 
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The  change  in  the  density  of  electron  states  per  unit  volume  is: 


dN  k  „ 
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(91) 


Electrons  are  fermions  and  can  assume  two  possible  spin  states  of  +Vi  and  -Vi.  This  property 
doubles  the  density  change  of  allowable  electron  states: 
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The  density  of  states  function  g  is  defined  for  an  infinitesimal  change  in  energy  or  wave-vector 
as: 


g(E)dE  =g{k)dk  = 
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dk 


(93) 


The  vector  magnitude  /ccan  be  solved  for  in  Equation  (84)  to  give: 
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Appendix  A:  Density  of  States  Derivation  (Continued) 
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Appendix  B:  Fermi  -  Dirac  Distribution  Derivation 


Blakemore230  presented  this  derivation  for  the  Fermi-D irac  distribution  Law.  This 
derivation  requires  several  assumptions: 

•  A  very  large  number  of  electrons  exist. 

•  The  electrons  at  a  given  energy  are  indistinguishable. 

•  The  number  of  electrons  in  the  system  and  the  total  energy  of  all  electrons  is  fixed. 

•  The  system  is  in  thermal  equilibrium. 

For  g,  electron  states  with  Ei energies,  n,-  represents  the  number  of  occupied  electron  states  and  g 
-  rii  represents  the  number  of  unoccupied  states.  The  number  of  distinguishable  distributions  is: 


In  addition  to  the  g,  occupied  electron  states  at  energies  Ei,  there  are  other  states  at  higher  and 
lower  energies.  The  combination  of  the  1/14  electron  states  and  these  other  states  comprises  the 
system.  The  total  number  of  distinguishable  system  distributions  wis: 


w 


=UWr  =n 


Sr 


nKg-n) 


(100) 


The  transitions  of  electrons  between  states  are  subject  to  these  conditions: 


n 


(101) 


nE  =V  n  E 

r  r 
r 


(102) 


The  most  probable  distribution  of  electrons  in  states  at  thermal  equilibrium: 
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Appendix  B:  Fermi  -  Dirac  Distribution  Derivation  (Continued) 


8w 

8n 


=0 


(103) 


The  factorials  of  very  large  numbers  in  the  probability  Equations  (99)  and  (100)  can  be  estimated 
with  Stirling's  approximation: 


In  y!  =ylny-  y  for  y  »  1 


Equation  (100)  becomes: 


In  w  =  lnSr  -  nr  ln  nr  -  U, .  -  nr)ln{sr  ~  Hr  ) 

r 

Rewrite  conditions  (101)  and  (102)  as: 

n-Yjnr  =° 


(104) 


(105) 


(106) 


nE  -JnE  =0 

4-U  r  r 


(107) 


The  maximum  value  of  In  w,  subject  to  conditions  (106)  and  (107)  is  found  by  the  method  of 
Lagrange  undetermined  multipliers.  At  a  maximum  for  In  w. 


_8_ 

8n 


f  \ 

r 

lnw  +  a 

n~yLnr 

+  P 

nE  - 

-Yn  E 

4—4  r  r 

K  r  7 

r  )  _ 

0 


(108) 


Differentiation  of  each  term  in  (108)  equation  reference  goes  here  gives: 

5 


8n 


-ln  w  =ln 


f  ^ 

g  -n 

°  r  r 


ii 

\  r  J 


(109) 
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Appendix  B:  Fermi  -  Dirac  Distribution  Derivation  (Continued) 
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_8_ 

8n 


P 


nE 


Y  n  E 

r  r 


=~PEr 


(111) 


Equation  (108)  becomes: 


In 


g  -n 
n  , 

V  r  J 


a  +  f5Er 


(112) 


gr  1  +  exp  (a  +  [3Er) 


(113) 


For  any  of  the  states  gr  within  the  system,  the  probability  that  an  nr  particle  at  energy  Er  occupies 
that  state  is: 


/(£)  = 


1 

1  +  exp(a  +  f5Er ) 


(114) 


If  we  allow  a  =  Ef/I<bT  where  EF  is  defined  as  the  Fermi  energy  and  ks  is  the  Boltzmann 
constant,  the  energy  of  the  highest  occupied  state  at  T  =  OK  and  /?  =  1/(I<bT)- 


f(Er)  = 


1 


f 

1  +  exp 

v 


Er  -  Er. 

KE 


(115) 


Q.E.D. 
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Appendix  C:  FORTRAN  Program  WaveShave.f 


PROGRAM  WaveShave 

INTEGER  NumPoints,  N,  EOF1,  EOF2 
REAL  ReadFreq,  ReadAbs,  LowFreq,  HighFreq 
CHARACTER*^  IRFileln,  IRFileOut,  Spacer 
Spacer  = '  1 

*  Frequency  Filter  Values 

LowFreq  =  1.9E3 
HighFreq  =  2.6E  3 

OPEN  (Unit=l,  File-FILELIST.TXT1,  Status='0 L D ') 

*  Read  in  the  data  files  to  integrate 
EOF1 =  0 

DO  WHILE  (EOF1  .EQ.  0) 

READ  (1,  FMT  ='(A12)',  IOSTAT=EOFl)  IRFileln 

*  Check  for  EOF 

*  IF  (EOF1  .EQ.  -1)  THEN 

*  GO  TO  700 

*  END  IF 
PRINT  *, 11 
PRINT  *,  IRFileln 
PRINT*  '************ 1 

*  Open  the  input  file 

OPEN  (Unit=2,  Fi le=l RFilel n,  Status='OL D ') 

*  Open  the  output  file  as  a  TXT  type. 

IRFileOut  =  I  RFilel  n(l:10)  //  'XT' 

EOF2  =0 

OPEN  (Unit=3,  F i I e=l R F i I eO ut,  Status='N E W ') 

DO  WHILE  (EOF2  .EQ.  0) 

ReadFreq  =  0 
Read  A  bs  =  0 

500  READ  (2,  *,  IOSTAT=EOF2)  ReadFreq,  ReadAbs 

*  Check  for  EOF 

IF  (EOF2  .EQ.  -1)  THEN 
GO  TO  550 
END  IF 

*  Check  for  data  points  under  the  low  wavenumber 

IF  (ReadFreq  .LT.  LowFreq)  THEN 

GO  TO  500 
END  IF 

*  Check  for  data  points  under  the  high  wavenumber 

IF  (ReadFreq  .GT.  H ighFreq)  THEN 

GO  TO  500 
END  IF 

PRINT  *,  ReadFreq,  Spacer,  ReadAbs 
WRITE  (3,  *)  ReadFreq,  ReadAbs 
550  CLOSE  (2) 
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Appendix  C:  FORTRAN  Program  WaveShave.f  (Continued) 


CLOSE  (3) 
600  END  DO 
700  CLOSE  (1) 

CLOSE  (3) 

END 
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Appendix  D:  FORTRAN  Program T rapint2.f 
PROGRAM  T rapint2 

INTEGER  NumPoints,  N,  EOF  1,  EOF 2,  Countl,  Count2,  N umDeltaFreq 

REAL  Freq,  FreqO,  Freql,  DeltaFreq,  SumDeltaFreq,  Abs,  LeftAbs,  RightAbs,  IntArea 

REAL  ReadFreq,  Read  Abs 

CHARACTER*^  IRFile,  Spacer 

Countl  =  0 

Count2  =  0 

Spacer  = 1  1 

*  Open  the  results  output  file 

OPEN  (Unit=3,  File-INTRESULT.TXT1,  Status='NEW') 

OPEN  (Unit=l,  File-FILELIST.TXT',  Status='OLD') 

*  Read  in  the  data  files  to  integrate 
EOF1 =  0 

DO  WHILE  (EOF1  .EQ.  0) 

READ  (1,  FMT  =  '(A  12)',  IOSTAT=EOFl)  IRFile 
PRINT 

PRINT  *,  IRFile 
PRINT  *,  '************' 

Countl  =  Countl  + 1 

*  Count  the  number  of  data  points  in  file 

EOF2  =0 
Count2  =  0 

OPEN  (U nit=2,  F i I e=l R F i I e,  Status='OL D ') 

DO  WHILE  (EOF2  .EQ.  0) 

Count2  =Count2  + 1 
ReadFreq  =  0 
Read  A  bs  =  0 

READ  (2,  *,  IOSTAT  =E 0 F 2)  ReadFreq,  ReadAbs 
IF  (EOF2  .EQ.  -1)  THEN 
GO  TO  1000 

ENDIF 
END  DO 

1000  NumPoints  =  Count2 
CLOSE  (2) 

*  Compute  average  DeltaFreq 
EOF2  =0 

Count2  =  0 
FreqO  =  0 
Freql  =  1 
DeltaFreq  =  0 
SumDeltaFreq  =  0 
N  umDeltaFreq  =  0 

OPEN  (Unit=2,  File=l RFile,  Status='OL D ') 

DO  WHILE  (EOF2  .EQ.  0) 
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Appendix  D:  FORTRAN  Program T rapi nt2.f  (Continued) 

1100  ReadFreq=0 

Read  A  bs  =  0 

READ  (2,*,  1 0  ST  AT  =E  0  F  2)  ReadFreq,  ReadAbs 
Count2  =Count2  + 1 

*  First  wavenumber 

IF  (Count2  .EQ.  1)  THEN 
Freql  =  ReadFreq 
GO  TO  1100 

END  IF 

*  Last  wavenumber 

IF  (Count2  .EQ.  N umPoints)  THEN 
GO  TO  1500 

END  IF 
FreqO  =  Freql 
Freql  =  ReadFreq 
DeltaFreq  =  Freql  -  FreqO 
NumDeltaFreq  =  N  umDeltaFreq  +  1 
SumDeltaFreq  =  SumDeltaFreq  +  DeltaFreq 
END  DO 

1500  DeltaFreq  =  SumDeltaFreq  /  NumDeltaFreq 
CLOSE  (2) 

*  Integrate  the  data  file 
EOF2  =0 

Count2  =  0 
I  ntA  rea  =  0 

OPEN  (U nit=2,  File=l RFile,  Status='OL D ') 

DO  WHILE  (EOF2  .EQ.  0) 

1550  Count2  =  Count2  + 1 

ReadFreq  =  0 
Read  A  bs  =  0 

READ  (2,  *,  IOSTAT  =EOF2)  ReadFreq,  ReadAbs 
IF  (Count2  .EQ.  1)  THEN 

I  ntA  rea  =  I  ntA  rea  +  R  ead A  bs 
GO  TO  1550 
END  IF 

IF  (Count2  .EQ.  NumPoints)  THEN 
I  ntA  rea  =  I  ntA  rea  +  R  ead  A  bs 
GO  TO  1600 
END  IF 

I  ntA  rea  =  I  ntA  rea  +  (2*  Read  A  bs) 

PRINT  '(15.4,  A2,  E15.7,  A2,  E15.7,  A2,  E25.10,  A2,  E  15.7)',  Count2, 
Spacer,  DeltaFreq,  Spacer,  I  ntA  rea.  Spacer,  ReadFreq,  Spacer,  ReadAbs 
END  DO 

*  Compute  the  trapezoidal  integral 
1600  I  ntA  rea  =  (I  ntA  rea/2 )* DeltaFreq 

PRINT  *,  =========================== 
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Appendix  D:  FORTRAN  Program T rapi nt2.f  (Continued) 


PRINT  '(E 20.10)',  IntArea 

WRITE  (3,  '(A  12,  E 20. 10) ')  I RFile,  I ntA rea 

CLOSE  (2) 

IF  (EOF1  ,EQ.  -1)  THEN 
GO  TO  2000 

END  IF 

1800  END  DO 
2000  CLOSE  (1) 

CLOSE  (3) 

END 
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Appendix  E:  Nanoparticle  Size  M  easurement  Data  Summary 


Nanoparticles 

Pellet  to 
Collector 
(cm) 

o2 

Press. 

(Torr) 

Number  of 
Measurements 

Mean 

Diam. 

(nm) 

Standard 

Deviation 

(nm) 

SPB013 

4.4 

1 

235 

1 

1 

SPB016 

4.4 

1 

105 

5 

1 

SPB017 

4.4 

1 

172 

4.8 

0.9 

SPB018  a 

4.4 

1 

88 

4.0 

0.6 

SPB024 

4.4 

5 

111 

5 

1 

SPB029 

4.4 

5 

97 

6 

1 

SPB022 

4.4 

10 

97 

6 

1 

SPB023 

4.4 

55 

120 

12 

3 

SPB025 

4.4 

90 

62 

12 

3 

SPB011 

4.4 

100 

334 

20 

5 

SPB019 

4.4 

200 

147 

22 

7 

SPB026 

4.4 

200 

156 

14 

4 

SPB037 

4.4 

200 

296 

27 

9 

SPB040 

4.4 

200 

226 

20 

6 

SPB027 

4.4 

310 

49 

24 

6 

SPB028 

4.4 

410 

183 

30 

10 

SPB020  b 

4.4 

500 

192 

30 

10 

SPB014 

6.3 

37 

369 

13 

3 

SPB015 

2.3 

37 

247 

12 

4 

a  Used  for  Fermi  - 
b  U  sed  for  Fermi  - 

Dirac  regressions,  vacuum  cell  sample  SI 

Dirac  regressions,  vacuum  cell  sample  SI 

Table  6:  TI02  nanoparticle  size  measurement  data. 
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Appendix  F:  SPB018  T ransmission  Electron  M  icrograph 


spb.tem.32.tif 

Ti02  Nanoparticles  100  nm 

Print  Mag:  290000x  @  8.  in  Direct  Mag:  160000x 

15:28  09/25/07 


Figure  50:  TEM  micrograph  for  SPB018  4.0  ±  0.6  nm  titanium  dioxide  nanopartides.  Inset  is  a 
histogram  of  manual  particle  diameter  measurements.  Labels  on  particles  are  diameter  measurement 
flags. 


Appendix  G:  SPB 020 Transmission  Electron  M  icrograph 


spb.tem.45.tif 

Ti02  Nanoparticles  100  nm 

Print  Mag:  290000x  @  8.  in  Direct  Mag:  160000x 

15:46  09/25/07 


Figure  51:  TEM  micrograph  for  SPB020  30  ±  10  nm  titanium  dioxide  nanoparticles.  Inset  is  a 
histogram  of  manual  particle  diameter  measurements.  Labels  on  particles  are  diameter  measurement 
flags. 
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Appendix  H:  SPB018  High  Resolution  Transmission  Electron  M  icrograph 


Figure  52:  High  resolution  micrograph  of  SPB018  titanium  dioxide  nanoparticles. 
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Appendix  I:  Conduction  Band  Electron  Population  Calculation 


This  appendix  outlines  the  details  of  the  calculation  for  the  number  of  electrons  excited 
into  the  conduction  band  for  a  4  nm  (4  x  10"9  m)  diameter  titanium  dioxide  nanoparticle 
composed  of  the  rutile  polymorph  (density  =  0.235  g  •  cm-3)84  with  Estate  =  0.3  eV.  The  density 
of  conduction  band  electrons  is  determined  from  Equation  (61): 


NrJE)=BTmExp 


CB 


-CE 


State 


>5/2 


where  B  = 


[mnk\  i 

V  e  *4  md  c  =  - 


h 


2kr 


The  aggregate  constants  B  and  C  are  evaluated: 

25/2(9.1094x  10  31^)(3.1416)(l.3807 x  10  23i  •  A'-1) 
(6.6261  xlO"34  7  »s)3 

1 J  =  IN*  m  =  I  ( kg  •  in  •  s  2  j  •  m  =  I  kg  •  m 2  •  ,v  2 


B 


4.8296  x  1021  K  3  •  m  3 


1  m 


v I OOcm j 


B  -  4.8296  x  1015  K  3  •  cm  3 


2(l. 3807  x  10^23^  •  m2  •  .r2  •  K _1) 
leV  =1.6022  x  10“19/ 


1 

^1.6022x10  19  7 " 

2(l.3807  xlO-23  7  •A”1) 

l  W  ) 

C  =  5.8021x  103  A  •  eV1 


For  temperature  T  =298  K,  Equation  (61)  is  evaluated: 
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Appendix  I:  Conduction  Band  Electron  Population  Calculation  (Continued) 


Ncb(E)  =  4.8296  x  lOl5K~>  •  cm~3  (298  Exp 


(5.8021xl03Z*ey1)(o.3ey) 
298  K 


N  cb{e)  =  1  22xW6  an 


The  number  of  conduction  band  electrons  for  one  gram  of  rutile  atT  =  298  K : 


electrons 

gram 


f 

(7.22xlOI6cm  3) 

V 


1  an:' 

426^, 


1.69x10 ‘V1 


ForT  =  695  K,  the  number  of  conduction  band  electrons  for  one  gram  of  rutile: 


Ncb(E)  =  4.8296  x  1015  •  cm^\(ei5K)m  Exp 


(5.8021xl03Z*ey1)(o.3ey) 
67 5  K 


7 


JVct(e)  =  6.43x10 


18  -3 

cm 


electrons 

gram 


(6.43xlOI8cm-3) 


^  lcm3  ^ 


1.51X1018,?  1 


Constants  and  Conversion  Factors?0 


Boltzmann  Constant: 

k  =  1.3807  x  10"23J  *K  1 

Plank  Constant: 

h  =  6,6261  xlO'34J-s 

Electron  Rest  M  ass: 

me=  9.1094  xlO'31  kg 

eV  toj  Conversion: 

1  eV  =  1.6022  x  10  19  J 
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